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COMPUTER  CODE  FOR  THE  ANALYSIS  OF  MULTILAYERED 
FIBER  COMPOSITES  -  USERS  MANUAL 
by  Christos  C.  Chamis 
Lewis  Research  Center 


SUMMARY 

A  computer  code  to  carry  out  the  multilevel  linear  analysis  required  to  efficiently 
design  structural  components  made  from  multilayered  fiber  composites  is  described. 

The  inputs  to  the  code  are  constituent  materials  properties,  some  factors  reflecting  the 
fabrication  process,  and  composite  geometry.  The  code  performs  the  micromechanics, 
macromechanics,  and  laminate  analysis  of  fiber  composites.  The  code  outputs  are  the 
various  ply  and  composite  properties,  composite  structural  response  (accounting  for 
bending-stretching  coupling  etc.),  and  the  composite  stress  analysis  results,  including 
the  results  of  the  combined-stress  strength  criteria.  The  code  is  in  FORTRAN  IV  com¬ 
piler  language  and  can  be  used  efficiently  as  a  package  in  complex -structural  analyses, 
finite -element  methods,  buckling  and  vibration  studies,  and  structural  syntheses.  The 
input-output  format  is  described  extensively.  Required  input  data  to  the  code  for  various 
fiber-matrix  composites  are  given.  The  FORTRAN  compiled  listing  and  sample  trial 
cases  are  included  to  aid  the  designer  or  analyst  in  using  this  code.  The  code  consists 
of  two  parts.  In  the  first  part,  the  mechanics  to  use  the  code  are  described;  in  the  second 
part  the  equations  programmed  are  described.  The  code  has  been  used  successfully  in 
the  analysis  of  various  fiber  matrix  multilayered  composites.  It  was  also  used  (and 
proved  to  be  efficient)  in  the  structural  synthesis  of  multilayered  thornel/epoxy  compos¬ 
ite  plates,  in  buckling  studies  of  simply  supported  multilayered  fiber-composite  plates, 
and  in  the  computation  of  lamination  residual  stresses  in  angle  ply  composites.  Selection 
of  correlation  coefficients  for  new  composite  systems  is  described.  Possible  extensions 
for  temperature -dependent  properties,  material  nonlinearities  and  failure  load  envelopes 

are  indicated. 


*A  part  of  the  work  described  herein  was  supported  by  the  Advanced  Research  Projects  Agency, 
Department  of  Defense,  through  contract  AF  33(615)- 3110  administered  by  the  Air  Force  Materials 
Laboratory  while  the  author  was  a  member  of  the  Engineering  Design  Center,  Case  Western  Re¬ 
serve  University,  Cleveland,  Ohio. 


INTRODUCTION 


The  importance  and  need  of  a  multilevel  analysis  in  designing  structural  components 
with  multilayered  fiber  composites  are  documented  in  reference  1.  A  multilevel  analysis 
which  was  found  to  be  efficient  in  predicting  the  structural  response  of  multilayered  fiber- 
composites  (with  the  constituent  materials  properties,  the  fabrication  process,  and  the 
composite  geometry  known)  is  also  documented  in  reference  1. 

The  multilevel  analysis  presented  in  reference  1  consists  of  (1)  micromechanics 
theories  for  the  thermoelastic  properties  and  the  stress -level  limit  of  the  single  ply  as 
functions  of  constituent  materia,ls  properties  and  the  particular  fabrication  process, 

(2)  the  combined-stress  strength  criterion  of  the  single  ply,  and  (3)  multilayered  com¬ 
posite  structural  response  and  analysis  (macromechanics  or  laminate  analyses)  where 
the  interply  layer  effects  are  taken  into  account.  The  computer  code,  to  carry  out  this 
multilevel  analysis  and  supplemented  as  noted  by  the  additional  references,  is  described 
herein. 

The  computer  code  has  been  programmed  in  FORTRAN  IV  and  has  been  extensively 
used  in  the  UNIVAC  1107,  1108,  and  IBM  7094.  Since  this  report  is  to  serve  as  a  user’s 
manual,  the  code  is  divided  into  two  parts.  In  the  first  part,  the  mechanics  of  using  the 
code  are  described  with  respect  to  program  format,  input -output,  sample  input  data 
sheets,  and  tables  of  input  data  for  several  fiber  composites.  In  the  second  part,  the  pro¬ 
gram  is  described.  Sample  case  runs  of  Thornel-50/epoxy  composites  with  unidirectional 
and  angle  plies  are  included  with  the  compiled  listing  in  appendix  B.  Sample  cases  for 
bending,  stretching,  coupling,  and  lamination  residual  stresses  are  also  included. 

The  format  of  the  program  is  described  in  the  section  MAIN  PROGRAIVI  and  follows 
the  FORTRAN  IV  program  format  for  the  7094.  The  subroutines  required  to  carry  out 
the  various  levels  of  analysis  are  described  individually  in  their  respective  sections.  In 
these  subroutines,  the  equations  programmed  are  given,  the  various  alternatives  for  es¬ 
tablishing  certain  properties  (such  as  strain  magnification  factors,  longitudinal  com¬ 
pressive  stress  Limit,  and  combined -stress  strength  criterion)  are  discussed,  and  the 
subroutine  input -outputs  and  the  global  storage  locations  (common  to  all  parts  of  the  pro¬ 
gram)  are  identified.  The  input-output  format  is  described  in  detail  separately  in  the 
sections  Input  Ply  Properties  and  Output.  These  descriptions  are  quite  extensive  so  that 
designers  and  analysts  with  little  or  no  programming  experience  as  v/ell  as  experienced 
programmers  can  easily  use  the  code. 

In  appendix  A  the  FORTRAN  symbols  are  defined.  The  definitions  include  such  in¬ 
formation  as  in  which  part  of  the  program  each  global  variable  is  generated.  The  input 
data  in  tables  IV  to  XllI  (currently  acceptable  for  the  analysis  of  several  fiber  compos¬ 
ites)  provide  for  immediate  use  of  the  code.  The  inclusion  of  the  compiled  FORTRAI4 
Isting  (appendix  B)  with  the  sample  trial  cases  (appendix  C)  should  further  amplify  the 
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detail  descriptions.  Filament  and  fiber  are  used  interchangeably  in  the  description  and 
in  the  discussion.  Ply,  unidirectional  laminate,  and  unidirectional  composite  are  also 
used  interchangeably. 

It  is  noted  that  the  global  storage  of  the  composite  and  ply  properties  is  very  impor¬ 
tant  when  this  multilevel  analysis  is  to  be  used  as  a  subroutine  package  to  generate  struc¬ 
tural  behavior  properties  for  structural  analysis  purposes. 

This  code  has  been  used  successfully  in  predicting  the  ply  thermoelastic  constants 
(refs.  1  to  3)  in  laminate  analysis  (boron,  graphite,  carbon,  and  glass -filament-epoxy 
composites;  refs.  1,  and  3  to  5)  buckling  analysis  (ref.  6)  and  structural  synthesis 
(ref.  1).  It  has  also  been  used  to  calculate  the  lamination  residual  stresses  in  angle -ply 
composites  (ref.  7). 

Mr.  Tom  Delivuk,  then  with  the  same  center,  converted  the  original  ALGOL  CODE 
to  the  initial  FORTRAN  IV  CODE  which  resulted  into  the  CODE  described  herein. 


SYMBOLS 


cx 

R 


cx 

BIDE 

^cx 

Cel 

Ce2 

CSANB 
D 


D 


cx 

R 

cx 


D. 


EpEcf 

^m’  ^cm 


E 


'fll,  etc. 
Ill,  etc, 


composite  axial  stiffness 
reduced  axial  stiffness 

Boolean -true  if  interply  effects  are  included 

composite  coupling  stiffness 

string  with  force  variables 

string  with  displacement  variables 

Boolean-true  if  membrane  and  axial  symmetry  exists 

composite  flexural  rigidities 

reduced  bending  rigidities 
displacement  vector 
filament  equivalent  diameter 
filament  elastic  constants 
ply  elastic  constants 
matrix  elastic  constants 
fiber  normal  modulus 
ply  normal  modulus 


mil,  etc. 


G, 


fl2,etc. 


G-, 


Zl2, etc. 


ml2,  etc. 


H, 


H 


kc 


c 

K 


K 


K 


cll, c22, c33 
cxy, cyy, cxy 
ill 


K 


Zll 


K 


mil 


m 


m 

^z 


sc 


M 

M 

m 

N 


cx 


c  ATx 


cx 


N 


c  ATx 


N. 


N. 


N 


Zc 


N 


pc 


matrix  normal  modulus 

fiber  shear  modulus 

ply  shear  modulus 

matrix  shear  modulus 

interply  distortion  energy  coefficient 

array  of  constituents  heai  conductivities 

composite  heat  capacity 

index,  generally  ply  or  interply 

composite  three-dimensional  heat  conductivities 

composite  two-dimensional  heat  conductivities 

fiber  heat  conductivity 

ply  heat  conductivity 

matrix  heat  conductivity 

apparent  filament  volume  ratio 

apparent  matrix  volume  ratio 

apparent  void  volume  ratio 

actual  fiber  volume  ratio 

actual  matrix  volume  ratios 

ply  apparent  fiber  volume  ratio 

ply  apparent  void  volume  ratio 

array  of  limiting  conditions 

applied  moment 

thermal  moments 

load  condition  index 

applied  membrane  loads 

thermal  force 

number  of  fibers  per  end 

number  of  plies 

number  of  load  conditions 

string  PROPC  length 
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Pc 

Pz 

Pep 

P^P 

%,i,P,  r,s 
R 

RINDV 

SmTetc. 

TLINP 

^cb 

x,y,z 

a 

c 


^csx 

H 

^cs 

^IV  ^Ic 
12,  etc. 
^Z12,etc. 
12,  etc. 


string  PROP  length 
composite  properties  array 
ply  properties  array 
string  PROPC 
string  PROP  main  program 
indices  to  print  out  string  PROP 
transformation  matrix 

Boolean-TRUE  true  if  displacements  are  known 
ply  limit  stresses 
ply  thickness 

FALSE  if  ply  thickness  is  calculated  internally 
composite  local  curvature  changes 
structural  reference  axes 
composite  coefficient  of  thermal  expansion 
fiber  thermal  coefficient  of  expansion 
ply  thermal  coefficients  of  expansion 
matrix  thermal  coefficient  of  expansion 
correlation  factors  for  ply  thermoelastic  properties 
correlation  factors  for  ply  heat  conductivities 
correlation  factors  for  ply  strength 

matrix  strain  magnification  due  to  ply  strain  in  the  presence  of  voids 
interply  layer  thickness 
reference  plane  membrane  strains 
ply  strains 

angle  between  composite  material  and  structural  axes 

angle  between  ply  material  and  composite  axes 

fiber  Poisson's  ratio 

ply  Poisson's  ratio 

matrix  Poisson's  ratio 

constant 
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Pj  ^  filament  and  matrix  weight  density 

0’^  ply  stresses 

ma.trix  strain  magnification  due  to  ply  strain 
1,2,3  material  reference  axes 


USERS  MANUAL 


The  mechanics  required  to  use  this  code  for  the  analysis  of  multilayered  fiber  com¬ 
posites  are  described  in  this  part  of  the  manual.  Here,  it  is  assumed  that  the  user  is 
interested  in  using  the  code  as  a  tool  only  and  that  he  has  available  to  him  a  FORTRAN  IV 
manual.  The  theory  on  which  the  code  is  based  is  described  in  the  second  part  of  the 
report. 

The  physical  representation  of  the  code  is  illustrated  in  figure  1.  The  geometry  of 
the  constituents,  the  ply,  and  the  composite  are  defined  in  this  figure.  The  required  in¬ 
put  properties,  correlation  coefficients,  and  computed  properties  are  summarized  in  fig¬ 
ure  1  in  symbolic  form. 

The  physical  arrangement  of  the  code  is  illustrated  in  figure  2 .  The  numbers  given 
in  each  block  of  cards  are  for  subsequent  discussion  and  do  not  appear  on  the  code.  Four 
steps  are  required  to  use  the  code  in  the  user’s  computer  facility: 

(1)  Obtain  the  code. 

(2)  Make  it  operational  in  the  user's  computer  facility. 

(3)  Supply  the  input  data. 

(4)  Interpret  the  code  output  results. 


Obtain  the  Code 

The  code  could  be  obtained  in  cards.  If  this  is  not  convenient  or  possible,  then  the 
cards  can  be  punched  from  the  compiled  listing  (see  appendix  B) . 


Make  It  Operational 

Making  the  program  operational  requires  the  availability  of  a  FORTRAN  compiler 
in  the  user's  computer  facility,  certain  control  cards  at  the  beginning  of  the  code,  and  the 
card  that  preceeds  each  subroutine.  Consult  your  computer  group  about  these  items.  The 
control  cards  present  in  the  code  are  only  for  the  Lewis  IBM  7044/7094  direct  couple 
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system.  Once  the  deck  of  cards  has  been  assembled  as  is  shown  in  figure  2  (except  Input 
Data)  with  the  proper  control  cards,  the  user  is  ready  to  compile  the  code  in  his  facility. 
The  compilation  will  indicate  whether  any  additional  modifications  are  needed.  Most 
modifications  will  be  minor  and  will  usually  deal  with  certain  logical  statements  peculiar 
to  each  compiler.  Consult  your  computer  group  for  these  modifications. 


Supply  the  Input  Data 

The  physical  arrangement  of  the  input  data  cards  is  illustrated  in  figure  3.  The  num¬ 
bers  in  the  group  of  cards  are  for  identification  purposes  in  this  description  and  do  not 
appear  on  the  cards.  Details  in  preparing  the  input  data  cards  are  summarized  in  table  I. 
A  detailed  description  of  these  cards  is  given  subsequently.  A  sample  for  preparing  input 
data  sheets  is  illustrated  in  table  n  for  the  Thornel~50/epoxy  composite  system. 

Listings  of  input  data  for  several  composite  systems  appear  in  tables  HI  to  XII. 

These  systems  are  shown  graphically  in  figure  4.  The  input  data  for  these  systems  can 
be  punched  from  the  listings,  and  the  cards  that  need  alterations  for  the  specific  problem 
can  be  modified  accordingly. 

Input  data  for  additional  composite  systems  may  be  easily  prepared.  This  is  done  by 
selecting  a  related  system  from  those  in  tables  IH  to  XH  and  modifying  those  entries  that 
need  modification.  Table  I  and  the  section  Detail  Description  of  Input  Data  explain  where 
and  how  each  entry  is  read  in. 

After  the  input  data  have  been  properly  assembled  (as  is  shown  in  fig.  3) ,  it  is  placed 
in  its  physical  position  (fig.  2),  and  the  code  is  ready  to  be  run  for  results. 


Detailed  Description  of  Input  Data 

The  card  group  numbers  referred  to  here  are  given  in  figure  3  and  table  I.  The  se¬ 
quential  order  of  the  entries  in  each  card  group  is  given  in  table  I. 

(1)  Composite  system  card.  -  The  composite  system  title  is  punched  on  this  card. 

The  title  can  be  55  characters  long  including  blanks. 

(2)  Data  control  card.  -  The  number  of  plies  N^,  number  of  ply  properties  N  ,,  num¬ 
ber  of  composite  properties  the  number  of  fibers  per  end  Nf,  and  the  number  of 

load  conditions  Np^  are  entered  in  this  card.  The  number  of  ply  properties  and  the 
number  of  composite  properties  are  always  the  same:  they  are  Np^  =  71  and  N  ^  =  54- 
The  others  have  to  be  entered  according  to  the  composite  system  and  the  load  conditions. 

(^)  Constituent  materials  elastic  properties.  -  The  constituent  elastic  properties  are 
entered  in  this  group  of  cards.  The  fiber  properties  are  entered  first  and  then  the  ma- 
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rix.  Enter  only  extension  moduli,  Poisson's  ratios,  and  zero  values  for  shear  moduli 
when  the  constituent  material  is  isotropic.  For  example,  in  a  glass/ epoxy  system, 

Ef33  =  ®f22  =  ^fll  "  ''flS  =  "^£12-  ^f23  =  ^fl3  =  ^fl2 

computed  internally. 

(4)  Correlation  coefficients  for  ply  elastic  constants,  expansion  coefficients,  and 
strain  magnification  factors.  -  The  correlation  coefficients  that  make  theory  agree  with 
experiment  are  entered  in  this  group  of  cards.  The  first  entry  in  this  group  is  It 

is  selected  so  that  predicted  extensional  moduli  and  Poisson’s  ratios  correlate  with 


measured  values.  The  procedure  for  selecting  /3^  is  iterative.  First  the  code  is  run 
with  jS,  equal  to  some  initial  vaTue.  Experience  has  proven  that  =  4  is  usually  a 
good  approximation  for  the  initial  value.  For  many  systems  this  is  also  the  terminal 
value.  Next,  obtain  values  for  i3^  greater  and  smaller  than  four,  and  select  the  proper 
value  for  |3^  by  interpolation.  The  aforementioned  description  for  selecting  applies 
to  the  selection  of  all  correlation  coefficients  in  this  code. 

The  second  entry  in  this  card  group  is  /3^,  v/hich  is  the  correlation  coefficient  for 
the  ply  shear  moduli  %3-  Thejhird  entry  is  /3”  which  is  the  correlation 

coefficient  for  G^g.  The  fourth  entry  is  /3^  which  is  the  correlation  coefficient  for  the 
ply  thermal  coefficients  of  expansion.  The  next  three  entries,  P'^,  a,nd  P”,  are  cor¬ 
relation  coefficients  for  strain  magnification  factors  ^'^/il2’  ®p23’  ^®spec- 

tively.  These  coefficients  are  entered  as  zeros.  Experience  v;ith  several  composite 
systems  ha,s  shown  that  the  correlation  coefficients  P^  are  not  needed.  However,  they 
are  provided  for  possible  future  use. 

The  coefficient  P^  is  the  ratio  of  the  thickness -to -width  of  the  rectangle  formed  by 
an  in -situ  end  or  tow  of  fibers.  Another  way  to  visualize  this  is  that  P^  is  the  ratio  of 
ply -thickness  per  end  or  tow -ply  width.  The  value  for  P^  is  obtained  from  electron 
photomicrographs  or  indirectly  as  described  in  reference  5.  Entries  9  and  10  are  enter¬ 
ed  as  zeros;  these  fields  are  empty  and  are  available  for  future  use.  Entries  11  to  14 
are  for  the  coefficients  these  coefficients  are  alternates  to  and  a.re  to  be  used 

if  the  P  coefficients  do  not  provide  the  desired  correlation.  Note  that  when  a  |3^  co- 
efficient  is  used,  the  corresponding  coefficient  is  entered  as  zero  and  vice  versa. 

Entries  15  to  17  are  for  the  coefficients  y^,  which  are  alternates  to  coefficients  P^. 
Note  that  when  a  P^  coefficient  is  used  the  corresponding  y^  coefficient  is  zero  and 
vice  versa.  The  y  coefficients  are  entered  with  values  of  one.  Entries  18  to  20  are 
entered  as  zeros.  These  are  empty  fields  and  are  available  for  future  use. 

Experience  with  the  code  thus  far  has  shown  that  all  the  correlation  coefficients  ex¬ 
cept  P^  are  approximately  the  same  for  several  composite  systems.  (See  tables  III  to 

xn.) 


a 


(5)  Fiber  thermal  coefficients  of  expansion.  -  The  coefficients  01^22’ 

jgg  are  entered  on  this  card.  If  the  fiber  is  isotropic,  then  =  ®f22  ~  *^fll' 
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(6)  Matrix  thermal  coefficients  of  expansion.  -  The  coefficients  ^m22’ 

^m33  entered  in  this  card.  When  the  matrix  is  isotropic,  =  ^xn22  ~  ‘^mlT 

(7)  Constituent  heat  conductivities  and  heat  capacities.  -  The  first  four  entries  in 

this  group  are  for  the  fiber  heat  conductivities  K^22’  ^f33’  capacity  h^^.. 

The  next  four  are  for  the  corresponding  matrix  properties.  The  next  three  are  zero 
entries,  and  the  last  one  is  the  heat  conductivity  for  air.  (See  card  group  7  of 
table  I. ) 

(8)  Correlation  coefficients  for  heat  conductivities.  -  The  four  entries  in  this  card 

are  for  the  correlation  coefficients  ^h3’  ^®sP®ctively.  These  co¬ 
efficients  are  as  follows:  |3j^  is  for  matrix  with  voids,  for  ^l22’ 

and  /Sj^g  for  K^33-  They  are  selected  as  was  described  in  /3^  in  card  group  (4). 

(9)  Constant  n.  -  The  value  for  ti  is  entered  in  this  card. 

(10)  Boolean  for  thickness.  -  The  letter  T  is  entered  in  this  card  if  the  ply  thickness 
is  supplied.  The  letter  F  is  entered  if  the  ply  thickness  is  computed  internally. 

(11)  Boolean  for  membrane  and  bending  symmetry.  -  The  letter  T  is  entered  in  this 
card  if  the  composite  has  both  membrane  and  bending  symmetry;  otherwise  the  letter  F 
is  entered. 

(12)  Boolean  for  interply  layer  contribution.  -  The  letter  T  is  entered  in  the  card  if 
the  interply  layer  contributions  on  the  composite  are  desired;  otherwise,  the  letter  F  is 
entered. 

(13)  Boolean  for  input  displacements.  -  The  letter  T  is  entered  in  this  card  if  the 
displacements  are  inputs;  otherwise,  the  letter  F  is  entered. 

(14)  Composite  angle,  constituent  densities,  and  fiber  eqviivalent  diameter.  -  The 
composite  angle  (angle  between  composite  material  1-axis  and  structural  x-axis  (fig.  5) 
is  the  first  entry  in  this  card.  The  fiber  and  matrix  densities  are  the  second  and  third 
entries.  The  fourth  entry  is  the  fiber  equivalent  diameter. 

(15)  Ply  void  volume  ratio.  -  The  void  volume  ratio  of  the  plies  is  entered  in  this 
group  of  cards;  the  first  entry  is  for  the  first  ply,  and  the  last  entry  is  for  the  last  ply. 
The  bottom  or  the  inner  ply  in  the  composite  is  selected  as  the  first  ply  for  convenience. 
The  number  of  entries  is  equal  to  the  number  of  plies  in  the  composite.  (See  tables  I 
and  II. ) 

(16)  Ply  fiber  volume  ratio.  -  The  ply  fiber  volume  ratio  is  entered  in  this  group  of 
cards.  The  first  entry  is  for  the  first  ply,  which  is  the  bottom  or  inner  ply  in  the  com¬ 


posite.  The  last  entry  is  for  the  last  ply.  The  number  of  entries  equals  the  number  of 
plies.  (See  tables  I  and  II. ) 

(17)  Ply  orientation  angle.  -  The  ply  angle  (measured  from  the  composite  material 
1-axis  to  the  ply  material  1-axis  (fig.  5))  is  entered  in  these  cards.  The  first  entry  is 
for  the  first  ply  which  is  the  bottom  or  inner  ply  in  the  composite.  The  last  entry  is  for 
the  last  ply.  The  number  of  entries  equals  the  number  of  plies.  (See  tables  I  and  II. ) 
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(18)  Ply  thickness.  -  The  ply  thicknesses  are  entered  in  this  group  of  cards.  Two 
options  are  available.  When  the  Boolean  TLINP  is  F,  the  ply  thicknesses  are  computed 
internally.  In  this  case^  the  values  entered  do  not  correspond  to  the  actual  ply  thick¬ 
nesses.  When  the  Boolean  TLINP  is  T,  the  ply  thicknesses  are  supplied  through  the  in¬ 
put.  In  this  case  the  values  entered  correspond  to  the  ply  actual  thicknesses.  The  first 
value  entered  is  the  thickness  for  the  bottom  or  inner  ply  of  the  composite.  The  last 
value  entered  is  for  the  last  ply  and  the  number  of  values  entered  equals  the  number  of 
plies  in  the  composite.  (See  tables  I  and  II.) 

(19)  Ply  temperaiure  difference.  -  The  ply  temperature  difference  AT,,  (the  dif- 

- : - - 

ference  between  cure  or  processing  temperature  and  i  ply  temperature)  is  entered  in 
this  group  of  cards.  The  first  entry  is  for  the  first  ply  (which  is  the  bottom  or  inner  ply);, 
and  the  last  entry  is  for  the  last  ply.  The  number  of  entries  equals  the  number  of  plies. 
(See  tables  I  and  11. )  There  are  three  special  cases  associated  with  the  temperature  dif¬ 
ference  in  addition  to  the  general  case  just  described: 

(a)  The  residual  stress  case  at  room  temperature  where  AT^^  equals  the  difference 

between  cure  or  process  temperature  and  room  temperature. 

(b)  The  zero  temperature  effects  case  where  AT,.  =  0. 

th 

(c)  The  no  residual  stress  case  where  AT^^  equals  the  difference  between  i  ply 

temperature  and  room  temperature. 

(20)  Correlation  coefficients  for  strength.  -  The  coefficients  that  correlate  predicted 

amd  measured  values  for  strength  are  entered  in  this  group  of  cards  (see  table  I).  These 
coefficients  are  selected  in  the  same  manner  as  was  described  for  in  card  group  (4). 
The  first  two  entries  are  the  coefficients  and  which  are  for  the  ply 

longitudinal -tensile  strength.  The  third  entry  is  P22T'  transverse - 

tensile  strength.  The  fourth  entry  /^j2s  intralaminar  shear  strength.  The 

fifth  entry  is  for  the  ply  transverse  shear  strength.  The  sixth  entry  is  for 

interply  delamination  limit  strain.  Entries  seven  and  eight  are  the  coefficients  K^i2TT 
and  KLc.rp-,?  which  a.re  for  ply  combined-stress  strength  in  the  tension-tension  and 
tension-compression  quadrants,  respectively.  Entries  9  and  10  ,8^^  and  are  for 

the  ply  longitudinal  compressive  strength.  Entry  11  /^22C  transverse  com_- 

pressive  stress.  Entries  12  and  13  and  ^2  coefficients  for  an  alternate  method 
to  coznpute  the  ply  longitudinal  compressive  strength  (see  section  Subroutine  GLLSC(J)). 
Entry  15  is  entered  as  unity.  This  field  is  allocated  for  possible  future  use.  Entries  15 
and  16  K^2CT  ^n2CC  combined-stress  strength  in  the  compression- 

tension  amd  compression-compression  quadrants,  respectively.  (See  tables  I  and  II.) 

(21)  Constituent  strength  properties.  -  The  constituent  strength  properties  are  enter¬ 

ed.  in  these  two  cards.  The  six  entries  are,  sequentiaily,  in-situ  fiber  bundle  strength 
Srr-,  in-situ  matrix  compressive  strengths  S  p,  in-situ  allowable  matrix  transverse 
tensile  strain  in --situ  allowable  m9.tr3J£  transverse  compressive  strain  e^^^pQ, 
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in-situ  allowable  matrix  shear  strain  e^^pg,  in-situ  allowable  matrix  torsional  strain 
^mpTor'  tables  I  and  n. ) 

(22)  Membrane  loads.  -  The  membrane  (in-plane)  loads  are  entered  in  these  cards. 

The  first  entry  is  the  value  for  for  the  first  load  condition.  The  second  entry  is 

the  value  for  for  the  second  load  condition,  and  so  on  until  N,  values  for  N 

have  been  entered.  Continue  with  N.  values  for  N  and  after  that  with  N,  values 

—  cyy  4C 

for  N^y.  A  total  of  values  are  entered  sequentially.  Note  that  no  empty  fields 

are  allowed  because  they  will  be  interpreted  as  zero  values  for  the  load  conditions  by  the 

code.  Note  also  that  zero  values  for  N_„,  N,„,„,  and  N„^„  have  to  be  entered  even  if 

CXX  cyy  cxy 

the  displacements  are  read  in.  This  is  the  case  when  RINDV  equals  T  (TRUE).  (See 
tables  I  and  H. ) 

(23)  Moments.  -  The  local  bending  moments  are  entered  in  these  cards.  The  des¬ 
cription  is  analogous  to  that  for  the  forces  (card  group  (22)). 

(24)  Displacements.  -  The  local  reference  plane  strains  (e  e  ,  and  e  ) 

-  csxx'  csyy^  cxy' 

and  curvatures  w^^yy,  and  w^^^y)  are  entered  in  these  cards.  The  first  six 

entries  are  for  the  first  load  condition  and  are  entered  sequentially  starting  with  e 

csxx 

and  ending  with  The  next  six  entries  are  for  the  second  load  condition,  and  so  on. 

(See  tables  I  and  n. )  No  blank  fields  are  permitted  in  the  displacement  cards.  Blank 
fields  are  interpreted  to  be  zero  values  by  the  code. 

Note  that  zero  values  are  to  be  entered  for  the  displacements  even  when  the  loads 
are  inputs.  This  is  the  case  when  the  Boolean  (RINDV)  equals  T  (TRUE). 


Input  Ply  Properties 

There  could  be  cases  when  the  user  would  prefer  to  supply  some  of  his  own  ply  prop¬ 
erties  instead  of  using  the  code  to  compute  them.  The  user  has  to  provide  his  own  for¬ 
mats  for  these  cases.  They  are  analogous  to  those  for  reading  in  the  ply  temperature  dif¬ 
ference  AT^j  (card  group  (19)).  The  physical  location  for  these  statements  is  described 
in  the  section  MAIN  PROGRAM  and  by  a  comment  (after  DO  loop  155)  in  the  compiled 
listing  (see  appendix  B). 


Output 

The  program  output  consists  of  printing  out  (1)  the  input  data,  (2)  the  composite 
three -dimension  strain-stress  and  stress-strain  relations  about  the  structural  axes, 

(3)  the  composite  properties  generated  in  array  PC,  (4)  the  composite  constitutive  equa¬ 
tions  about  the  structural  axes,  (5)  the  reduced  bending  and  ajdal  stiffness. 
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(6)  displacement -force  relations,  (7)  the  current  load  or  displacement  condition,  and 
(8)  the  ply  properties  generated  in  array  PL. 

The  printout  of  the  input  data  is  preceded  by  its  code  name.  The  first  and  second 
lines  of  printout  (see  table  XIII  for  corresponding  FORMATS)  are 

THORNEL  -  50/EPOXY 

NL,  NPL,  NPC,  NFPE,  NLC 

8  71  54  1420  1 

The  output  of  the  composite  three-dimensional  strain-stress  temperature  relations 
and  composite  stress -strain  relations  about  the  structural  a,xes  are  printed  under  the 
hes-dings 

3-D  COMPOSITE  STRAIN  STRESS  RELATIONS  ~  STRUCTURAL  AXES 
The  matrices  \EV  and  in  the  equation 

{^c)s  ^t^cls^^'^c)  ■ 

are  printed  out  in  FORMATS  454  and  457  of  subroutine  GACD3. 

3-D  STRESvS  STRAIN  RELATIONS  -  STRUCTURAL  Al^ES 
The  matrix  [E  ]  in 

L*  o 


(^ds  =  t^cL  t^c)s 

is  printed  out  FORMATS  456  and  458  in  GACD3.  The  subscript  s  in  the  preceding 
eqviations  indicates  that  the  relations  are  written  about  the  structural  axes.  It  is  noted 
that  these  properties  are  only  local  to  subroutine  GACD3.  They  can  be  made  global  if 
needed. 

The  output  of  the  composite  properties,  generated  in  array  PC  are  printed  under 
the  heading 

COMPOSITE  PROPERTIES  -  VALID  ONLY  FOR  CONSTANT  TEMPERATURE  THROUGH 
THICKNESS 
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LINES  1  TO  31  3-D  COMPOSITE  PROPERTIES  ABOUT  MATERIAL  AXES 
LINES  33  TO  54  2-D  COMPOSITE  PROPERTIES  ABOUT  STRUCTURAL  AXES 


Fifty -four  entries  are  printed  under  this 

heading  as  follows: 

PC(1) 

weight  density 

PC(2) 

tc 

thickness 

PC(3)  to  PC(ll) 

[^cJ 

three-dimensional  stress -strain  relations 

about  material  axes 

PC(12)  to  PC(14) 

{«c} 

three-dimensional  coefficients  of  expansion 
about  material  axes 

PC(15)  to  PC(18) 

three-dimensional  heat  conductivities  and  heat 
capacity  along  material  axes 

PC(19)  to  PC(30) 

®cll’  ^cl2’  ^cl2 

three-dimensional  constants  about  material 

axes 

PC(31) 

z 

distance  to  reference  plane  from  bottom  of 
composite 

PC(32) 

blank 

PC(33)  to  PC(38) 

two-dimensional  stress -strain  relations  about 

structural  axes 

PC(39)  to  PC(47) 

^cll’  *^cl2’  ^cl2 

two-dimensional  elastic  constants  along  struc¬ 
tural  axes 

PC(48)  to  PC(54) 

(“c)’  Kc  ’  “c 

two-dimensional  coefficients  of  thermal  expan¬ 
sion,  heat  conductivities,  and  heat  capacity 
along  structural  axes 

Array  PC,  its  corresponding  string,  and  headings  are  controlled  by  the  following 
formats  in  subroutine  GOCFD2:  Headings  FORMATS  225,  226,  and  227;  and  string  and 
array  PC  FORMAT  320. 

The  output  for  the  composite  constitutive  equations  are  printed  under  the  heading 
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FORCES  FORCE  DISPLACEMENT  RELATIONS  DISPL  THERMAL  FORCES 


>  = 

[Cox]' 

< 

{^csx} 

>  - 

{^cATx^ 

> 

i“cx}. 

JCckI 

^t^cATx)_^ 

The  elements  of  matrices  A^,  C^,  ^cATx’  ^cATx  printed  out.  The 
FORMATS  are  220  and  330  in  GPCFD2  and  STRING  RESF  in  BLOCK  DATA. 

The  output  for  the  reduced  bending  rigidities  is  printed  under  the  heading 

REDUCED  BENDING  RIGIDITIES 

R. 

The  elements  of  D^  are  printed  out  in  one  line.  The  corresponding  FORMATS  are  355 
and  360  in  GPCFD2. 

The  output  for  the  reduced  axial  stiffness  A^  is  printed  out  under  the  heading 

cx 

REDUCED  STIFFNESS  MATRIX 

The  corresponding  FORMATS  are  364  and  360  in  GPCFD2. 

The  inverse  of  the  constitutive  equations  is  printed  out  under  the  heading 


DISP  DISPLACEMENT  FORCE  RELATIONS  FORCES 


r  ^ 

{^caxl 

’[^cx] 

[Ccxf 

-1 

< 

'{Ncx)^ 

[■’cxJ. 

The  elements  of  this  inverse  are  printed  out.  The  FORMATS  are  682  and  683  in  COMSA 
and  STRING  DISP  in  BLOCK  DATA. 

The  output  for  the  current  load  condition  is  printed  next  to  the  headings 
FOR  THIS  CASE  NBS(X,  Y,  XY-M)  IS 
and 

FOR  THIS  CASE  MBS  (X,  Y,  XY-M)  IS 
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The  current  values  of  N^y,  N^y,  M^,  Mj,y,  and  M^y  are  printed  out  under 

these  headings.  The  FORMATS  are  161  and  162  in  the  main  program. 

The  output  for  the  current  displacement  conditions  is  printed  under  the  heading 

FOR  THIS  CASE  THE  DISPLACEMENTS  DISV(ECSXX,  ECSYY,  ECSXY,  WCBXX, 
WCBYY,  WCBXY)  ARE 

The  FORMAT  is  163  in  MAIN  PROGRAM. 

The  output  of  the  ply  properties  generated  in  array  PL  are  printed  out  under  the 
heading 

LAYER  PROPERTIES,  ROWS -PROPERTY,  COLUMNS -LAYER 
according  to  FORMAT  20  in  MAIN  PROGRAM. 

Seventy -one  entries  are  printed  out  under  this  heading  as  follows: 


PL(1,I) 

^v 

ply  void  content 

PL(2,I) 

kf 

ply  apparent  fiber  content 

PL(3,I) 

kf 

ply  actual  fiber  content 

PL(4,I) 

l^m 

ply  apparent  matrix  content 

PL(5,I) 

ply  actual  matrix  content 

PL(6,I) 

ply  weight  density 

PL(7,I) 

h 

ply  layer  thickness 

PL(8,I) 

h 

ply  and  interply  layer  thickness 

PL(9,I) 

interply  layer  distortion  energy  coefficient 

PL(10,I) 

z 

distance  from  bottom  of  composite  to  ply 

centroid 

PL(11,I) 

^cg 

distance  from  reference  plane  to  ply 

centroid 

PL(12,I) 

^cs 

angle  from  structural  axes  to  composite 
material  axes  (same  for  all  plies),  fig.  2 

PL(13,I) 

angle  from  ply  material  axes  to  composite 
material  axes  (fig.  2) 
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PL(14,I) 

^Is 

PL(15,I)  to  PL(23,I) 

PL(24,I)  to  PL(26,I) 

(“il 

PL(27,I)  to  PL(29,I) 

(K,) 

PL(30,I) 

PL(31,I)  to  PL(32,I) 

''a2’ 

PL(43,I)  to  PL(48,I) 

^IJ.22’  ^/il2'  ^/ii: 

PL(49,I) 

PL(50,I) 

AT 

PL(51,I)  to  PL(60,I) 

^aiT’ 

PL(61,I) 

^a2Q!/3 

PL(62,I) 

PL(63,T) 

- - - 

PL(64,I)  to  PL(69,I) 

PL  (70, 1) 

PL(70,I) 

angle  from  ply  material  axes  to  composite 
structural  axes  (fig.  2) 

ply  stress -strain  relations 

ply  thermal  coefficients  of  expansion 

ply  heat  conductivities 

ply  heat  capacity 

ply  elastic  constants 

ply  strain  magTiifi cation  factors 

interply  delamination  factor 

ply  temperature 

ply  limiting  stresses 

coefficient  in  combined -stress  -  strength 
criterion 

combined -stress  -  strength  criterion 
interply  delarnination  criterion 
ply  applied  strains  and  stresses 
adjacent  ply  relative  rotation 
Hoffmanns  failure  criterion 


The  FORMAT  for  this  output  is  25  and  STRING  PROP  in  MAIN  PROGRAlvI. 


PROGRAM  DESCRIPTION 

The  main  program  (or  control  program.)  and  theoretical  equations  programmed  in 
the  code  are  described  in  this  portion  of  the  report.  The  main  control  program  is  des¬ 
cribed  first,  followed  by  descriptions  of  the  various  subroutines  in  their  physical  se¬ 
quential  order  (fig.  2).  It  is  assumed  that  the  user  of  this  portion  of  the  code  has  a  work¬ 
ing  knowledge  of  computer  programming  and  that  he  is  familiar  with  the  terminology, 
such  as,  micromechanics,  macromechanics,  and  laminate  analysis  of  multilayered  fiber 
composites. 

The  assumptions  and  details  leading  to  the  derivation  of  the  equations  programmed 
in  the  code  are  not  included  here.  However,  they  are  described  in  the  references  cited. 

It  is  suggested  that  the  interested  user  have  these  references  available  to  him. 
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The  information  provided  in  this  portion  of  the  code  together  with  the  compiled  list¬ 
ing  should  be  sufficient  to  enable  the  user  to  modify,  implement,  and  extend  the  code 
according  to  his  needs. 


MAIN  PROGRAM 


The  main  program  contains  the  global  variables,  the  various  subroutines,  the  input 
data  and  format,  the  various  program  control  statements,  and  the  output.  These  are 
discussed  subsequently.  The  flow  chart  of  the  main  program  is  shown  in  figure  6. 

The  global  variables  are  given  in  the  following  list  (for  substitution  and  definition. 


see  appendix  A) : 

Boolean 

TLINP,  CSANB,  BIDE,  RINDV 

Integers 

Np^,  Np^,  Nf,  M,  Q.,  Q^,  Qp,  Q^,  Qf 

Real 

^cs’  ^f’  ^m’  ^f^^’  ‘^)f,m’  ^ 

Real  arrays 
(maximum  dimensions) 

\l’  %  ^Zc’  tr  (1^50),  P^(71K50),  P^(1X54),  E^^,  E^f,  E^^, 

-^CX’  ^cx’^cx’  ^cx’  “^cx  “f’  “m’  “z’  ^cATx’ 

^cATx’  ^csx’  ^cbx(1^3),  ^g(2x8),  ^^(2x10),  ^^(1x4), 
Lg^(lx6), 

String  arrays 

C  (55-  •  spaces  per  field,  composite  system  title)  Read  in. 

P^  (eight  spaces  per  field,  Np^  fields) 

Cgj  (six  spaces  per  field,  six  fields) 

Cg2  (six  spaces  per  field,  six  fields) 

Pep  (six  spaces  per  field,  Np^  fields) 

Current  dimensions 

"r  Npr  Npe,  N, 

Real  arrays 

Sr  Kjr  p/”"  •‘P’ 

(current  dimensions) 


The  subroutines  are  as  follows: 


INVA 

GLLSC 

GACD3 

BLOCK  DATA 


inverse  of  an  array 

generates  ply  stress -limit  conditions 

generates  composite  three-dimensional  elastic  and  thermal  properties 
and  the  two-dimensional  thermal  properties 

DISP  (string)  and  RES  F  (string) 
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GPCFD2 


generates  composite  two-dimensional  elastic  constants  and  constitutive 
equations 

generates  heat  conductivities  of  the  ply 

generates  some  ply  basic  properties  and  the  ply  thermoelastic  con¬ 
stants 

generates  ply  strain  magnification  factors 

generates  the  ply  strain  and  stress  states  due  to  applied  loads  and 
check  for  ply  failure  and  interply  delamination 

These  subroutines  are  described  in  detail  in  the  next  section.  The  strings  of  code  iden¬ 
tifier  DATA  are 

Cg  Read  in  according  to  format  4  in  MAIN  PROGRAM 
P^p  PROP  internally  defined;  PLHD,  PLF,  PLL  output 
Cgj  RESF  internally  defined;  FDRHD,  FDRF,  FDRL  output 
Cg2  DISP  internally  defined;  FDRF,  MDRL  output 
P^,p  PROPC  internally  defined;  PCHD,  PCF,  PCL  output. 

The  strings  and  arrays  P^p  are  printed  out  in  the  main  program,  and  0^2? 

and  P^p  are  printed  out  in  subroutine  GPCFD2.  All  other  input-outputs  are  operated 
by  standard  FORTRAN  formats. 


GPHK 

GECL 

GSMF 
COM  PS  A 


Input 


Control  program 


composite  system  title,  N^,  Np^,  Np^,  Np  (E,  G)^  a^; 

Hkc>  TLINP;  CSANB;  BIDE;  RINDV;  9.  p’  p^,  d.;  k  ,, 

V  '’ic"  'r  D, 

See  portion  of  flow  chart  after  CONTROL  PROGRAM  block  in  fig¬ 
ure  6. 


SUBROUTINE  DESCRIPTION 
Subroutine  INVA(N,  A,  C) 

This  procedure  computes  the  inverse  of  a  square  matrix  A  by  Gauss  elimination 
and  stores  it  in  C.  The  check 

IaM  0 
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is  made  and,  if  satisfied,  the  program  continues;  otherwise,  the  message  "SINGULAR 
MATRIX”  is  displayed.  The  subroutine  inputs  are  N,  A  order  and  array,  respectively. 
The  output  is 


A"^  -  C 


Subroutine  GLISC  (A) 


This  subroutine  generates  the  simple  limit  stress  of  the  single -ply.  The  limit 
stresses  for  the  i^^  ply  are  generated  from  the  following  equations: 


S7I  irr  —  Sfrr  I  ^-frpkf  + 


a  IT  ■  “^fT  '^fT^  ^  '^mT"‘m 


SfllC  =  +  ^fC^ 


^  %1 


(1  -  kj)  +  kj 


V-'H 


1  -  k 


The  second  part  of  the  preceding  equation  was  proposed  in  reference  8 

®aiCD  ""  h®Zl2S  ^  ^2 
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^122C  ~  ^22C 


~nipC 

^v^y.22 


E 


122 


q  _  q  _  fl  J  'mpS 

“^ZISS  ■  '^123 


'112 


v^lil2 


^^233  “  ^233 


mpS 


^,.0 


'm 


v‘^^p23, 


The  transverse  shear  limiting  conditions  for  the  interply  layer  are  not  generated  here. 
However,  provisions  for  them  are  made  in  PL(58,I)  and  PL(59,I)  (where  I  denotes  the 
column  (ply)  index).  The  limiting  stresses  3^^^,^.  -  3^233  <?'pdel  scored  in 
PL(51,I)  to  PL(57,I)  and  in  PL(60,I),  respectively.  The  required  input  to  the  procedure 
is  global  and  is  stored  in  the  following  arrays: 


LSC  ^rnpC’  ^mpS’  ^mptor^ 


BET 


r 

^fT’  ^mT’  ^22T’  ^123’  ^233’  ^del’  ^a2TT^a2TC 
.^fC’^mC’^22C’  ^2’  ^3’  K^12CTKa2CC. 


The  fiber  and  matrix  moduli  are  input  data.  The  ply  moduli  E^22’  ^123 

products  of  are  storedin  PL(32,I),  PL(36,I),  PL(34,I),  and  PL(43,I)  to  PL(48,I), 

respectively.  The  ply  moduli  and  the  strain  magnification  factors  are  generated  in  sub¬ 
routines  GECL  and  GSMF. 


Subroutine  GACD3(C) 

This  subroutine  generates  the  three-dimensional  thermoelastic  properties  of  the 
composite  about  its  structural  (x,y,z)  and  material  (1,2,  3)  a:ces.  The  angle  6  is  meas¬ 
ured  from  X  of  the  structural  axes  system.  (3ee  fig.  5.)  In  figure  5  replace  xx  etc. 
by  11  etc. ,  and  measure  9  from  the  material  axes  for  properties  about  the  material 
axes.  These  composite  properties  are  generated  from  the  following  equations: 
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The  arrays  {a  }  and  {a  }  in  the  preceding  equations  are  given  by 


t 

(a!\=|a!  a  a  a  a  a  I 
I  cJ  L  cxx  cyy  czz  cyz  czx  cxyj 


and 


For  all  practical  purposes  the  two-dimensional  thermal  coefficients  of  expansion  about 

the  composite  structural  axes  are  the  same  as  a  .a  ,  and  a  in  the  array 

cxx’  C5^’  cxy  •' 

for  the  three-dimensional  case. 

The  matrices  [E^],  [E^J,  given  by 
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Note  that  for  the  case  of  an  anisotropic  material,  the  elements  (1,6)  (2,6)  (3,6)  (4,  5),  and 
their  symmetric  parts  will  not  be  zero. 


cos^e  sin^d  0  0  0 

sin^e  cos^0  0  0  0 

[R^j]  =0  0  10  0 

0  0  0  cos  6  sin  9 

0  0  0  -sin  9  cos  6 

-sin  26  sin  20  0  0  0 


where  0  =  0,.  for  properties  about  the  composite  material  and  0  =  0  .  +  0  for  prop- 

61  ^1  CS 

erties  about  the  composite  structural  axes  (see  fig.  5). 
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(sin  20J  -  sin  -(sin  20.  -  sin  20. 


-(sin  20.  -  sin  20. (sin  20.  -  sin  20^_^] 


-[{sin  20.  -  sin  20^  j) 

X  (cos  20.  -  cos  20._^)] 


[(sin  20.  -  sin  20^ 

X  (cos  20.  -  cos  20^  j^)] 


0  0 


0  0 


0  0 


0  0 


0  0 


0  0 


[-(sin  20.  -  sin  20. 

x(cos  20.  -  cos  20.^^)] 
[(sin  20.  -  sin  20^  j) 

x(cos  20.  -  cos  20.  _j^)] 


(cos  20j^  -  cos  20._j)' 


where  i  >  1  and  denotes  the  ply  index.  The  angles  0.  and  0.  ^  (fig.  5)  are  given  by 


^i  =  'a  "  ^cs 


1-1  ^fi-l  ^cs 


The  composite  heat  capacity  is  the  same  for  both  the  three-  and  the  two-dimensional 
cases.  It  is  given  by 


c  i=l 


and  t  is  given  by 


Yj  *Zi 

i=l 


The  composite  three-dimensional  heat  conductivities  along  the  composite  material  axes, 
assuming  an  orthotropic  composite,  are  given  by 


cll  ^  Y  ^122  ®z)i 

c  i=l 
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^022  =  f  I 


]hi(hn 


1 


The  angle  6^  is  measured  from  the  material  axes  (fig.  5). 

The  composite  two-dimensional  heat  conductivities  along  the  composite  structural 
axes  are  given  by  (see  ref.  9  for  the  transformation  equations) 


==i-Vt,.  K,,,  cos^  0  +  Kjoo  sin^  0), 
cxx  I  Zi  Zll  Z22  /I 

C  i:=l 


cyy  t  ^ 
’^c  i=l 


IN? 

-E‘a( 


2  2  ' 
sin  9  +  K^22  ®/i 


Kcyx  =  K^yVL*Zi(,^Z22-%lji--  20. 

c  i=l 


^czz  ^c33 


The  angle  0  in  the  last  set  of  equations  is  measured  from  the  composite  structural  axes 
and  is  equal  to  0gg  +  The  inputs  to  the  subroutine  are  N^,  ^Zi+1’  ^IV  ^cs’  ^Zi’  l-^i^’ 
H.,  h?.,  and  {K?.}  which  are  all  global.  The  variable  N?  is  input  data.  The 

remaining  quantities  are  either  generated  or  are  transferred  from  information  stored  in 
PL(11,I),  PL(13,I),  PL(15-23,I),  PL(8,I),  PL(24,I)  to  PL(26,I)  PL(30,I),  and  PL(27,I) 
PL(29,I).  The  outputs  are  t^  and  the  arrays  are  {a^},  [E^],  h^,  and  (K^). 

The  composite  thickness  t^  is  stored  in  PC(2).  The  arrays  [E^]"^?  and  [E^] 

for  both  composite  material  and  structural  axes  are  printed  out  under  the  headings: 


3-D  COMPOSITE  STPtAIN  STRESS  TEMPERATURE  RELATIONS  -  STRUCTURAL  AXES 
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and 


3-D  COMPOSITE  STRESS  STRAIN  RELATIONS  -  STRUCTURAL  AXES 

The  composite  material  axes  properties  [E  ]  and  {a  }  are  stored  in  PC(3)  to 
PC(14)  as  global  variables.  The  corresponding  moduli  are  stored  in  PC(19)  to  PC(30). 
The  three-dimensional  heat  conductivities  and  heat  capacity  along  the  material  axes  are 
stored  in  PC(15)  to  PC(18).  The  two-dimensional  thermal  coefficients  of  expansion  along 
the  structural  axes  are  stored  in  PC(48)  to  PC(50).  The  two-dimensional  heat  conduc¬ 
tivities  and  heat  capacity  along  the  structural  axes  are  stored  in  PC(51)  to  PC(54).  Note 
that  the  heat  capacity  is  a  scalar  quantity  and  is  independent  of  the  reference  axes. 
Therefore,  PC(54)  equals  PC(18). 


Subroutine  BLOCK  DATA 

In  this  block,  the  strings  C^j  and  C^g  which  are  printed  out  with  the  composite 
constitutive  equations  are  defined.  The  string  C^j  contains  the  resultant  force  notation 

^cx’  ^cy’  ^cxy’  ^cx’  ^cy’  ^cxy'  string  Cg2  contains  the  notation  for  the 
corresponding  displacements. 


Subroutine  GPCFD2  (RES F,  DISP,  PROPC) 

This  subroutine  generates  the  required  section  properties  and  the  force-deformation- 
temperature  relations  for  a  two-dimensional  multilayered  composite.  It  also  generates 
the  plane -stress  elastic  constants  for  the  composite.  The  force -deformation -temperature 
relations  generated  in  this  procedure  are  defined  in  the  following  equation; 


ffNcxl^ 

1 

[Ccx]' 

^csx'' 

rINoATx)' 

/ _ 

1 

[I^cx]_ 

''^cbx 
^  J 

The  generic  equations  for  the  elements  in  the  arrays  [A  ],  [C  ],  [D  ],  (N  .  ™  ),  and 

OX  ox  ox  I  0  ^  JL  X  i 

{^cATx} 
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1^-1 

=  t  HjISj] 


N;  N^l 

tCex)  -  E  ^^a<"a-M  -  "?i)[Ral’'l>=ar'[Rfl] "  L  -rpjHjPjJ 


E''VLi-4>[«al’'['=ar*[Rfll*  E  "rpjH)[Sil 

i=l  j=l 


ly 

cATxl  =  E  ^’'a(“a+i  -  "a)['^al[®ar'  (“a) 

i=l 


(“cATxl  =  E  ^^a<"?i+i  -  =^a>t'*al'^l®ar'  (“al 

i=l 

The  arrays  and  [S-]  are 


l“a)  ' 

1^11  “22 

Oil 

2o 

COS  6 

sin^d 

-  sin  29 

2 

sin^0 

cos^6 

-  -  sin  29 

2 

-sin  26 

sin  29 

cos  20 

26 


1 


0 


[E„]  = 


Sj22  =  S3ii=i(sin  20.  -  sin  20._j)2 


Si21-Sjl2-  -Sjll 


Sj32  =  Sj23  =  ^  (sin  20.  -  sin  20.  _^)(cos  20^  -  cos  20^^^) 


^j31-®jl3-  -Sj23 


S-oQ  =  -  (cos  20.  -  cos  20.  J 

300  A  1  1-i 


Here  0.  equals  the  0  ^  0,  in  figure  5.  The  reduced  bending  rigidities  (ref.  6)  are 

1  CS  6 

generated  in  this  procedure  according  to  the  equation 


>C  1 

X  cxj 


The  reduced  axial  stiffnesses  are  generated  in  the  procedure  according  to  the  equation 


=  ^cx  - 


ic  1 

X  cxj 


The  two-dimensional  composite  elastic  constants  are  generated  from  the  following  equa 
tion  (assuming  AT^.  =  AT  for  i  =  l(l)Np: 
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i=i 


The  inputs  to  this  subroutine  are  t^^  (relative  to  composite  structural  axes), 

H-,  and  the  ply  elastic  constants.  These  quantities  are  global  and  are  located,  respec¬ 
tively,  in  PL(7,I),  (50,1),  (14,1),  (9,1),  and  (31,1)  to  (42,1).  The  arrays  [^^1^ 
and  [Sj]  and  the  dimensions  are  generated  within  subroutine. 

The  outputs  are  the  force-deformation-temperature  relations,  which  are  stored  in 
the  global  arrays  ACX  =  RAC  =  A^^,  CPC  =  C^^,  FLX  =  RDC  =  NSDT  = 

CX  CX  CA  CX  CA. 

^cATx’  MSDT  -  These  a,re  printed  out  under  the  heading 

FORCES  FORCE  DISPLACEMENT  RELATIONS  DISPL  THERMAL  FORCES 
The  reduced  bending  rigidities  are  printed  out  under  the  heading 
REDUCED  BENDING  RIGIDITIES 

The  reduced  axial  stiffnesses  are  printed  out  under  the  heading 

REDUCED  STIFFNESS  MATRIX 

The  inverse  of  the  constitutive  equations 


are  printed  out  under  the  heading 

DISP  DISPLACEMENT  FORCE  RELATIONS  FORCES 
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The  distances  z^,  z^,  z^.  are  stored,  respectively,  in  PC(31),  PL(10,I),  and  PL(11,I). 
The  two-dimensional  composite  stress-strain  relations  is  stored  in  PC(33)  to  PC(38)  and 
the  two-dimensional  composite  moduli  and  Poisson's  ratio  are  stored  in  PC(39)  to 
PC(47).  The  two-dimensional  thermal  properties  are  stored  in  PC(48)  to  PC(54)  as  is 
described  in  the  section  Subroutine  GACD3. 


Subroutine  GPH(CF,  CM,  R,  Q,  CP) 


This  subroutine  is  used  to  calculate  the  ply  conductivities 
specific  equation  programmed  in  this  subroutine  is 


K^22  ^^33' 


The 


K"  =  TC 
laa  mao! 


where  (x  takes  the  values  2  and  3.  The  subroutine  is  called  from  subroutine  GECL. 

The  subroutine  input  variables  CF,  CM,  R,  and  Q  and  the  output  variable  CP  are  defined 
in  the  call  statement  in  GECL.  They  denote,  respectively,  fiber  conductivity 
matrix  conductivity  (modified  for  void  effects),  actual  fiber  volume  ratio  kj,  cor¬ 
relation  factor  computed  conductivity  which  is  the  subroutine  output. 


Subroutine  GECL(KV,  KF) 

The  thermoelastic  properties  of  the  single  ply  are  generated  in  this  procedure.  In 
addition  the  actual  fiber  and  matrix  volume  content,  the  ply  thickness,  density,  and  the 
interfiber  spacing  are  generated.  The  equations  programmed  to  generate  basic  ply  prop¬ 
erties  are 


kf  =  (1.0  -  kY)kj 


Pi  =  Pfh  +  ^m^m 
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The  arrays  in  the  last  two  equations  are  given  by 

1  _  ^'i2i  _ 

^122  ^Z33 

1  _ 

^122  ^Z33 

^in  _  ^122  1 

^Zll  ^Z22  ®Z33 

and 

fa,  j,  )  =  [a, ,  a„,  a,],  „ 

I-  Z,  f,  mJ  I-  V  2’  3JZ,  I,  ni 

The  arrays  [E^],  [Ej],  and  [E^^]  are  generated  locally  in  the  arrays  ECL,  ECF,  and 
ECM,  respectively.  The  arrays  [C^^]  and  [C^^]  and  the  constants  in  them  are  given  by 
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i3j=  1.0 


1/VCF(1, 1) 


VCF(25  1) 


if  VCF(1,1)  0 


if  VCF(1,1)  =  0 


The  variables  VCF(1, 1)  and  VCF(2, 1)  are  empirical  (adjustment)  factors  and  are  read  in. 
Here  and  subsequently,  the  elements  in  the  array  VCF  constitute  experiment-theory  cor¬ 
relation  (semiempirical)  factors  and  are  selected  so  that  the  predicted  and  experimental 
results  for  a  particular  fiber-matrix  system  from  a  particular  fabrication  process  are  in 
good  agreement.  The  variable  /3^  could  be  selected  to  be  different  from  unity  if  addi¬ 
tional  adjustment  is  needed. 

The  elements  in  the  arrays  [C^^]  and  are  generated  by  substituting 


1.0 


l/VCF(l5  4) 


if  VCF(1,4)  0 


VCF(2,4) 


if  VCF(25  4)  =  0 


The  equations  programmed  to  generate  the  shear  moduli  are 


^ml2  T-  ,  ^m 

- -  K  c  +  - 

^t2^  ^r2 

'-'f  ^fl2  ^m 


+1^ 

^t2p,  f~,<2 

"^f  ^fl3  ’^m 
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where 


_^!22L-k,  +  ^2!- 

ptT2p  ptt2 

423  *^111 


c’  = 


|3* 


c;- ^ 


P»»  =1  ^  1  i3” 

Vk^/ 

\  m/ 


The  variables  /3|,  jS;^,  13' and  respectively,  are 

/3|  =  1.0 


'^V/VCF(1.2)  „^cF(l,2),iO 


VCF(2,2) 


if  VCF(2,2)  =  0 


|3’’  =  1.0 
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r 


P''  =  < 

m  ^ 


0\1/VCF(1,3) 


m 


VCF(2,  3) 


if  VCF(1,3)  0 

if  VCF(2,3)  =  0 


The  equations  programmed  for  the  ply  heat  capacity  and  the  ply  heat  conductivities 


are 


H 


cZi  -  —  V“cf^A  +  ^cm^m^m; 


K  =  K 
mao;  m^a 


+  Kv  -  21S,(K„„<,  -  KJ 


and 


^Zll  “  +  *^m^mll 

The  subscript  O'  takes  the  values  (1,  2,  and  3).  The  remaining  variables  are  read  in 
globally  in  the  arrays 


BTA- 


and 


%1 

%2 

^f33 

^cf 

CHK  = 

Kmll 

^m22 

^m33 

Hem 

0 

0 

0 

Hv 

The  small  subroutine  GPHK(CF,  CM,  R,  Q,  CP)  preceding  subroutine  GECL  is  used  for 
programming  convenience  to  compute  the  variables  K^2  ^Z33- 

Inputs  to  subroutine  GECL  are  the  fiber  and  matrix  material  properties  and  the  cor¬ 
relation  factors.  These  properties  are  read  in  globally  and  are  (E,  y,  G,  P,  H  ,  K, 

^^f,  m’  ^f’  *^f’  BTA,  TLINP,  and  (k^,  k^^^  (where  i  =  1(1)N^  and  is  the  number 

of  layers).  For  the  corresponding  code  identifiers,  see  appendix  A. 
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The  outputs  of  subroutine  GECL  are  the  basic  ply  properties  (kp  p^,  t^  (if 
TLINP  =  FALSE),  and  6^  which  are  stored  in  PL(3,I)  and  in  PL(5,I)  to  PL(8,I);  the 
ply  stress-strain  relations,  which  are  stored  in  PL(15,I)  to  PL(23,I);  the  ply  thermal 
coefficients  of  expansion,  heat  conductivities,  and  heat  capacity,  which  are  stored  in 
PL(24,I)  to  PL(30,I),  and  the  ply  moduli  and  Poisson's  ratios,  which  are  stored  in 
PL(31,I)  to  PL(42,I). 

Subroutine  GSMF(SL11,  SL22,  SL12,  SL23,  KV,  KF,  J) 

The  strain  magnification  factors  from  which  the  ply  unidirectional  limiting  stresses 
are  constructed  are  generated  in  this  subroutine.  These  factors  are  <p^22’  *<^^12’ 

^p23  constructing  8^22’  ®U2’  ®Z23’  ^’^spectively. 

Three  methods  are  employed  to  compute  <p^22*  two-dimensional,  Daniel’s 

indirect,  and  Kies’s  one -dimensional  method.  Filament  and  matrix  orthotropicity  and 
the  effects  of  voids  are  included  in  all  of  these  methods  as  is  described  in  reference  1. 
Kies’s  two-dimensional  method  is  selected  to  construct  S^2  “  current  program. 
However,  either  of  the  other  methods  and  even  new  ones  (as  they  become  available)  could 
be  chosen  if,  at  some  future  date,  they  are  found  to  be  more  appropriate.  In  addition, 
optional  degrees  of  freedom  for  adjusting  these  factors  can  be  read  in  globally.  The  op¬ 
tions  are  given  with  the  appropriate  equation.  The  input  and  output  subroutine  informa¬ 
tion  is  discussed  at  the  end  of  this  section. 

The  equations  programmed  in  this  subroutine  are 


k,  =  (1  -  ^>4 


k^  =  (i 


and 
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(1  -  kf) 


strain  magnification  factor  <p^22'  ”  The  three  methods  used  to  compute  the  strain 
magnification  factor  are  given  in  the  following: 

(1)  Kies's  two-dimensional  method: 


if  /3^  0 


i;\i/2 


if  /3^  =  0 


__  (1  -  ^fi2^f2l)^m(g^m22 
■  %12%2p^f(p%22 


1  +  p(A  -  1)  E 


1  -  i^,2iP(^fi2  -  B)  +  J—  p(i..,-  .  B)  -  J  o 


Zl2rZll 


°'Z22  _  ^Z12^Z11 
^Z22  ^Zll 


^IJ.22 


■t  ^m22 


>1.0 


^!Ji22^  -  PL(45,  J) 


(2)  Daniel’s  indirect  method: 
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‘^M22  = 


2  \  ^122 


k  =  0.  83  [/-iV  2  ^  -  1.  35  (l\  ^ 


2  +1.78 


0.  35  <  kj  <  0.  75 


'^M22^v  ^  PL(44,J) 


(3)  Kies's  one -dimensional  method: 


^[i22 


1  -  p/l  - 

\  ‘^f^%22 


•*  PL(43,  J) 


Note  that  PL(46,  J)  is  blank  for  any  other  method  that  might  be  of  interest, 
Strain  magnification  factor  " 


1.0/^’ 


if  /3^  ^  0 


4kAl/2 


y\  if  ^^  =  0 


^IJ.12 


1  -pi 


^m«)^ml2 


^f<p^fl2 


^n2^v  -  PL(47,J) 


Strain  magnification  factor 


^1122- 
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if  jSy  ^  0 


if  -  0 


1 

^li2Z - ^ 

2(1  -  p)  +  (2p  -  1) 

^f(p^f23 


^P23^-PL(48,J) 


Inputs  to  subroutine  GSMF  are  the  ply  applied  stresses  i^-^W  ^122'  ^l\2’  *^f23^’ 

the  void  and  apparent  fiber  content,  the  ply  index,  and  ply,  fiber  and  matrix  elastic  con¬ 
stants.  The  stresses  0-^22’  ^^12  are  transferred  from  PL(67,J)  to  PL(69,J), 

respectively.  (J  denotes  ply  index  in  this  case.)  The  stresses  ^^£3  assigned  the 
value  of  unity.  The  void  and  fiber  contents  are  transferred  from  PL(1,  J)  and  PL(2,  J). 
The  ply  elastic  properties  are  transferred  from  PL(31,J),  PL(32,J),  PL(37,J),  and 
PL(38,  J).  The  fiber  and  matrix  properties  are  read  in  globally.  The  coefficients 
are  in  VCF  as  follows: 


/3 


e 


VCF  = 


L^m 


'S 

m’ 

'm’ 

yft 

^e’ 

y^,  y^*,  0.0, 0.0,  0.0, 


The  outputs  of  subroutine  GSMF  are  the  magnification  factors  stored  in  PL(43-48,  J) 
as  previously  described.  It  is  important  to  note  that  the  magnification  factor  (Pp^22 
pends  on  the  applied  stress  level;  therefore,  GSMF  is  called  from  the  stress  analysis 

subroutine  COMPBA. 


Subroutine  COMPSA(M) 


In  this  subroutine  the  stress  and  strain  state  of  each  ply  are  computed  given  the  edge 
membrane  forces,  the  ply  temperature  and  the  changes  in  curvature.  In  addition,  two- 
ply  combined-stress  strength  criteria  and  the  interply  delamination  criterion  are  gen¬ 
erated.  The  equations  programmed  for  the  i^^  strain  and  stress  states  are 
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+  ^^cATx^  +t^cxJ^^cbx^  >  “^[^a^^^cbx^ 
^''a^  =  t®aJ'^  tR^i][A^^]‘^  ^{N^}  +  {i^cATx^  t^cxJ  {"^cbx^ 

The  reference  plane  strains  and  the  changes  curvatures  are  computed  from 


when  one  or  both  of  the  membrane  force  and  the  moments  are  given. 

The  strains  are  generated  locally  in  EPSL  and  SIGL,  respectively,  and  are  stored  in 
PL(64,I)  to  PL(69,I).  The  matrices  [R^.]  and  [E^.J  are  generated  locally  from  informa¬ 
tion  transferred  from  PL(14, 1)  and  PL(31, 1)  to  PL(42, 1).  The  distance  z^.  and  the  ply 
temperature  AT^j  are  transferred  from  PL(11,I)  and  PL(50,I),  respectively.  The  re¬ 
maining  matrices  are 

-  ACX 
Ccx  -  CPC 

NcATx-NSDT 

Ncx  - 

“cx  - 

^cbx  '^^^m  (local  curvature  from  bending  analysis) 
where  m  denotes  the  load  condition. 
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It  is  important  to  note  that  the  stress  analysis  in  the  coded  form  also  handles  the  case 
where  both  the  reference  plane  membrane  strains  and  the  local  curvatures  are  given.  In 
this  case  the  ply  strains  are  given  by 

^^cxi^  “  '^^csx^  "  ^  ^’''^cbx^ 

where  are  the  i**^  ply  strains  along  the  structural  axis,  are  the  reference 

plane  merSkne  strains,  z  is  the  distance  from  the  reference  plane  to  the  centroid  of 
the  i^^  ply,  and  local  curvatures.  They  are  read  in  the  array  Dyj^ 

where  m  denotes  the  load  condition. 

The  corresponding  i^^  ply  stresses  are  given  by 


where  are  the  i*^^  ply  stresses  along  the  material  axes,  are  the  i*’^  ply  elas¬ 
tic  constants  about  the  ma,terial  axes,  is  the  transformation  matrix  of  the  i  ply, 

{e  .}  are  the  i^^  ply  strains  along  the  structural  axes  as  given  by  a  previous  equation, 

AtT  is  the  temperature  of  the  i^*^  ply,  and  {a,.}  are  the  thermal  coefficients  of  expan- 
IX  fU 

Sion  of  the  i  ply  along  the  material  axes . 

The  displacement  force  relations  are  printed  out  under  the  title 


Displacement  Displacement  force  relations  Forces 

VcxH  t-^cx^l-^cxl 

L  ^ 

Two  similar  sets  are  printed  out.  In  the  first  set  the  displacement  and  force  vectors  are 
in  symbolic  form.  In  the  second  set  the  displacement  and  force  vectors  have  their  nu¬ 
merical  values.  See  outputs  of  trial  cases  (appendix  C). 

The  failure  criterion  may  be  determined  by  either  of  the  following  methods. 

(1)  Modified  distortion  energy 


PL(62,I) 
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The  parameters  ot  and  /3  are  specified  as  follows: 


C  <  0 


T  (T^22  -  ® 


F  ^Z22  ® 


O'  =  T 


min(S^llC.  «  -  C 


®Z22T  ^ 


^Z22C  ^  ^ 


®^Z12Q'^  ■  ^il2Q!/3 


[^Z11®Z22(^  +  '"^12  ^  ^  ^12^. 


lUap 


BET(1,7)  a,l3  =  T 

BET(2,7)  O!  =  C,  iS  =  T 


BET(1,8) 
BET(2,  8) 


O'  =  T,  ^  =  C 

a,  fi  =  C 


The  multiplyer  of  K^j2a/3  generated  in  subroutine  GLLSC  and  is  stored  in  PL(61,I) 
The  constants  '^i2ap  constitute  theory-experiment  correlation  factors. 

(2)  Hoffman’s  criterion  (ref.  9): 


^ZllC  *■  ®Z11CD^ 
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or  as  given  by  the  displacement  force  equation  described  previously. 

The  inputs  to  the  subroutine  are  the  ply  angle  measured  from  the  structural  axes 
from  PL(14,I)),  the  distance  from  the  reference  plane  to  centroid  of  the  ply  (z^j  from 
PL(11,I)),  the  ply  temperature  from  PL(50,I)),  the  interply  delamination  limit 

(A^deij  from  PL(60,I)),  and  the  ply  thermoelastic  properties  stored  in  PL(24  to  26,  -I) 
and  PL(31  to  42,1).  The  ply  extensional  and  coupling  rigidities  =  ACX  and  = 
CPC;  the  local  curvatures  =  WXX;  the  adjustment  constants  K^2TT  “  BET(1,7), 

K’i2ct  =  BET(2,  7),  kJj2tc  "  ^^^(1,  8),  and  K’^2cc  "  BET(2,  8);  and  the  load  condi¬ 
tions  =  NBS(m). 
cx 

The  subroutine  outputs  are  the  modified  distortion  energy  PL(62,I),  Hoffman's  cri- 
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terion  PL(71,I),  the  interply  delamination  criterion  PL(63,I),  and  the  adjacent  ply  rel¬ 
ative  rotation  (A(^j  from  PL(70,I)). 

IMMEDIATE  EXTENSIONS 

The  code  can  be  modified  and  supplemented  to  handle  nonlinear  material  response, 
temperature  dependent  properties,  and  load  envelopes  for  various  angle  ply  composites. 
The  details  of  these  modifications  will  become  apparent  once  the  user  has  some  experi¬ 
ence  in  using  this  code. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  October  7,  1970, 

129-03. 
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APPENDIX  A 


LIST  OF  CODE  IDENTIFIEPvS 


Engineering 

symbol 

FORTRAN 

symbol  code 

Comment 

■^CX 

ACX 

composite  axial  stiffness;  generated  in  sub¬ 
routine  GPCFD2 

aS. 

cx 

RAC 

reduced  axial  stiffness;  computed  in  sub¬ 
routine  GPCFD2 

BIDE 

Boolean 

TRUE  if  interply  effects  are  included;  input 

*"cx 

CPC 

composite  coupling  stiffness;  generated  in 

subroutine  GPCFD2 

Cel 

RESF 

string  with  force  variables  in  BLOCK  DATA 

Ce2 

DISP 

string  with  displacement  variables  in 

BLOCK  DATA 

CSANB 

Boolean 

TRUE  if  membrane  and  bending  s^^mmetry 
exists;  input 

C’s 

S'tring 

composite  title;  MAIN  PROGRAM  formab  4 

FLC 

composite  flexural  rigidities;  generated  in 

subroutine  GPCFD2 

D^ 

cx 

RDC 

reduced  bending  rigidities;  computed  in 
subroutine  GPCFD2 

DISV,DISVI 

displacement  vectors;  DISVI  is  either  read 
in  MAIN  PROGRAM,  or  is  generated  in 
subroutine  COM  PSA 

d, 

DIAF 

filament  equivalent  diameter;  input 

EpEcf 

ECF 

filament  elastic  constants;  generated  in 
subroutine  GECL 

ECL 

ply  elastic  constants;  generated  in  sub¬ 
routine  GECL 

ECM 

mairix  elastic  constants;  generated  in  sub¬ 
routine  GECL 

44 


Engineering 

FORTRAN 

symbol 

symbol  code 

^fll,ai,mll 

EF11,EL11,EM11 

%2,a2,mll 

EF12,EL12,EM12 

PL(9,I) 

“kc 

CHK 

he 

HHC 

i>  j 

I,J 

^cll,c22,c33 

HK11,22,  33 

cxy,  cyy,  cxy 

HK11,22,33 

%1,  Zll,mll 

CHK 

KF,V 

^,m 

KFB,MB 

hfZ,  vl 

KFL,  VL 

Lsc 

LSC 

^cATx 

MSDT 

Mex 

MSB 

m 

M 

Ncx 

NBS 

^cATx 

NSDT 

Comment 

filament,  ply,  and  matrix  normal  moduli; 
filament  and  matrix  moduli  input 

filament,  ply,  and  matrix  shear  moduli; 
filament  and  matrix  shear  moduli  input 

interply  distortion  energy  coefficient;  gen¬ 
erated  in  MAIN  PROGRAM 

array  of  constituents  heat  conductivities; 
input 

composite  heat  capacity  stored  in  PC(18) 
and  PC(54) 

index,  generally  ply  or  interply 

composite  three-dimensional  heat  conduc¬ 
tivities  along  the  material  axes  in  PC(15 
to  17) 

composite  two-dimensional  heat  conduc¬ 
tivities  in  PC(51  to  53) 

see 

apparent  fiber  and  void  volume  ratios;  input 

actual  fiber  and  matrix  volume  ratios 

ply  apparent  fiber  and  void  volume  ratios; 
input 

array  of  limiting  conditions;  input 

thermal  moments;  generated  in  GPCFD2 

applied  moment;  input 

load  condition  index 

applied  membrane  loads;  input 

thermal  force;  generated  in  GPCFD2 
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Engineering 

symbol 

FORTRAN 

symbol  code 

Comment 

NFPE 

number  of  filaments  per  end;  input 

NL 

number  of  plies;  input 

Nzc 

NLC 

number  of  load  conditions;  input 

NPC 

string  PROPC  length;  input 

V 

NPL 

string  PROP  length;  input 

Pc 

PC 

composite  properties  array;  generated  in 
GACD3  and  GPCFD2 

Pz 

PL 

ply  properties  array;  portions  generated  in 
all  parts  of  the  program 

^cp 

PROPC 

string  PROPC;  composite  properties  identi¬ 
fiers  in  GDCFD2 

PZP 

PROP 

string  PROP;  ply  properties  identifiers  in 
MAIN  PROGRAM 

QF,I,P,R,S 

indices  to  print  out  string  PROP 

R 

R 

transformation  matrix;  GACD3,  GPCFD2, 
COMPSA 

RINDV 

Boolean 

T(TRUE)  if  displacements  are  read  in;  in¬ 
put 

Sziixetc. 

PL(51  to  59,1) 

ply  limit  stresses;  generated  in  GLLSC 

*Z 

TL 

ply  thickness;  input  if  TLINP  =  TRUE,  gen¬ 
erated  in  GECL  if  TLINP  =  FALSE 

TLINP 

Boolean 

F(FALSE)  if  ply  thickness  calculated  in¬ 
ternally;  input 

w  u 
cb 

w 

XX 

composite  local  curvatures  relative  to  the 

structural  axes 

0! 

c 

CTE 

composite  coefficient  of  thermal  expansion; 
three-dimensional  in  PC(12  to  14),  two- 
dimensional  in  PC(48  to  50) 
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Engineering 

FORTRAN 

Comment 

symbol 

symbol  code 

Q!x  , 

f,  Z,m 

VAF,AL,AM 

filament,  ply,  and  matrix  thermal  coeffi¬ 
cients  of  expansion;  input  and  VAL  gen¬ 
erated  in  GECL 

VCF 

correlation  factors  for  ply  thermoelastic 
properties  and  strain  magnification  fac¬ 
tors;  input 

BTA 

correlation  factors  for  ply  heat  conduc¬ 
tivities;  input 

BET 

correlation  factors  for  ply  strength;  input 

PL(8,I) 

interply  layer  thickness;  generated  in 

MAIN  PROGRAM 

^csx 

UX 

reference  plane  membrane  strains;  solved 

in  terms  of  N  or  input 

cx 

EPS,PL(64  to  66,1) 

ply  strains;  generated  in  COMPSA 

^cs 

THCS 

angle  between  composite  material  and 
structural  axes;  input 

^li’  ^Ic 

THLC 

angle  between  ply  material  and  composite 
axes;  input 

*'fl2,Zl2,ml2 

NUF12,L12,M12 

filament,  ply,  and  matrix  Poisson’s  ratio; 
input 

71 

PIE 

constant;  input 

^f,m,  Z 

RHOF,M,L 

filament  and  matrix  weight  density;  input 
and  generated  in  GECL 

SIGL,PL(67to  69,1) 

ply  stresses;  generated  in  COMPSA 
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APPENDIX  B 


COMPILED  LISTING 


C  MULTILAYERED  FILAMENTARY  COMPOSITE  ANALYSIS  IS 

C  A  COMPUTER  CODE  FOR  THE  LINEAR  ANALYSIS  OF  MULTILAYERED  FIBER 

C  COMPOSITES.  THE  ANALYSIS  UTILIZES  MICROMECHANICS,  MACROMECHANICS, 

C  AND  LAMINATE  THEORY,  THE  ANALYSIS  IS  RESTRICTED  TO  MEMBRAN E, PL  AT E 

C  AND  THIN  CALLED  SHELL  TYPE  STRUCTURES.  THE  INPUTS  ARE  CONSTITUENT 

C  MATERIAL  PROPERTIES,  CORRELATION  COEFFICIENTS  AND  COMPOSITE  GEOMET 

C  RY.  THE  LOAD  CONDITIONS  ARE  EITHER  FORCES  OR  DISPLACEMENTS  AND 

C  TEMPERATURE  AT  THE  DESIRED  SECTION,  THE  OUTPUTS  ARE  STRESS/STR AIN 

C  /TEMPERATURE  RELATIONS  AND  THEIR  INVERSE,  OTHER  THERMAL 

C  PROPERTIES,  STRENGTH  PROPERTIES,  STRESS  ANALYSIS  RESULTS  AND 

C  THE  MARGIN  OF  SAFETY. 

C  MFC A  -  MAIN  PROGRAM 

LOGICAL  TLINP  ,CSANB,BIDE  ,RINDV 
INTEGER  QI ,QS,OP,aR,QF 

REAL  NLF12,NUF2  3,NUF13,NUF21,NUF32,NUF31  , 

NUL12,NUL23,NUL13 ,NUL21 ,NUL32 ,NUL3I , 

NLM12,NUM23,NUM13 ,NUM21  ,NUM32 ,NUM31  , 

L  SC,  MLR  ,NBS  ,MBS,  l<VL,KF  L,  NSDT,  MSDT  ,Ky  L  J,  KFL  J 
COMMON  /MAGE/ J 
COMMON 

EM22,EM 11,EM23,EM12 ,NUM21 ,NUM12,NUM23 , 

EF2  2  ,EF  lL,EF23,Er-12  ,NUF21  ,NUF12  ,NUF23  , 

EM33,NUM1.3,RH0M,ECM(3  ,3  )  ,eM13  ,  VAM(3  )  ,AXC(3  s  3  )  ,FLC(3,  3) , 
EF33,NUF13,RHOF ,FCF  < 3 ,3 ) ,EF 13 , VAF { 3 5  , 

bET{  2,8  ).,NBS(3, 10)  ,PL(  71  ,50)  ,  WXXOi  ,LSC  (6>  ,  PCiSA  >  ,CPC(3,  3», 
CHK(3,4) ,BTA(4) ,TLINP,OIAF,NFPE,PIE, 

2  MBS (3, 10)  ,RAC  (3,3)  ,DI S VI ( 10  ,6 )  , 

2  CSANB,NPL,NL ,NSDT<3)  ,VCF(2,10)  ,  NUM32 , BI DE , MS DT ( 3 ) , R I  NOV 
DIMENSION  KVL(50)  ,KFL(50)  ,THLC(50}  jTLlSO)  ,MLR(3,10  >,PR0P(71) 

4  FORMAT! 55H  ) 

3  COMTINLE 

READ!  5,4! 

WRITE!  6,15) 

WR  I  TE  !  6,4  ) 

WR  ITE ! 6, 301 

5  FORMAT! 5! 5) 

READ! 5,5)  NL ,NPL , NPC , NFPE , NLC 

10  FORMAT! 51 5) 

11  FORMA T(20H  N L , NP L , NPC , NF PE , NLC ) 

WR  i  TE  !  6,  1 1 ) 

WRITe(6,10)  NL,NPL, NPC, NFPE, NLC 
READ! 5,35)  EFll,EF22,eF33,NUF12  ,NUF23 , NUF 13 , EF 12, EF2  3, EFl 3, 

2  EMI 1,EM22 ,EM33 , NUM12 , NUM2  3 ,NUM13 , EM12 ,  EM23,  EMI  3 

WR ITE ! 6,70) 

WR ITE ! 6,37)  EF 1 1 , EF 22 ,EF33 , NUF 12  ,NUF23  ,NUF13 , EF 12 , EF23 , EF 1 3, 

2  EM11,EM22 ,EM33 ,NUM12 ,NUM23  ,NUM13, EMI 2, EM23, EM 13 

DATA!PROP!I)  ,  I  =  1,71)/6HKV  ,6HKF  ,6HKFB  ,6HKM  , 

2  6HKMB  ,6HRH0L  ,6HTL  , 

2  6HDELTA  ,6HILDC  ,6HZB  ,6HZGC  ,6HTHCS  , 

2  6HTHLC  ,6HTHLS  ,6HSC11  ,6HSC12  ,6HSC13  ,6HSC22  , 

2  6HSC23  ,6HSC33  ,6HSC44  ,6HSC55  ,6HSC66  ,6iHCTEll  , 

2  6!CTE22  ,6HCTE33  ,&HHKl  1  ,6HHK22  ,6HHK33  ,6HHCL  , 

2  6HELI1  ,6HEL22  ,6HEL33  ,6HGL23  ,6HGL13  ,6HGL12 
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2  6HNUL12  ,6HNUL21  ,6HNUL13  f6HNUL31  ,6HNUL23  ,6HNUL32  » 

2  6HSMFK22f  6H  SMFD22  »6H  SMF  S22  1 6H S KFC22  »6HS MFS 1 2  1 6HS MFS  23 f 
2  6HILMFC  f6HTEMPD  ,6HL SCI IT 16HLSCII C t6HLSCl 1 Df6HLSC2 2T, 

2  6H.  SC22C,6HLSC12  ,6HLSC23  ,6HLSCC23 ,6HLSCC13t6HLSC0F  , 

2  6HKL 12ABf 6HMDEIE  ,6HRELR0Tf 6HE PSl 1  f6HEPS22  ,6HEPS12  ,6HSIG11  . 

2  6HSIG22  ,6HSIG12  tGHOELFI  ,6HHFC  / 

15  FORMAT(lHl) 

20  F0RMAT(//47H  LAYER  PROPERTIES,  ROWS-PROPE  RTY ,  COLU  MNS -L  AYE  R ) 

25  F0RMAT(I3,3X,A6,2X, 8614.41 

30  FORMAT!//) 

35  FORMAT! 5E 15. 8) 

37  FORMAT! 10E13. 5) 

41  F0RMAT!4H  VCF ) 

WRITE!  6,41) 

READ!5,35)  !!VCF!I,J),J  =  1,10), I  =  1,2) 

WRITE!  6,37)  !!VCF!I,J),J  =  1,10),  I  =  1,2) 

40  FORMA T!4H  VAF ) 

WRI TE! 6,40) 

READ! 5,35)  ! VAF! I ) ,I  =1,3) 

WRITE!6,37)  {VAF!I),I  =  1,3) 

45  FQRMAT!4H  VAM) 

WR ITE! 6,45) 

READ!*, 35)  !VAM!I),I  =1,3) 

WRITE!6,37)  !VAM!I),1  =1,3) 

50  F0RMAT!59H  THERMAL  CONDUCTIVITIES  AND  HEAT  CAPACITIES  OF  CONST  ITUE 
2NTS) 

55  F0RMATI4H  CHK) 

WRITE! 6,55) 

READ!5,35)  (!CHK!I,J),J  =  1,4), I  =  1,3) 

WRITE!6,37)  !!CHK!I,J),J  =  1,4), I  =  1,3) 

60  FORMAT !4H  BTA ) 

WR ITE ! 6,60) 

READ! 5,35)  IBTA! I ) ,I  =  1 ,4) 

WRITE!6,37)  IBTAII) ,I  =  1,4) 

65  F0RMAT!4H  PIE ) 

WRITE!  6,65) 

READ!5,35)  PIE 
WRITE(6,37)  PIE 

70  FORMA T!/96H  EF 11 ,EF22 ,EF33,NUF12  ,NUF23 ,NUF13 , EF12, EF23, EFl 3, EMI  1, E 
2M22,EM33,NUM12,NUM23  NUM13,EM12 ,EM23 ,EM13) 
eO  F0RMAT1/6H  TLINP) 

WRITE!  6,80 
75  F0RMATIL6) 

READ!*, 75)  TLINP 
WRITE! 6,75)  TLINP 
65  F0RMATI/6H  CSANB) 

WR  ITE! 6,85) 

READ! 5,75)  CSANB 
WRITE!  6,75)  CSANB 

87  F0RMAT!/5H  BIDE) 

WRITE16,87) 

READ!  5,75)  BIDE 
WR1TE!6,75)  BIDE 

88  F0RMAH/6H  RINDV) 

WR  ITE! 6,88) 

READ! 5,75)  RINDV 
WRITE  16,75)  RINDV 
90  F0RMAT!/20H  THC S,RHOF ,RH0M,DI AF ) 

WRITE!  6,90) 
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REAO(5,35)  THC  S.RHOF  ,RHOM,DI  AF 
WRITE(6,37)  THC S.RHOF ,RHOM,DI AF 
S5  FORMAT (AH  KVL ) 

WRITE(fi,  95) 

REAO(5,35)  {KVL(I),I  =i,NL) 

WR  ITE( £,37)  (KVL( I ) ,I  =  1 ,NLI 
100  FQRMAT( AH  KFL ) 

WRITE) £,100) 

READ(  £,35)  (KFLd  )  ,I  =  1  ,NL) 

WRITE(fc,37)  (KFLd  )  ,I  =  1  ,NL  ) 

105  FORMA T(5H  THLC ) 

WR  ITE( £,105) 

R£AO( 5,35)  ( THLC (I ) , 1  =  1 ,NL) 

WRITE(£,37)  (THLC(I),I  =  1,NL) 
lie  FORMA T(3H  TL ) 

.WRITE(£,110) 

READ)  5,35)  (TLd  )  ,I  =  1  ,NL) 

WRITE(£,37)  (TLd), I  =  1,NL) 

111  F0RMAT(6H  PTEMP) 

WRITE!  £,111) 

REAO(5,35)  (PL)50,I),  1=1, NL) 

WRITE(£,37)  (PL(50,I),  1=1, NL) 

115  F0RMAT(/4H  BET) 

WRITE! £,115) 

READ(5,35)  ((BET(I,J),J  =  1,8), I  =  1,2) 
WRITE(£,37)  ((BET(I,J),J  =  1,8), I  =  1,2) 

12C  F0RMAT(/4H  L SC  ) 

WRITE) £,120) 

READ(5,35)  (  L  SC  d  )  ,I  =  1,6) 

WR  ITE  (  £,37  )  (LSCd  )  ,I  =1,6) 

13C  F0RMAT(/4H  NBS) 

WRITE!  £,  130). 

ReAD(5,35>  ((NBS(I,J),J  =  1,NLC),I  =  1,3) 
WRITE(6,37)  ((NBS(I,J),J  =  1,NLC),I  =  1,3) 

131  FDRMAT(/4H  MBS) 

WRITE) £, 131) 

READ(5,35)  ((MBS(I,J),J  =  1,NLC),I  =  1,3) 
WRITE(£,37)  ((MBS(I,J),J  =  1,NLC),I  =  1,3) 

132  F0RMAT)/6H  DISVl) 

WRITE)  £,132) 

READ(5,35)  ) ) 0  I SVl ) I , J) ,  J  =  l,6),  1=  1,NLC) 
WRITe)6,37)  )  (DlSVld  ,J)  ,  J=  1,6),  1=  1,NLC) 
140  CONTINLE 
142  00  145  J  =  1,NL 
PL)  1,  J  )  =  KVL  (  J) 

PL(  2, J  )  =  KFL ( J ) 

PL( 7, J  )  =  TL) J) 

PL) 12, J)  =  THCS 
PL) 13, J)  =  THLC(J) 

PL)  14,  J)  =  THCS+THLC)  J) 

PL)13,J)  =  PL)13,J)*PIE/180.0 
PL(14,J)  =  PL(14,J)<‘PIE/180.0 
KVLJ  =  KVL(J) 

KFLJ  =  KFL(J) 

145  CALL  GECL  (KVL  J, KFLJ) 

DO  155  J  =  2,NL 
PL) 9, J  )  =  0.0 
INE  =  J-1 

IF  (.NOT.  BIDE)  GO  TO  155 
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o  o 


PL(9,J)  =  PL(  8f  J)^PL(8,I  NE) 

PL(  9,  J  )  =  PL(9,J)*PL<9,  J) 

PL(9,J)  =  (2.0^PL(8,J)^PL(8,INE)  )/PL(9,J> 

PL(9,J)  =  0.0186  «(1.0-“PL(9f  Jn  *EM12 
155PL(49,J)  =  0.0093/(PL<8,J>  ^  PL(8,INE)) 

READ  IN  DESIRED  PLY  PROPERTIES  HERE.  SEE  FORMAT  111  AND  THE 
FOLLOWING  THREE  CARDS  FOR  SAMPLE  INPUT. 

CALL  GACD3(3.0) 

CALL  GPCFD2 
DO  1S5  M  =  1,NLC 
CALL  CCMSA(M) 

161  F0RMAT(//33H  FOR  THIS  CASE  N8  S(  X  ,  Y ,  XY-M)  IS  ,  3F10  .0  ) 

162  F0RMAT(//33H  FOR  THIS  CASE  MB  S(  X ,  Y ,  XY-M)  IS  ,  3F10.D) 

163  F0RMAT(//79H  FOR  THIS  CASE  THE  DISPLACEMENTS  DI SV< ECSXX , EC SYY, ECSX 
2Yf WCBXX,WCBYY,WC8XY)  ARE  ,/lH  ,  6E15.5) 

WR  ITE( E, 15) 

IF  (  RINDV)  GO  TO  165 
WRITE(6,161)  (NBS(I,M),  I  ==  1,3) 

WRITE(6,162)  (MBS(I,M),  I  =  1,3) 

GO  TO  166 

165  CONTINLE 

WRITE(^,163)  <DISV1(M,J),  J  =  1,6) 

166  CONTINLE 
WRITE( 6,20) 

WR  ITE( 6,30) 

QF  =  0 

01  =  0 
OR  =  0 
OP  =  0 
OP  =  NL/8 
OR  =  MGD(NL,8) 

IF  (QP  .LT.  1)  GO  TO  185 
DO  175  OS  =  1,QP 
01  =  (QS~1)*8+1 
OF  =  QS*8 
DO  17C  I  =  1,NPL 

170  WRITe(6,25)  I,PROP(I),(  PL(I,J),J  =  QI,OF) 

175  WRITE(6,  15) 

IF(QR  .LE.  0)  GO  TO  185 
01  =  NL--QR+1 
QF  =  NL 

DO  laC  I  =  1,NPL 

180  WRITE(6  ,25)  I  ,PROP  (I  )  ,  <P  L  (I  ,  J)  ,  J  =  01, QF) 

185  IF  ((  QP  .NE.  C)  .OR.  (OR  . LE .  0))  GO  TO  195 

QI  =  1 
QF  =  QR 

00  19C  I  =  1,NPL 

19C  WRITE(6,25)  I  ,  PROP ( I )  ,  ( P L « I , J ) , J  =  QI,QF) 

195  CONTINLE 

WR  ITE( 6,15) 

GO  TO  5 
END 


51 


SUBROUTINE  INVA(N,A,C) 

C  CALCULATE  S  INVA  IN  C 

DIMENSION  A(N,N»  ,C(N,N)  ,6(656)  pD(6,6) 

LOGICAL  TLINP,CSANB,BI0E,R1NDV 

REAL  NIF12,NUF23,NUF13,NUF21 ,NUF32  ,NUF31 , 

2  NLL12,NUL23,NUL13,NUL21 , NUL32 , NUL3 1 , 

2  NLM12,NUM23,NUM13,NUN21 ,NUR32 ,NUM31 , 

2  LSC,  MLR ,NBS,MBS,KVL,KFL,NSDT,MSDT 

COMMON 

2  EM22,EM 11,EM23,EM12 ,NUM21 , NLM12 , NUM23 , 

2  EF22,EF11,EF23,EF12,NUF21 ,NUF12,NUF23 , 

2  EM3  3  ,NUM13,RHOM,ECM (3,3) ,E  M13 , VAM(3 )  ,AXC(3,3),FLC(3,3), 

2  EF33,NUF13,RHOF  ,ECF ( 3,3 ) ,EF13 , VAF (3)  , 

2  BET(  2,  3  )  ,NBS(  3,10)  ,PL  (71  ,50  ,  WXX  (3)  ,LSC  (6)  ,  PC(54  )  ,CPC(3,  3), 
2  CHK( 3,4) ,BTA (4) ,TLINP,DIAF,NFPE,PIE, 

2  MBS (3,10) , RAC (3,3)  ,D I S VI ( 10 ,6 )  , 

2  CSANB,NPL,NL ,NSDT (3 ) ,VCF (2,10) ,NUM32 ,BIDE,MSDT (3 ) ,RINDV 
DO  365  I  =  1,N 
DO  365  J  =  1,N 
B(  I,J)  =  Ad,  J) 

C(  I ,J )  =  0.0 
IF  (  I  .NE  .J)  GO  TO  365 
C(  I, J  )  =  1.0 


365 

CONTINUE 

N  1 

=  N-l 

DO 

39  5  I 

=  1,N1 

DO 

380  K 

=  I  ,N 

IF 

(  B(K, 

I)  .EQa  0. 

SI 

=  B(K, 

I  ) 

DO 

37C  J 

=  I  ,N 

37C 

B( K, J  )  = 

B(K, J)/Si 

DO 

37  5  J 

=  ItN 

275  C(K,J  )  =  C(K, J)/S1 
380  CONTINUE 
IPl  =  I  +  l 
DO  395  K  =  IP1,N 
IF  (B(K,I)  .EQ.  0.0)  GO  TO  395 
00  385  J  =  1,N 
385  8(K,J  )  =  B(K, J)-B(I  ,J) 

DO  39 C  J  =  1,N 
390  C(K,J  )  =  C(K, J)-C (I ,J) 

295  CONTINUE 

C  I  LOOP*** 

SI  =  B(N,N) 

IF  (SI  .EQ.  0.0)  GO  TO  405 
B(N,N  )  =  B(N,N)/S1 
DO  40 C  J  =  1,N 

4CC  C(N,J  )  =  C(N, J)/S1 
405  IF  (SI  .NE.  0.0)  GO  TO  415 
41C  F0RMAT(16H  SINGULAR  MATRIX) 
WRITE(6,410) 

GO  TO  430 
415  DO  420  11 =2,N 
I  =N  +  2-  1 1 
IMl  =  I-l 
DO  420  KK=1, IMl 
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K=I-KK 

DO  421  J  =  1,N 

421  C(K,J)  =  C(K,  J)-C  (I  K,n 

420  8<K,  I  )  =  0.0 

C  END  UPPER  TRIANGLE  REDUCTION 

DO  425  I  =  1,N 

DO  425  J  =  1,N 
D(  I,J  )  =  0.0 

DO  425  K  =  1,N 

425  D(  I,J)  =  D(I,J)+A(I,K)*C(K,J) 
43C  CONTINUE 
RETURN 
END 
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SUBROUTINE  GLLSC(J) 

C  GENERATES  LIMIT  STRESS  CONDITION'S  FOR  SINGLE  LAYER 

LOGICAL  TL  INP  fCSANB  ,B  I  DE  ,  RI  NO  V 
REAL  NLFl 2,NUF23,NUF13,NUF21  ,NUF32 , NUF31 , 

2  NIL12,NUL23,NULL3 ,NUL21 ,NUL32  ,NUL31 » 

2  NLM12,NUM23,NUM13 ,NUM21  ,NUM32 ,NUM31 , 

2  LSC, MLRfNBS.MBS, KVL,KFL,NSDT,MSDT 

COMMON 

2  EM2  2,EM 11,EM23,EM12 ,NUM21 ,NLM12 ,NUM2  3  , 

2  EF22,EF 11,EF23,EF12 ,NUF21 ,NUF12,NUF23 , 

2  EF3  3  ,NUF 13fRHOF ,ECF ( 3 ,3)  fEF13  tVAF (3 )  , 

2  £M3  3,NUM13,RHOM,ECM (3 .3 )  ,E Ml  3 , VAM ( 3 )  ,AXC(3 ,3  )  » FLC(3, ^  ), 

2  BET(2f3),NBS(3fl0),PL(71.50)fWXX(3),LSC(6»,PC(54)fCPC(3f3), 

2  CHK<3,4)  ,BTA (4)  ,TLINP,DIAF,NFPE,PIE, 

2  MB S( 3,10)  ,RAC (3,3)  ,0ISV1(10,6)  , 

2  C  SANB,NPL,NL ,NS0T{3)  ,VCF (2,10)  ,NUM32  ,BI OE,MSDT (3  ) ,RINDV 
PL(51,J)  =  LSC(1)*(BET(1,1)*PL(3,J)+(BET(1,2)*PL(5,J)*EM11/EF11)) 

PL(52,J)  =  LSC(2)*(8ET(2,2)*PL(5,J)+(8ET{2,1)*PL(3,J)*EF11/EM11)) 

PL(54,J)  =  BET(1,3)«(LSC(3)/PL{45,J))*PL(32,J) 

PL(55,J)  =  BET(2,3)*(LSC(4) /PL{45,J)  )<'PL(32,  J) 

PL(56,J)  =  BET(1,4)*(LSC(5)/PL(47, J) ) *PL(36, J) 

PL(53,J)  =  8ET(2,4)*PL(56,J)  +  BET(2,5) 

PL(57,J)  =  BET(1,5)*(LSC(5)/PL(48, J) )*PL(34,J) 

C  FOYE'S  LONGITUDINAL  COMPRESSIVE  STRENGTH  METHOD 


SI 

= 

PL(2,J)*(-1.0  +  EM12 

/EF12)  + 

SI 

EM12/S1 

S3 

= 

PL(1,J)/(1.0  -  PL(2, 

J)  ) 

S2 

= 

1.0  +  S3 

S3 

= 

1.0  -  2.0’S‘S3  +  S3’«'S3 

S4 

= 

S  1+S3/S2/3.0 

END  FOYE'S  METHOD 

I  = 

S2 

PL( 

I 

,  J)  =  AMINKPK  I  ,  J) 

,  S4) 

IF 

(  J 

.LE.  1)  GO  TO  445 

PL{  60 

,J)  =  BET(1,6)4<(LSC(6)/PL(49 

JMl 

J-1 

Si 

= 

PL(lC,JMl)  +  {0.5=f=PL(7 

,  JMl)  ) 

S2  =  C  .25*(PL  (  8,J)-PL  (8,  JMl  )  )<-PC  (31 ) 

ZJ  =  S1+S2 

IF  (ZJ  .GE.  0.0)  GO  TO  435 
S4  =  PC( 31) 

435  IF  (ZJ  .LT.  0.0)  GO  TO  440 
S4  =  PC(2)-PC(31) 

44C  S3  =  ( S4*  S4)-( ZJ*ZJ) 

445  SI  =  (  1.0  +  (  4.  0*PL  (37,  J)  ) -PL  (39,  J)  )  =>=PL  (32  ,  J! 

52  =  ( l.O-PL (41,J) )*PL(31 , J) 

53  =  2.0+PL( 37,J)+PL(37,J) 

S3  =  S2*(  2.0+PL (38,J)  +  PL(41  ,J) ) 

S3  =  S2*PL(31,J)*PL(32,J) 

53  =  SQRT(S3) 

54  =  ( Sl+S2)/S3 
45C  PL( 61, J)  =  S4 

RETURN 

END 
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SieFTC  GACD32  OEbUGtDECK 


SUBROITINE  GAC03(C) 

C  GENERATES  3-D  AXIAL  AND  THERMAL  CONSTANTS 

DIMENSION  EL(6,6»  ,R(6,6)  ,RT(6,6)  ,S  ( 6  ,6  )  ,D1  (6 , 6  )  ,  D2  (  6, 6  ) , 

2  EC(  6,6)  ,EC  I  (6,6>  ,CTL(6)  ,CTC  (6)  ,CTD(6) 

LOGICAL  TLINP,CSANB,BI0E,RINDV,BWECI 
REAL  NLF12,NUF23,NUF13,NUF21 ,NUF32 ,NUF31 , 

2  NIL12,NUL23,NUL13,NUL21 ,NUL32fNUL31  , 

2  NLM12,NUM23,NUM13,NUM21 ,NUN32 ,NUM31 , 

2  L  SC,MLR,NBS,MBS, KVL,KFL,NSDT,MSDT 

COMMON 

2  EM22,EM 11,EM23,EM12 ,NUM21 ,NUM12,NUM23 , 

2  EF22,EF 1 1 ,EF 23 ,EF 12 , NUF2 1 , NUFl 2, NUF23 , 

2  EM33,NUM13,RHOM,ECM(3,3) ,EM13,VAM(3) ,AXC(3,3 ) ,FLC(3,3), 

2  EF33,NUF13,RHOF,ECF(3,3) ,EF13,VAF(3) , 

2  BET<2,8),NBS(3,10),PL(71,50), WXX{3 ) ,LSC (6) , PC(54 ) , CPC( 3,  3), 

2  CHK(3,4)  ,BTA(4)  ,TLINP,OIAF,NFPE,PIE, 

2  MBS(3,10),RAC(3,3) ,DISV1(I0,6) , 

2  CSANB,NPL,NL,NSOT(3)  ,VCF (2,10)  ,NUM32 , BI OE, MSOT ( 3 ) , R I  NOV 
454  F0RMAT(//27X,69H  3-D  COMPOSITE  STRAIN  STRESS  TEMPERATURE  RELATION 
2S  -  STRUCTURAL  AXES//) 

456  FORMAT(//33X, 56H  3-D  COMPOSITE  STRESS  STRAIN  RELATIONS  -  STRUCTURA 
2L  AXES//) 

457  FORMAT (/I  IX, 6E 1 4. 4 , 5X , 1E14. 4 ) 

456  F0RMAT(/21X,6E14.4) 

459  FORMAT(lHl) 

BWECl  =  .TRUE. 

461  DO  435  I  =  1,6 

CTL(  I  )=0.0 
CTC(  I  )=0.0 
CTO(  I  )=0.0 
DO  455  J  =  1,6 
45  5  EC(  I, J  )  =  0.0 
SRC  =0.0 
ST  =  C.O 

DO  462  I  =  15,18 

462  PCI  I )  =0.0 

IF  (.NOT.  SWECl)  GO  TO  464 
DO  463  I  =  51,54 

463  PCI  I)  =0.0 

464  CONTINUE 

DO  5CC  J  =  1 ,NL 
ST  =  ST+PL(7,J) 

SRC  =  SRC+(PL(6, J)*PL(7, J) ) 

ELI  1,  1)  =  PLI15, J) 

ELI  1, 2)  =  PLI16, J) 

ELI  1,  3),  =  PL  (17,  J) 

£LI2,2)  =  PLI18,J) 

ELI2, 3)  =  PL(19,J) 

ELI3,  ')  =  PLI20,J) 

EL  14, 4 )  =  PL (21, J) 

EL(5, 5)  =  PL(22,J) 

ELI6, 6)  =  PL(23,J) 

CTLI  1  )  =  PL I24,J) 

CTLI 2)  =  PL(25,J) 

CTLI  3  )  =  PL (26,  J) 

TH  =  PLI 13, J ) 
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IF  (B\^EC1)  TH  =  PL(14,J) 

R( 2, 2  )  =  COS( TH)«COS( TH) 

R(  1,  1  )  =  R{ 2,2) 

R{  2,  1  )  =  SIN  (  TH)<'S[H(  TH) 

R(  1, 2  )  =  R(2,l) 

R(  3,  2  )  =  1.0 
R(4,4)  =  CUS( TH) 

R  {  b,  4  )  =  -  SIN  (  TH) 

R(  5, 5  )  =  CQ5( TH ) 

R(4,5  )  =  S1N( TH) 

R(l,6)  =  0.5  «SIN(2.0*TH) 

R(  2, 6  )  =  -R ( 1  ,5) 

R  (  6 ,  2  )  =  S I  N  (  2 . 0*  TH  ) 

R(6,  6)  =C0S(2.0*TH) 

R(  6,1)  =  -R ( 6 ,2 ) 

IF  (J  .Lt.  i)  GO  TO  465 
THl  =  PL (  13, J) 

JM 1  =  J-  1 
TH2  =  PLI  13, JMl  ) 

51  =  (  SIN  (  2.  ''*TH1  )-SI  N(2. 0^TH2)  ) 

52  =  (  COSI  2.  THl  )-C0  S(  2. 0«TH2  )  ) 
S(  2,  2  )  =  S1*S] 

S(  1,  I  )  =  S1*S1 
S(  2,  1  )  =  -Sl^Sl 
S(  1,  2  )  =  -S1<=S1 
S(  6,  1  )  =-Sl*S2 
S(  1,  t  )  =-Sl<'  S2 
S(  6,  2  )  =  SI*  S2 
S(  2, 6  )  =  SI*  S2 
S(  o,  6  )  =  S2*  S2 
OD  46C  K  =  1,6 
DO  4  6  C  L  =  1,6 
46C  S(K,L)  =  0.25*S(K,L) 

46  5  on  47 C  K  =  1,6 

DG  47C  L  =  K,6 

47 C  EL(L,K)  =  EL(K,L) 

DO  4  7  5  K  =  1,6 

00  4  7  5  L  =  1,6 

475  kT(K,L)  =  R(L,K) 

00  4  8  C  K  =  1,6 
DO  4  3  C  L  =  1,6 
01(R,L)  =  0.0 


DO  48 C  M  =  1,6 

48Q  01(K,L)  =  01  (  K  ,  L  )  +  (  RT  (  K  ,  M)  *E  L  (  M  ,  1 )  ) 
DO  48  5  K  =  1,6 
00  485  L  =  1,6 
D2( R, L  )  =  0.0 
DO  485  M  =  1,6 

485  D2(K,L)  =  02 ( K , L )  +  ( 01 ( K , H ) *R (  R ,  L  )  ) 


51  =  PL{7,J) 

52  =  PL(9,J) 

DO  490  K  =  1,6 
00  49C  L  =  1,6 

49C  EC(K,L  )  =  EC(K,L) 
PCI  15 )  =  PC ( 15)  + 
PCI  16)  =  PCI  16)  * 
PCI  17 )  =  PCI  17)  + 
PCI  18)  =  PCI18)  + 


+  IS  1*02 (K,L))+IS2*SIK,L)) 
PL (7, J)* I PLI27, J) *Rtl  ,1  )  + 
PL  I  7,  J  )  *  (  PH27,  J  )  *R  12 ,1  )  + 
PL  17, J) /PL{29  ,  J) 
PLI7,J)*PH30,J) 


PLI28,  J)*RI  2,  I )  ) 
PLI28,  J)*RI  1,  1)  > 
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c 


PC(53)  =  PC(53)  PL(7,J)A(PL(28,  J)  ~  PL  (  27  ,  J  )  I  *R  (  6^  2  ) /2  •  C 

DO  5CC  K  =  1, 6 
$3  =  0,0 
DO  49  5  L  =  1,6 

495  S3  =  S3  +  (  S1=^D1  (K,L)«CTL(L)  ) 

5CC  CTD(K  )  =  CTD(KH’S3 
END  J  LOOP 
00  501  K  =  1,6 
DO  501  L  =1,6 

501  EC(K,LI  =  ( l.C/ST)^EC (K,L) 

SRC  =  SRC /ST 
CALL  INVA (6, EC  ,ECI  ) 

DO  510  K  =  1,6 
S3  =  C.O 
DO  505  L  =1,6 

5C5  S3  =  S3+(£CnK,L)^CTD  (D) 

51C  CTC(K)  =  CTC  (K)  +  (  1.0/ST)^S3 
IF  {  .NOT.  BWECi)  GO  TO  506 
PC(  48  )  =  CTC  (  1) 

PC(49)=  CTC{2) 

PC(  50  )  =  CTC(  6) 

PC( 51  )  =  PC( 15)/ST 
PCI  52)  =  PCI 16>/ST 
PCI  53)  =  PCI 53)/ST 
PCI  54)  =  PCI 18)/ST 
5C6  CONTINIE 

IF  I  .NOT.  BWECI)  GO  TO  511 
WRI  TE  I  6,459) 

WR ITEI 6,454) 

WRITe(  6,457)  ((ECItI,J),J  =  1  ,6  )  ,CTC  I  I  )  ,  I  =  1,6) 

WRITE!  6,456) 

WRITE(  6,458)  IIECII,J),J  =  1,6), I  =  1,6) 

WRITE!  6,459) 

BWECi  =  .FALSE. 

GO  TO  461 

TRANSFER  COMPOSITE  PROPERTIES  IN  PC 
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PC(  25) 
PC(  26) 
PC( 27) 
PC(  28  ) 
PC(  29) 
PC(  30  ) 
RETURN 
END 


-Eci  (2a)/ECi  aa) 
-ECI  (ia)/ECI  (2.2) 
-ECI  (3,i)/ECI  (la) 
-ECI  (ia)/EC  I  (3.3) 
-ECI  (3.2)/ECI  (2.2) 
-ECI (2,3)/ECI (3.3) 


$I8FTC  BLOCKl  DECK, LIST 
BLOCK  CATA 

CGPMGN/GPCOM/RESF  (6)  ,DI  SP(6) 


DATA( D  ISP (  I  )  a  = 
6HWYV  ,6HWXY 

1,6) /6HUX 
/ 

,6HVY 

,6HVXPUY  ,6HWXX 

? 

DATA(R£SF(  na  = 

'  6HNY  ,6HMXY 

END 

1  ,6>  /6HNX 

/ 

,6HNY 

,6HNXY  ,6HMX 

9 
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IlM  tM  I'NJ  OU  CM 


SUBRUUIINE  GPCFD2 

C  GENERATES  THE  REQUIRED  SECTICN  PROPERTIES  FOR  LINEAR  BENDING 

C  THEORY  OF  MULTILAYERED  FILAMENTARY  COMPOSITE 

REAL  Ml, NT 

DIMENSION  EL(  3,3)  ,R(3,3)  ,RT(3,3)  tSOfO)  fEC{3,3),CC(3t3)t 
2  FC( 3,3) ,D1(3,3)  ,02(3,3)  ,D3(3,3) ,04(3)  ,MT(3)  ,NT{3), 

2  CTL ( 3),PR0PC (54) ,RDC (3,3) 

LOGICAL  TLINP  fCSANB, BIDE, RIND  V 

REAL  NIFL2,NUF23,NUFI3,NUF21 ,NUF32,NUF31 , 

2  NUL12,NUL2  3,NUL13,NUL2I,NUL32  ,NUL3l  , 

2  NL'M12,NUM23,NUM13,NUM21  ,NUM32,NUM3l  , 

2  L  SC, MLR ,NBS,MBS, KVL,KFL,NSOT,MSOT 

COMMON 

2  EM22,EMli,EM23,EM12 ,NUM2i , NUMi2 , NUM23 , 

2  EF22,EF li,EF23,EF12 ,NUF 2 1 ,NUF 1 2 , NUF23 , 

2  EM3  3,NUM13,RHOM ,ECM(3 ,3)  ,EM13 , VAM(3)  ,AXC ( 3 , 3 ) , FLC ( 3 , 3 ) , 

2  EF33,NUF13,RH0F ,£CF ( 3 , 3 ) , EF 13 , VAF ( 3 ) , 

2  BEK  2,3), NBS(  3,10)  ,  PL  (71, 50)  ,WXX(3)  ,  LSC  (6  )  ,  PC  ( 54  ) ,  CPC  (  3,  3  ) , 
2  CHK(3,4)  ,6TA(4)  , TLI NP, DI AF , NFPE ,PI E, 

2  MBS  (3, 10)  ,RAC  (  3 , 3  )  ,0  I  S  VI  ( 1  0  ,6  )  , 

2  C^AN3,NPL,NL,NSDT(3 )  ,VCF (2  ,10)  , NUM32 , BI DE , MS DT ( 3 ) ,RINDV 
COMMON /GPCOM /RE  SF  ( 6  )  ,  0 1  SP  ( 6  ) 

□ATA{ PROPC ( I ) , I  =  1,54)/6HRH0C  ,6HTC  ,6HCC11  ,6HCC12  , 


2 

6HCC  13 

,6HCC22 

,6HCC23 

,6HCC33 

,6HCC44 

,6HCC55 

,6HCC66 

2 

6HCTE11 

,6HCTe22 

,6HCTE33 

,  6HHK11 

,6HHK2^ 

, 

2 

6HHK33 

, 6 THHC 

,6HEC11 

t6HEC22 

,6HEC33 

, 

2 

6HEC23 

,6HEC31 

,6HEC12 

,6HNLC12 

,6HNUC21 

,6HNUCi3 

f 

2 

6HNUC31 

,6HNUC23 

,6HNUC32 

,6HZCGC 

,6HB2DEC 

,6HCC11 

, 

2 

6HCC  12 

,6HCC13 

,6HCC22 

,6HCC23 

,6HCC33 

,6HEC11 

,6HEC22 

2 

6HEC  12 

,6HNUC12 

,6HNUC21 

f6HCSN13 

,6HCSN31 

,6HCSN23 

, 

2 

6HC  SN32 

,6HCT£il 

,6HCTE22 

,6HCTE12 

,6HHK11 

,6HHK22 

,6HHK12 

2 

6HHHC 

/ 

F0RMAT(A6, 

4X,3E14.4 

,1X,3E14. 

4  ,A6,1E14 

.4) 

1C  FORMAT(/) 

15  FOKMAT(//) 

16  FORMAT(lHl) 

2C  F0RMAT(//7H  FORCE S , 34 X ,2 9H  FORCE  DISPLACEMENT  RELATIONS, 

2  29X,  6H  DISPL,7X,15H  THERMAL  FORCES) 

225  F0RMAT{//77H  COMPOSITE  PROPERTIES  ~  VALID  ONLY  FOR  CONSTANT  TEMPER 
2ATURE  THROUGH  THICKNESS) 

226  F0RMAT{60H  LINES  1  TO  31  3-D  CQ^POSITE  PROPERTIES  ABOUT  MATERIAL 
2AXES) 

22?  F0RMAT(63H  LINES  33  TO  54  2-D  COMPOSITE  PROPERTIES  ABOUT  STRUCTUR 
2AL  AXES) 

CZ  =  0,0 
DO  23C  J  =  1,NL 

C  CG  IS  TAKEN  AT  THE  GEOMETRIC  CENTER 

CZ  =  C2+PL(7,J) 

23CPL(10,J)  =  CZ-(0.5*PL(7,J)} 

PC(  31 )  =  CZ/2.0 

ZBC  =  PC( 31) 

DO  235  I  =  1,NL 

235  PL(  11,  I  )  =  PL ( 10,1 )-ZBC 
C  END  CALCULATIONS  FOR  CG 

DC  24  C  K  =  1,3 
MT(  K  )  =0.0 

240  NT(K)  =  0.0 

00  24  5  K  =  1,3 
00  24  5  L  =  1,3 
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EC( K, L 

)  - 

0*  0 

FCIK.L 

)  - 

0.  0 

CC(K,L 

)  = 

3o  0 

S(K.L  )  = 

L-  •  0 

P J  ^  t  L 

)  = 

a .  J 

F  L  (  K  t  L 

}  = 

0 «  C' 

6  EG 

J  L 

i_)  1 )  p  ^  ^ 

[)u  P9C 

J 

=  1  ,ML 

FL  (1,1 

)  - 

1. 0/PL ( 31  ,  J) 

EL  {  2t  2 

)  - 

i.  3/PL(3.P  .  J) 

EL (  3v  ? 

)  = 

I.'-  /PL  (36,  J) 

EL{  1,  2 

)  ^ 

-PL ( 38, J) /PL (32 

EL (  29  1  )  = 

-PL  (  37, J)  /PL  (31 

CALL  l!\VA  (  3,r:  L  ,eL  ) 

TH  =  P 

L(  i 

4,2  ) 

!<  (  1  ,  1  ) 

= 

cn5(  TH)^Ci)S(  TH) 

E  {  2  9  2  ) 

= 

COS(  )S(  TH ) 

P  (  1 , 2  ) 

= 

SIN(  TH}^S[  N(  TH) 

K  {  2,  1  ) 

= 

SIH  (  TH)^SI  H(  TH) 

E  { i ,  9 ) 

= 

C,3^SIN(2.C4TH) 

W  (  3  9  2  ) 

SIH( 2.0  4TH) 

^(393) 

= 

cns(  2.0  4Th) 

9(^,3) 

-R ( 1 , 3) 

K  (  3  ,  i  ) 

-P (3,2) 

on  23  c 

K 

-  .1,3 

A  p  i5  j; 

L 

1,3 

P  T{ K9  L 

>  = 

K  (  L  ,  K  ) 

I  F  (  J  . 

L  E  0 

I)  GU  TO  233 

Si  =  S 

IN( 

2c-04TH) 

JMl 

J-L 

52  -  S 

INf 

2.0^PL(13, JMl  ) ) 

S3  =  C 

L  S( 

2.0^^  TH  ) 

54  =  C 

CS( 

2.0^:^  PL  (  13,  JMl  )  ) 

5(2,2) 

= 

(  Si-  S2  (  51-52) 

S  (  ]  ,  i  ) 

= 

(  5i-S2)^'(  S1-S2) 

5(2,1) 

= 

-  S  (  1  , 1) 

S  (  1 9  2  ) 

= 

5(2,1) 

S(  3,3) 

= 

(  S3-  54  )  M  S3-  S4) 

S  (  3 ,  i  ) 

- 

-  (  Si-S2  ) {  S3-  54) 

S  (  i  ,  ) 

= 

-  (  SI-  S2  )  4'  (  S3-  54) 

5(3,2) 

= 

{  Si~S2)^  (  S3-S4) 

5(2,  .3  ) 

_ 

5(3,2) 

235 

S4 

= 

C 

.  5 

PL 

(  7,  J) 

51 

= 

P 

L  (  IC, 

J  )-PC  ( 

31)4- 

54 

52 

P 

L  (  i 

0, 

J)-PC  ( 

31  )- 

54 

S3 

= 

c 

«  ^  ' 

(  SI*  Sl-S 

2-  52 

) 

IF 

( J 

.LE 

0 

i  )  GO 

TO  2  63 

DO 

26 

K 

.= 

1,3 

DO 

26 

c 

L 

= 

1  ,  3 

2  6  C 

S(  K 

9  L 

) 

= 

c. 

25>?^PL  ( 

9,  J) 

^  S{K,L) 

56 

= 

p 

•  25 

*  ( 

PL  (  8  ,  J )-PL 

(8, JMl) ) 

56 

5 

6 +  5 

4+PL (10, 

JMl  ' 

-PC  (31  ) 

265 

IF 

(  J 

•  GT 

1)  GO 

TO  2 

TO 

S6 

= 

C  oO 

27C 

QG 

27 

5 

K 

= 

1,3 

DO 

27 

c; 

L 

1,3 

DUK 

y 

L  ) 

= 

0.  0 

on 

2 

7 

f  M 

1  ,3 
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275  OKKtL)  =  OL  (K,L)  4- (  RT(  ’{‘EKM,  L)  ) 

DG  2B0  K  =  1  ,  3 
00  28C  L  =  1,3 
02( K, L  )  =  0.0 
DO  28 C  M  =  1,3 

28C  02(K,L)  =  D2(  K,L)^-(DL  (K,^>1)^R{  M,L)  ) 

00  285  K  =  1,3 
OU  285  L  =  1,3 
S7  =  Q.O 

S7  =  (  51-  S2)^-=02{K,L)^‘S(K,L) 

EC(  K,L  )  =  EC  (K,L)  4-  S7 
S7  =  C.O 

S7  =  S5*D2(K,L)  +  S6=«'S(K,L) 

CC( K, L  )  =  CC (K,L) +S7 
S7  =  C.O 

S7  =  {  1  .  J/3.0  (  (  Si=«^Sl^Sl  )  -  (  S2^'S2«SZ  )  )  *02  (  K,L)  +  (  (56  >5^56  )  ^S(  K,l  )  ) 

285  FC(  K,  L )  =  FC ( K,L)  +  S7 
CTL( 1  )  =  PL(24, J) 

CTL<2)  =  PL(25,J) 

CTL<  3  )  =  0.0 
DO  29C  K  =  1, 3 
D4<K)  =  u.O 

DU  291  L  =1,3 

291  D4{K)  =  D41K}  +  (D1(K,L)*CTL{L)  ) 

NT{K)  =  NT(K)+{PL(50,  J)*D4(K)*PL(7,J)  I 
29C  MT(K)  -  MT(K)+( S5*PL(50, J)*D4(K) ) 


END  J 

LOOP***  * 

DO  29  5 

K 

=  1,3 

DO  29  5 

L 

=  1,3 

295 

OK  K,L  ) 

=  EC(K,L)  /CZ 

CALL  INVA ( 3,Dl ,D2 ) 

PC(  33  ) 

DK  1  ,  1  ) 

PC(  34  ) 

= 

01(1,2) 

PC(  35  ) 

= 

01(1,3) 

PC{  36  ) 

= 

01(2,2) 

PC(  37  ) 

- 

01(2,3) 

PC(  38  ) 

= 

01(3,3) 

PC(  39  ) 

= 

1.0/02(1  ,1) 

PC( 40  ) 

= 

1.0/02(2,2) 

PC(  41  ) 

= 

1.0/02(3,3) 

PC( 42  ) 

-□2(2,I)/D2(1,1) 

PC( 43) 

= 

“D2(  1,2)  /D2(2,2) 

PC( 44  ) 

= 

-02( 3,1) /n2(l ,1) 

PC<  45 ) 

= 

-02(  1 , 3)  /D2(3  ,3) 

PC( 46  ) 

= 

-02(3,2)  /02(2 ,2) 

PC( 47  ) 

= 

-02(2,3)  /D2(3 ,3) 

DO  30  5 

I 

=  1,3 

DO  30  0 

J 

=  1,3 

AXC(  I  , 

J) 

=  EC ( I  ,J) 

CPC(  I, 

J  ) 

=  CC  (  I  ,  J  ) 

BCG 

FLC(  I  , 

J) 

=  FC  (  I  ,  J  ) 

NSDTC  I  ) 

=  NT(  I  ) 

305 

NSOT{  I  ) 

=  NT( I ) 

WR ITE( 6, 

225) 

WRI  TE(  6, 

220) 

WRITE(  e, 

227) 

DO  31 C 

I 

=  1,54 

31C 

WRITE(  6, 

320)  I  ,PROPC (I  )  ,PC(I  ) 
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32C  F0RMAT(I3,3X,A6,E14.4) 

WR  I TF  (  6 , 2  1  6 ) 
wRITE(6,220) 

WR I TE ( 6,210) 

33 C  FQRMAT{2X,A6,4X,3E14.4tlX,3EI4.4,3X,A6 ,4XtE14.4) 

DO  335  I  =  1,3 
WRITEC  6,210) 

33  5  WR  ITE(  6,  3  30  )  RE  SF  ( I  )  ,  (  EC  ( I  ,  J  )  ,  J  =  I  , 3  )  ,  (  CC  (  I  ,  J  )  ,  J  =  1,3), 

2  L)  ISP  {  1)  ,NT(  I  ) 

WR  I  TE( 6,2  10) 

DG  34 C  I  =4,6 
IM3  =  1-3 
WRITER  6,210) 

34C  WRITE(  6,330)  RE  SF  ( I  )  ,  (CC  (  J  ,  I  M3  )  ,  J  =  1  , 3  )  ,  (  FC  (  I  M3  ,  J  )  ,  J  =  1,3), 
2  D ISP ( I  ) ,MT( I  M3) 

CALL  INVA ( 3, EC  ,EC  ) 

DO  345  I  =  1,3 

DO  345  J  =  1,3 

FC(  I,  J  )  =  0.  0 
DO  345  K  =  1,3 

DU  345  L  =  1,3 

345  FC(I,J)  =  FC(  I  ,J)+CC(I  ,K)  >5^EC(K,L)^CC(L,J) 

DO  350  I  =  1,3 

DO  35C  J  =  1 , 3 

35C  RQC(I,J)  =  FLC{I,J)  -  FC{I,J) 

WR  ITE  (  6,2  15) 

355  FGRMAT(27H  REDUCED  BENDING  REGIDITIES) 

WR  ITE  (  6,  35  5) 

WRiTE{  6,2  10 

36C  F0KMAT(9t 1 3. 5 ) 

WRITb{6,360  ((ROC(i,J),  J  =  1,3)  ,  I  =  1,3) 

WR  I  TEC  6,215) 

00  361  I  =  1,3 

00  361  J  =  1,3 

36  1  FC{  I, J  )  =  FLC (  I  ,  J  ) 

CALL  INVA  (3,FC  ,FC  ) 

DO  362  [  =1,3 

DO  362  J  =  1,3 

D3{  I,  J  )  =  0.0 
DO  362  K  =  1,3 

DO  362  L  =  1,3 

362  D3(I,J)  =  D3(I,J)  CC  ( I  ,  K  )  C  (  K  ,  L  )  ^C  C  (  L  ,  J  ) 

DO  363  I  =  1,  3 

DO  36  3  J  =  1,3 

363  RAC(I,J)  =  AXC(I,J)  -  D3(I,J) 

WRITE!  6,215) 

36^  FORMAT(25H  REDUCED  STIFFNESS  MATRIX) 

WR ITE ( 6,364) 

WRITE(  6,360)  ((RAC(I,J),J  =  1,3),!  =  1,3) 

WRITE! 6,215) 

RETURN 

END 
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ro  M  NJ 


SUBROUTINE  GPHK (CF ,CM , R, Q,C P) 

REAL  R 

LOGICAL  TLINP,CSANB,BIDE,RINDV 
REAL  NLFl2tNUF23,NUF13 ,NUF2l,NUF32 ,NUF31 , 
NLLl2fNUL23,NUL13,NUL21  , NUL32  , NUL31 , 

NUM12»NUM23 »NUM13 ,NUM21 ,NUM32 ,NUM31, 

LSC,MLR,NBS,MBS,  KVL,KFL,NSDT,MSOT 
COMMON 

2  EM22 ,EM  U,EM23,EM12 ,NUM2l ,NUM12iNUM23  , 

2  EF22,EF ll,EF23,EF12 , NUF2 1 , NUF 1 2, NUF23 , 

2  £M3  3,NUM13,RHaM,ECM(3,3) ,EM13,VAM(3)  , AXC (3 , 3 ) , FLC( 3, 3  ) , 

2  EF33,NUF13,RHOF,ECF ( 3 ,3 ) ,EF 13 , VAF (3 ) , 

2  BET(  2,8  )  ,N8S( 3, 10) , PL ( 71 ,50 ) , WXX (3 ) , LSC (6 ) , PC ( 54 ) , CPC (  3,  3 )  , 
2  CHK( 3,4)  ,BTA (4)  ,TLINP,DIAF,NFPE,PIE, 

2  M8S( 3, 10) ,RAC (3 , 3 )  ,0  I S VI ( 1 0 ,6 )  , 

2  CSANB,NPL,NL  ,NSDT(3  )  ,VCF(2,10)  ,N1JM32  ,  BI  OE ,  MS  DT  <  3  ) ,  R  INDV 
C  GENERATES  PL  V  HEAT  CONDUCTIVITIES 

51  =  (1.0-CM/CF) 

52  =  0*S0RT(R ) 

SI  =  (  ).0/S2)-Sl 

IF(  SI  .GT.0.0)  GO  TO  525 
52C  FORMAT!  23H  BETA  TOO  LARGE  IN  GPHK) 

WRITE! E, 520) 

525  IF  (SI  .LE.  0.0)  GO  TO  530 
CP  =  1  .0-S2+(  1.0/Sl) 

530  CP  =  CP*CM 
RETURN 
END 
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SUr'KCL  1  INiE  GtLL(KV,Kr^) 

C  GE'lFHAft'S  hCL  FRG'^  Cnf:STITlJt:M  £3><np::-!rifiS 

i  <  E  A  L  A.  V  >  K.  F  t  K  E  lS  t  K  M  b  »  I  -i  j  » I  /  2  »  1  N  V  E  C  L  «  K 

0  IMtfv  S  inN  CFL  (  3 , 3  )  ,  Vm  L  {  A)  ,  C L  ( 3  ,  3  )  ,  C  |V  L  T  (  ?  ,  3  )  ,  I  1  (  3  ,  3  I ,  I  'i  2(  ^ ) . 

2  rCL  (  2,  3  )  ,2FLT  (  3  ,  5)  ,  !  'JVr.CK  i  ,3  ) 

LOGICAL  TL  ,t  SAN3  .3  lOb  ,  K  I  OP  V 

REAL  '■JLFi2,AUr23,'NUFl3,NUF2i  ,  KOF  32  .N.)F3l, 

2  0  LL  1  2  ,N  UL2  3  .NULL  3  ,NUL2  1.  .  t\OL32  ,  N:.JL3  1  . 

2  NL-U2,-MUM23,NU''113,NU.-^2i  ,\0''V;il  , 

2  L  SC  ,  :’^LR  ,\'FS  , -'il  S,  KVL.KI-  L,  KSuT.’-'o  }T 

CliF  ■•lu  N  /MA  Gt  /  J 
Cl.L'Ai'L  \ 

2  r;’''2  2  ,FA  1  1  ,L«2  3,L:  112  ,rvU''121  ,  N C N' 1  2  ,  N  ) 2  , 

2  t:;F22  ,EF  1 1,EF23,F  F  12  ,  'vJF2  1  ,NUF  1  2  ,N'JF2  3  , 

2  tM3  3  .N!)F  l3  ,  RFin  3  ,LC’-i  (  3  ,3  )  , E  ‘-U  3  ,  V A ■'  (  3  !  ,  Ax,C  (  b  ,  :  )  ,  rLC  (  ' ,  . 

2  F  F3  3  ,NUF  li.RHOF ,ECF  (3,3)  ,EF 1 3  ,  VAF  (3 )  , 

2  ;3ET  (  2,  d  )  ,N8  S(  3, 1'))  ,  PL  (  71.  ,3  w  )  ,  .V  XX  (  3  )  ,  LSC  (t.  )  ,  PC  (  V.  )  ,  C  Pt  (3,  -  )  . 

2  C H Ki3,A)  illTAlA)  ,rLiNL'.3IA^  ,)'iFPi:,Pl{;9 

2  M'3S(  3  ,  1; )  ,r<AL  (  3  ,3  )  ,i)  I  S VI  (  1  3  ,6  )  , 

2  C.  SA  j  li ,  \  P  L  r  \  L.  .  N  3(3  1  ( -''  )  ,  Vb  F  (2.1  )  .  0  ij  P  3  2  .  3 1  lI  -  .  R  i .  T  (  ^  )  .  1  i  . 

IF  (FFll  .NF.  EF22)  GO  TD  333 
NUF2i  =  NUF12 

EFi2  ^  Ff  U/(  2.  >*  (1  .  +NUF13  )  ) 

S"’?  IF  (tFll  .EJ.F.F22)  GJ  TO  C'-l' 

N  UF  2  ]  =  \|  LF  1  2<^  (  EF22/i.  F  1  1  ) 

£  A  C  [  F  (  t.  F  1 1  .  NF  .  F  F  3  3  )  G  O  T  0  3  4  G 

NUr31  =  OLF13 

EF13  =  EF  1  1/ (  2  .  >1=  (  i  ,  +NIJF  13)  ) 

34  3  IF  (FFll  .1.;.  t.F33)  1.0  TO  53' 

>]  UF  3  1  =  C'F  1  3*  (  6F  3  3  /E  F  1  1  ) 

EEC  IF  (1.212  .M,  b.Fi3)  GO  TO  ^55 

'■j  U  2  ?.  2  -  0  L  F  2  3 

F  F  2  3  =  tF  2  2  /  (  2  .  (  i  .  +  N  G'3  2  )  ) 

55  5  IF  (LF22  .  t .  tF33)  GO  Tp  bo.) 

M  OP  3  2  -  0  LF  2  3*  (  E F  3  2./-t  2  2  ) 

5tC  IF  (t'^Il  .NE  .  EH22)  ilLJ  TO  Cob 

MU-121  =  \iUVl2 

cv,l2  =  F'^  11/ I  2.  'UM  1  .  L  +  MU^-^iF  )  ) 

56  5  IF  (ton  .22,  2^22)  ^u.)  TO  57./ 

M  UM  21  =  A  CM  1  2  <=  (  F  ,-l  2  2  /r  MU) 

57C  IFIOvil  .ME,  F:K33)  GG  TO  575 
NOM  3i  =  MLM13 

FM13  =  E'-i  11 /(  2.  G-:' (  L.  d-NUMl  b)  ) 

57  5  IF  (  L M  1 1  .  E  J  .  t  P  ? 3  !  G 0  72  5  H'i 

M  J  ■)  3 1  =  1  U  '-1 1 3  *  (  F  M  3  3  /F  '-1 1  I  ) 

5cC  IF  (fcM22  .Mb.  10  13)  :.r  70  5'.c 

MUM  3 2  =  NLM23 

bM  2  3  =  tM  2  2  /  (  2  .  M=  (  !  .  .  +  .'JM2  3  )  ) 

56  5  IF  Or  M2  .EG.  t  o  3  1  G;.  T,)  5 '7 

MUM32  =  M  LM2r<‘  (  EM33 /L ''G'2  ) 

59l  KFG  =  (  1.  G'-KV  )MKF 

KM6  =  (  1. C-KF  )* (  1. 0-KV) 

CF  =  KFB/Kb 
KM  =  l.C-KF 
PL(  3,  J  )  =  KFG 
PL(4,J)  =  1.  .-KF 

PL (  5, J  )  =  KMB 

PL(6,J)  =  RHOM*KTB  +  fHGF>FKFb 

IF  {.MOT.  TLINP  )  PL  (  7  ,  J)  =  SCk  T(  PIE  *(  F  LOAT  ( ..:Fp  i:  ))/(  A  .  -U,  F- )  )  *  J  I  a  F 
IF  (.NOT.  TLINP)  PL(7,J)  =  PL ( o , J ) / S GP T ( V CF I i , A ) ) 

SI  =  P  IE/4.0/KFB 
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CM 

£9  5  IF 
CM 

fcCC  CM 


SI  =  SQR r< SI )-l. j 
PL(  8, J  )  =  S1*UI AF 

IF  (VCKi.L)  .EQ.  0.0)  Gu  TO  595 
CM  =  (  1 .0 /KMB )** ( 1 . 0/VCF ( 1  , 1 )  ) 

IF  (VCF(i,l)  .Nfc.  O.J)  GO  TO  600 
CM  =  VCF(2,I) 

CM  =  (  ]  .  O-KV  )  <‘CM 

A  =  (  (  1.  .'/FFl  l)  +  (  1. 0/tMli  )  *(  KFB/KMB)  ) 
L'.  =  (  (  I.v./lMI  1  )  <■(  1.  C/FFIL  )  *(  KMB/KF8)  ) 


ECF(  1,  ]  1 

L.O/FFU 

ECM (1,1) 

.= 

l.C/k-'ni 

ECF ( 1 »  i  ) 

= 

-M  IfZl/tFZZ 

ECM  (  i  ,  c  ) 

EC  F  (  1 »  3  ) 

= 

-H  LF  31  /FF  3 

fCM(  1,3) 

-MLM3i  /FM33 

ECF(  2, 3  ) 

1.  :)/[:f22 

EC.M  (2,3) 

U’:'/EM22 

ECF (2,1) 

= 

-raFi2/£Fi  i 

ECM (2,  1  ) 

= 

ECF( 2,  3) 

-•N  tF32/£F33 

ECF (2,3) 

= 

-N0M32/hM33 

ECF  (.3,3) 

= 

I.  .)/EF33 

CCM (  ; , 3  ) 

i.  0/£M3'^ 

ECF(  -i,  1  ) 

-  M  EF  L  3  /l  F  i  1 

CCF( 3,  I) 

= 

-NU^^ili/EMI  I 

ECF(  3,  3  ) 

-NLF23/EF22 

EC'-U  3,3) 

-^a^a3/tM22 

DO  61  5  I 

1  ,  3 

DC’  5  L 

= 

I,  3 

CMC (  1  , L ) 

= 

0 .  i  j 

CFL(  1,  L  ) 

= 

*  ' 

CFL(  ! ,  1  ) 

= 

L  .  0  /  (  A  ^  K  MB  ^  £  M 1  1  ) 

CFL (2,3) 

= 

laVCF 

CFL  (  3-^,  3  ) 

r: 

i.O/CF 

CFL  (  1,  3  ) 

— 

(  i,0/A  {  (NUF^i/  (CF^ 

C  FL  (  1 ,  3  ) 

= 

{  i  .  {;  /  A  )^  (  (  NUF  3  i  /  (CF 

CML (1,1) 

= 

1  .  )  /  ( r3’JaFt."M-  F  1  1  ) 

CML (2,3) 

i- G/CM 

CMC  (  .-!,  3  ) 

= 

i. O/CM 

CML  (  1 ,  ^  ) 

= 

(  (  (  NUM2  1  /  (Ci*^ 

CMC  (  1 ,  3  ) 

= 

(  1 .  ')/b  (  N'JM  M  /  (C  M 

DO  6 1 C  I 

= 

1  f  3 

no  6  ^  L  L 

= 

i,  3 

CFL  T(  I  ,L  ) 

=  CFL(L,I  ) 

CML  T(  1  ,L  ) 

=  CML  (L,n 

DO  62 C  I 

- 

1  »  3 

00  62 C  L 

1,  3 

IM2(  1  ,L  ) 

(J  •  0 

IMK  I,L) 

a: 

0.  G 

DO  t  ,1 0  K 

1  t  3 

IMK  (  ,L  ) 

=t 

1  M  1  (  I  ,  L)  +  (  C  F  L  T  (  I  t  K  ) 

IM>(  i,L) 

=i 

IM2(I  a)-^(CMLT(l  ,K) 

DO  6 .5  5  I 

1,3 

DO  635  L 

= 

It  3 

S2  =  C.O 
SI  =  3.0 
DO  63C  K 

1,3 

I  CP^F'''^' 
(  C  3 


[  CF^^EF3 


=  S  l+(  IMI (  I  ,K)*CF  L(K,L)  ) 
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63{  S2  =  S2  +  (  IM2(  I  .Kj'-^^CMKK.L)  ) 

05  ECL(1,L)  =  Sl’i'KFB  +  S.’-K'IB 
ELil  =  l.C'/ECKl.l) 

EL  2 2  =  l.C/ECL(2,2} 
tL'O  =  l.:/cCL(3,3) 

LIJL  12  =  -ECL  (2,1)  /ECL  (1,1) 

NJL21  =  -FCL(  1,2)  /tCL(2,2) 

NUL  13  =  -tCL( 3,1) /ECL (1,1) 

\UL  21  =  -tCL (  1  ,  3) /ECL ( 3, 3) 

!\|JL23  =  -ECL  (  3 . 2) /f-CL  (  2 ,2  ) 

NUL  22  =  -ECL(2,3)/ECL (3,3) 

IE  (VCE(i,2)  .EQ.  ).r,)  gCI  TO  640 

CM  =  {].0/K1  )=J-'4  (1 . 0/vCF  (  1  ,2  )  ) 
t.4C  IF  (VCF(1,2)  .r,F.  ).  0)  GO  TO  645 

CM  =  VCF(2,2) 

64  5  CM  =  (1  .  ,;-KV)  «C  M 

51  =  ( fMl 2/(CF4CF*bFl 2) ) ^KrB 

52  =  KNe/tCM^CM) 

EL  12  -  EM  li;/{  SH-S2) 

51  =  ( fMl 3/(CF-CF*EF13) )*KFy 

52  =  KM:>/(CM-CM) 

EL  1  3  =  E  H  2/(  S1  +  S2) 

IF  (VCK1,3)  .EQ.  O.M)  g  j  to  650 

CM  =  (l.(/K'‘<  )«>!'(  i.  C/VCF  (  1  ,3  )  ) 

c  5  (  I  F  (  VC  F  (  1 ,  2  )  .  N  F  .  0.0)  G  Q  T  C  655 

O’  =  V  C  F  (  2  ,  .3  ) 

655  CM  =  (  1 .0-KV) *CM 

51  =  {  t3M23/(  CF’«’CF*EF23  )  )»KF8 

52  =  K'-i/ (O'K-CM) 
tL23  =  FM  22/(  SI  +  S2) 

IF  (VCK1,4)  .LO.  .'.  ,)  GO  TO  b5o 
CM  =  [J.  ’/KM  )*4 ( 1. C/VCF  (  1  ,  t  )  ) 
c  5  fc  i  F  (  VC  i-  (  1 ,  4  )  .  f:  .  n .  - 1 )  ;;  (_,  j  n  ,,  5  7 
C«  =  VCF( 2,4) 

65?  CM  =  (  i  .  .'-KV  )  4CH 
00  659  I  =  1,3 
00  o  5  F  L  =  1,3 
CFL(  I,L)  =  :  .  n 
6  5  6  C  M  L (  I , L )  =  0 . 0 

CFL(1,])  =  1 . 0/ (  A*KMb-E  Ml  1 ) 

CFl ( 2, 1 )  =  X.  J/CF 

CFL ( 3, 1 )  =  i. 3/CF 


66 


CFL(1,2)  =  ( 1.C/A)*((NUF21/ (CF*EF22) )-(NUM21/ (CM*EM22 n ) 
CFL(1,2)  =  (1.0/A)*((NUF31/((CF*EF33>  )-(NUM31/<CM*EM33»)  ) 
CML(l,l)  =  1.0/(B*KFB*EFli ) 

CML{2, 2)  =  l.O/CM 
CML( 3,  -  )  =  l.O/CM 

CML(1,2)  =  (1.0/B>*((NUM21/(CM*EM22))-(NUF21/<CF*EF22) ) ) 
CML(i,3)  =  (1.0/B>*((NUM31/(CM*E«33))-(NUF31/(CF*Ei-'33)  )) 
DO  660  I  =  lf3 
DO  660  L  =  1,3 
CFLT(  I,LJ  =  CFL(L,I ) 

660  CMLT(  I,L)  =  CMLIL,!) 

DO  661  1=1,3 

VAL{ I  )  =  0.0 

51  =  0.0 

52  =  C.O 

DO  661  L  =  1,3 

51  =  S]+(CFLT( I ,L)*VAF (L) ) 

52  =  S2  +  (CMLT< I  ,L)*VAM(L) ) 

661  VALCn  =  S1*KFB  +  S2*KMB 
CALL  INVA(3,ECL,INVECL) 

PL( 15, J>  =  INVECL(1,U 
PLC 16, J)  =  1NVECL<1,2) 

PL{ 17, J)  =  1NVECL11,3) 

PL( 18, J)  =  INVECL(2,2» 

PL( 19, J)  =  INVECL(2,3) 

PL(20,J)  =  INVECL(3,3) 

PL(21,J)  =  EL23 

PL( 22, J)  =  EL13 
PL{ 22, J)  =  EL23 
PL( 23, J)  =  EL12 
PL<  24, J)  =  VAL<1» 

PLC  25, J)  =  VAL(2) 

PLC  26, J)  =  VALC2) 

PLC  31, J)  =  ELll 
PLC  32,  J )  =  EL22 
PLC  33,  J»  =  EL33 
PLC  34, J)  =  EL23 
PLC  35, J)  =  EL12 
PLC36,J)  =  EL12 
PLC37,J)  =  NUL12 
PLC  39, J)  =  NUL13 
PLC  38, J)  =  MUL21 
PLC  40, J)  =  MUL21 
PL<41,J)  =  NUL23 
PLC42,J)  =  NUL23 

PLC30,J)  =  (CHKCl ,4)*RHQF*KFB+CHKC2 ,4» *RHOM  *KMB > / PL C 6 ,  J ) 
DO  665  L  =  1,3 

51  =  2.0*  BTAC1»*CHKC2,L)*CHKC3,4» 

52  =  KV*C CHKC2,L)-CHKC3,4» ) 

665  CHKC3,L)  =  CHK C 2 , L )* C Sl-2 . 0*S2 ) / C-S2+ SI / BT AC  1 ) ) 

PLC27,J)  =  BTAC2)*KFB*CHKC1 ,1)+C1.0-KF)*CHKC3,1) 

DO  67 C  L  =2,3 
L26  =  26*L 
CHK  IL  =  CHKC 1  ,L) 

CHK3L  =  CHKC3,LI 

LPl  =  L+1 

BTAL  =  BTA(LPl) 

CALL  GPHKCCHKIL ,CHK3L , KFB , B TA L, PLL26 J I 
670  PLC  L26,J>  =  PLL26J 
RETURN 
END 
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SUaRfjLTIME  GSMF(SLli,SL22jSL12,SL?3,KV,KEfJ) 

C  GE'.iEK.AIES  TRANSVEKSt  ANQ  TrtU  SHEAR  M  AGNI  F  It  AT  I  CK  FAC  I'RS 

R  fc  A L  i<  V  *  K  F  ,  K F  a  ,  K HR  ,  MF  2  2  »  HF  1  2  ,  MF  2  3 
DIHFNSIOri  VMFCIOI 
LOGIC  AL  TL  I^Ji^  ,C  SAN3  ,D  I  06  ,R  INJ  V 
REAL  ''iLFi2,NUF23,NUF13  ,IMUF21  ,NiUF32  ,NUF3I  , 

2  ''1  LL  i  2  fNUL23,  NULi3  ,NUL21  ,  N'UL32  ,  NUL31  , 

2  NLMl  2  ,iMUM2  3,iMUf<l  3  ,  NUM2  1  ,  NUf'32  ,NUM31  , 

2  L  SC  ,  MLR  ,NB  S  ,  M!3  S  »  KVL  t  KF  L  ,  KSO  T,  ^'S0T 

COMMON 

2  EM72  ,EM  I  i,EM2  3  ,EM12  ,  NU32I  ,  NUMl  2  ,  i\iUM23  , 

2  F- F2  2  ,  EF  1  1,6  F23  jEF  12  ,  NUF  21  ,  NLF  i  2  ,  NUF2  3  , 

2  EF33,NUF  13,RHnE  ,EtF  (3,3)  ,tF13,VAF(5)  , 

2  EM  3  2  ,N JM  13,RH0H  ,ECM { 3  ,3 )  ,E  Ml  3  ,  VAM ( 3 )  ,  AXC (3 , 3  )  ,  FLC ( 3 , 3  ) , 

2  oET (  2, 8  )  ,NB  S ( 3,  10 ) , PL ( 71  ,50 ) , XX ( 3 )  ,  LSC (F )  , PC ( 54  ) , CPC ( 3,  j ) , 

2  CHK  (  3 , 4  )  ,8  TA  (  4)  ,  TLI  NP  ,  D  I  AF  ,  iNF  Pt ,  P  I  E, 

2  M6S( 3,10)  ,RAC  (3  ,3  )  ,01 SVl ( 10  ,6 )  , 

2  C  SAN  B  ,N  PL  ,NL  ,  NSD  T  ( 3  )  ,  VCF  (  2 , 1  0  )  , NO M32  ,  B I  OE  ,  -iS  DT  (  3  )  ,  K  I N OV 
EL  1  1  =  PL  (  31  ,  J) 
bL22  =  PL (32, J) 

NUL  12  =  PL (37,J) 

NUL  21  =  PL(  38,  J  ) 

VC  =  4.C*KV/(PIE4=  (1.0-KE)  ) 

VC  =  1 .0/ ( 1. 0- SORT (VC ) ) 

KFB  =  (i.O-KV)«KF 

KMB  =  ( 1  .0-KF ) * ( 1 . 0-KV) 

CF  =  KFB/KF 

CM  =  KMB/(  1.0-KF) 

C  KIES  EQUATION 

IE(VCF(1,5)  .EO.  0.0)  GO  Ttj  675 
P  =  SJRT(  4.0/PIE)x<(KFBM<«=(  VCF(1  ,5))  ) 

67  5  IF(VCF(1,5)  .NE.  0.0)  GO  TO  680 

P  =  S JKT(  4.0«KF6/PIE) *VCF  (2  ,5 ) 

E8C  A  =  CM*EM22*(  i.0-iNUFi2>^NUF2i) /(CF*EF22-(1.--  -N1JM12*NUM21  )  ) 

B  =  NUM12--i-'C.M*FM2  2*  (1.  0-NUF  1  2  ^NUF2l  )  /  (  CF*EF22T=  (  1  . f- NUM 1  2  *iNJ  M  2  1 )  ) 

51  =  (  1 .0-(NUL?  14P* (NUF12-B ) )  ) 4 (  1.0/H L22 ) *S L22 

52  =  (  I  .0 /EL  1 1  )  *  (  P*  (  NUFi2-B  ) -NUL12  ) -‘SLII 

53  =  (  1.0/ ( 1 .0  +  P* { A-1 . 0) )  ) * ( Sit S2) 

54  =  (  SL22/EL22)- (NUL12*SL1  1/ELl  I) 

IF  (ABS{S4»  -  C.O)  681  ,  631  ,  632 

68  1  HF22  =1.0 

GO  TO  6  83 
682  MF22  =  S3/S4 

IF  (MF22  .LT.  1.0)  GO  TO  664 
GO  TO  68  3 
664  MF22  =  l.C 
682  VMF (  1  )  =  MF22«VC 
C  DAMN  1  ELS  FUUATIJN 

SI  =  SQRKPIF /KFB  ) 

SI  =  ( SI- 2.0) 

Si  =  0.83^''S1*S1-1. 35*51  +  1.73 

MF22  =  Si*(1.0-NUiM23*NUM32)*(EL22/EM22) 

MF2  2  =  MF22*VC 
VMF( 2  )  =  MF22 

C  KIES  1-D  MAGNIFICATION  FACTOR 

P  =  SQRT(  4.0*KFB/PIE ) * VCF  (2  ,5 ) 

SI  =  CM*FM22/(CF*FF22) 

SI  =  1 .0-P*(  1  .  C-Sl ) 

MF22  =  1.0/Sl 
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FF2  2  =  MF  22*VC 
VMF( 3  )  =  MF22 

IF  (VCF(1,6)  .EU.  0.0)  tiG  TO  695 
p  =  kFB*« ( 1. 0/VCF (1,6)) 

695  IF  (VCF(1,6)  .KF.  0.0)  GO  TO  7'1C 
2  =  SJRT(  4.U*KF3/PIE)*VCF (2  ,6) 

7C-  SI  =  (  CM*tM12/(CF<=EF).2  )  ) 

S3  =  I .0-  (P* ( 1. 0- Si  ) ) 

MF12  =  1.0/ S3 

IF  (VCF(1,7)  .EQ.  0.0)  GO  TO  701 
P  =  KFB*« ( 1 . o/VCF (1,7) ) 

7C1  IF(VCF(1,7)  .we.  0.0)  GU  TO  702 
P  =  SGRT(  A.  j*KFB/PIt )* VCF  (2  ,7 ) 

7C2  SI  =  Cw*tM23/(CF#EF23) 

S3  =  2  .0*  (  1. 0-P  )  +  (  (  2.  0*P-1 . 0) 'i'Sl  ) 
MF23  =  1.0/S? 

VMF(6)  =  MF12*VC 
VMF(7)  =  MF23*VC 
PL(43,J)  =  VMF(5) 

PL( 44, J)  =  VMF (2) 

PL  (  45  ,  J  )  =  V»!F  (  1  ) 

PL  (  46,  J  )  =  ViBF  (  3  ) 

PL(  47,  J  )  =  V'/F  (6) 

PL  (  48,  J  )  =  VMF  (  7) 

f<  t  T  UR  W 

END 
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SUBRCLHNE  COMSACM) 

C  COMPUTES  STRAIN  ANO  STRESSES  IN  THE  LAYERS  GIVEN  THE  PRESCRIBED 

C  EDGE  FORCES, LAYER  TEMP. AND  CURVATURES 

REAL  KL12,NS,LV,MS 

DINENS  ION  AI N(3,3)  ,NS(3 I ,RL (3 ,3)  ,AL(3 ) ,  L( 3, 3  ), 

2CE( 6, 6  ),CE IN ( 6,6)  ,L V(6 )  ,01  SV(6)  ,MS  (  3 )  , 

2  TSK  3),  TS2(3)  ,SIGL(3)  ,EPSL(3) 

LOGICAL  TL  INP  ,C  SANB  ,8  IDE  ,  R  I  NO  V 
67S  FORMAT!  IH  1  ) 

6  8':  FORMAT!/) 

681  FORMAT!//) 

682  F0RMAT!//6H  D I  SP .  ,  34X  ,  2 9H  DISPLACEMENT  FORCE  R EL AT  I CNS , 4 oX , 

2  7H  FORCES) 

68  3  FC!RMATI2X,A6,4X,3E  14.  4 , 1  X  ,  3E  1  4. 4  ,1  6  X  ,  A6  ) 

684  FORMAT!  E 1 2 . 4 , 8X , 3E 1 4. 4 , 1 X , 3E 1 4. 4 , 3 X , Fi 4 . 4 ) 

REAL  NlFl  2,-mUF?  3  ,NUF1  3  ,NUF^1  , NUF3 2  , NUF 3 1  , 

2  NLL12,NUL2  3,NUL13  ,NUL,il  ,  NU  L32  ,  NU  L  3 1  , 

2  NLM12,NUM23,NUM13,NUM21 ,NuM32 ,\UM31  , 

2  L  SC  ,  MLR  ,  NB  S ,  MO  S  ,  K VL  ,  KF  L  ,  NSO  T ,  MSD  T 

(.□'■’HON 

2  FM22  ,EM  1  i,LM2  3,EM12  ,NLM2i  ,  NUM 1 2  ,  NU ''2  3  , 

2  EF2  2,EFll,EF^3,tF12  ,  NUF  2  I  ,  NUF  12,  NUF  23  , 

2  tM32 ,NUM 13, RHUM.ECM  !  3  ,3 )  ,E  M  3 , VAM ! 3 )  ,  AX C ! 3 , 3  )  ,  FLC ! 3 , 3  ) , 

2  cF  3  3  ,NUF  13,i<,Har  ,ECF  (  3 ,3  )  ,EF  13  ,  VAF  |  3  )  , 

2  sET  !  2,  8  )  ,NB  S  I  3  , 1 .')  ,  PL  (  71  ,5  J )  ,  WXX  !3  )  ,  LSC  (6  )  ,  PC  !  54  )  ,cr  C  !  3,  3  )  , 

2  CHK  !  3  ,4)  ,0  TA  (4)  ,  TLI  NP,D!  Ai-  ,  NFPE  ,  P  I  E  , 

2  MBS! 3, 10)  ,RAC  !3 ,3)  ,ni SVl ( 10  )  , 

2  CSANB  ,NPL,NL,NSL)T  !3)  ,  VuF  !2  ,1'' )  ,NCM32  ,0!  Otn  ,  MS  DT  !  3  )  ,K  INOV 
C Of'  ML N  /GP  C  UM  /R E  SF  I  6  )  ,  D  I  SP  !  6  ) 

CALL  INVA  !  3,4.  XC  ,AIN) 

00  685  I  =  1,3 
M  S I  I  )  =  M  B  S !  I  ,  w  ) 

tEE  NS!  1)  =  N3S!  I,-'') 

CC=:, .  : 

DO  673  1=  1,3 

678  CC  =  AMAX 1  ICC ,ABS1MS!  I  )  )  ) 

DC!  68  8  I  =  1,3; 

D  :j  6)6  J  =  1,3 
CE!  I,  J  )  =  AXC  !  I  ,  J) 

JP  3  =  J  +  3 
r  +  3 

Cf-1  I,  JP3)  =  CPC  !  I  ,  J) 

CE(  IP  3,J )  =  CPC ! I , J) 

688  CE!  IP  3 ,JP  2  )  =  FLC  !  I  ,  J) 

CALL  INVA  I6,CE  ,CEIM) 

IF  { .NCT.  RINOV)  GO  TO  677 
DO  675  J  =  1,6 
67  5  OISVI J  )  =  DI SVl !M, J) 

GO  TO  655 
67  7  CUNT  IN  L£ 

DO  689  I  =  1,3 
L  V(  I  )  =  N  SCI  )  +  NS[)TI  I  ) 

IP3  =1+3 

689  LVIIP3)  =  MS!I)  +  MSDTII) 

OD  691  I  =  1,6 

D ISV!  I  I  =  C. ^ 

DO  691  K  =  1,6 

691  DISV(I)  =  DISVII)  +  CEINII  ,K)  x^LVIK) 


70 


69  5  CONTINLE 

DO  692  I  =  1,3 
IP3  =  I  +  3 

652  WXX( I )  =  DISV{ IP3) 

WRITE(  i,679> 

WRITE( 6,682) 

WR  ITE (  6,680) 

DO  686  I  =  1,3 
WRITE( 6,680) 

686  WRITE(  6,683)  D  I  SP  { I )  ,  (CE  I  I  ,  J)  ,  J 
WRITE( 6,680) 

00  687  I  =4,6 
WR ITE (6,680) 

687  WRITE(  6,683)  0  I  SP  ( I  )  ,  ( CE  I  N(  I  ,  J)  ,  J  = 
IF  (R  INOV )  GO  TO  696 

DO  69C  I  =  1,6 
DISV(  I  )  =  0.0 
DO  69C  J  =  1,6 

69C  DISV(I)  =  DISV(I)  CE  IN(  I  ,  J)  ■f'LVIJ) 
GO  TO  693 

696  CONTINLE 

CALL  INVA ( 6,CE IN,CEIN  ) 

00  697  I  =  1,6 
LV(  I  )  =0.0 
DO  697  J  =  1,6 

697  LV(I)  =  LV(I)  +  CEINd  ,  J)  =s=DI  SV(  J) 
CALL  INVA (6,CEIN,CEIN) 

698  CONTINLE 

WR  ITE( 6,681 ) 

WR ITE ( 6,682) 

WRITE! 6,680) 

DO  693  I  =  1,3 
WRITE!  6,680) 

693  WRITE(  6,684)  D  1  SV!  I  )  ,  !CE  I  N  !  1  ,  J)  ,  J  =  1 
WRITE!  6,6  80 

DO  694  I  =  4,6 
WR ITE! 6,680) 

694  WRITE(6,684)  01SV!I),  (CEIN!I,J),J  = 

C  BEGIN  I  -  LOOP 

DO  79C  I  =  1,NL 
TH  =  PL! 14,1 ) 

RL(  2,  2  )  =  CaS!TH)*COS!TH) 

RL(1,1)  =  COS! TH)*COS! TH) 

RL(  1,  2  )  =  SIN!TH)*SIN!TH) 

RL!2,1)  =  SIN!  TH)*SIN!TH) 

SI  =  SIN!  2.0*TH) 

RL  !  2,  3  )  =  -0.  5*Sl 

RL  !  1,  3  )  =  0.  5*51 

RL!  3,  1  )  =  -SI 

RL!  3, 2  )  =  SI 

RL!  3,  3  )  =  CUS!2.0*=TH) 

AL( 1)  =  PL!24, 1 ) 

AL! 2)  =  PL!25,I ) 

AL! 3)  =0.0 
DO  715  J  =  1,3 
00  715  K  =  1,3 
715  EL!J,K)  =  0.0 

EL!  1,  1  )  =  1.0/PL!  31,1) 

EL ! 2, 2  )  =  1.0/PL! 32 , 1 ) 


1 ,6) ,RESF( I ) 


1 ,6)  ,R£SF(  I ) 


,6)  ,LV! I ) 


1 ,6)  ,LV !  I ) 


EL(  3,  3  )  =  1.0/PL{36,I  ) 

FL(1,2)  =  -PL  (38,1  ) /PL(32,I) 

EL(  2,  1  )  =  -PL(37,n  /PL(3i  ,I  ) 

CALL  INVA (3,EL,EL) 

IF  (  R  INDV»  GO  TO  6  99 
IF  (CC  .NE.  C.C)  GO  TO  699 
DO  725  J  =  1,3 
SI  =  G  .0 
DO  720  K  =  1,3 
72G  SI  =  S  1+CPC ( J  ,K 1*WXX( K) 

725  TSKJ  )  =  NS(  J)  +  S1+NSDT(  J) 

00  735  J  =  1,3 
SI  =  C.O 
DO  730  K  =  1,3 
73C  SI  =  SHAIN(J,K)*TS1(K) 

735  TS2(J)  =  S1-PL(11,I)*WXX(J» 

DO  745  J  =  1,3 
SI  =  c.o 
DO  74 C  K  =  1,3 
74C  SI  =  S  H-RL(  J,K)*TS2(K) 

745  EP  SL( J  )  =  SI 

DO  748  J  =  1,3 

748  TSKJ)  =  EPSL  {  J  )- (AL(  J)  >i‘PL(50  ,I  )  ) 

DO  755  J  =  1,3 
SI  =  G.O 
DO  75 C  K  =  1,3 
7  5C  SI  =  S]+FL(J,K)*rSl(K) 

755  SIGL( J  )  =  SI 
GO  TO  700 
699  CONTINLE 

DO  701  J  =  1,3 
TSKJ)  =  C.O 

701  TSKJ)  =  TSKJ)  +  DISV(J)  +  P  L  ( 1 1  ,  I  )  *D  I SV  (  J+3  ) 

DO  703  J  =  1,3 
EPSL(  J  )  =  0.0 
DO  702  K  =  1,3 

7C3  EPSL(J)  =  EPSL(J)  +  R  L  (  J  ,  K)  *  T  SI  (  K) 

DO  705  J  =1,3 
SIGL( J  )  =  0.0 

TS2(J)  =  EPSL(J)  -  PL(50, I ) *AL( J) 

DO  705  K  =  1,3 

705  SIGL(J)  =  SIGL(J)  *  E L ( J , K) ^ T S2 ( K) 

70C  CONTINUE 

CONSTRUCT  LAYER , COMBINED  STRESS  LIMIT  STRENGTH  CRITERION 
SIGLl  =  SIGL( 1) 

SIGL2  =  SIGL(2) 

SIGL3  =  SIGL( 3) 

PL  1 1  =  PL ( 1, I  ) 

PL2I  =  PL ( 2,1 ) 

IF  (A8S( SIGLl)  .EQ.  0.0)  SIGLl  =  .0001 
CALL  GSMF(SIGL1,SIGL2,SIGL3  ,1.0,PL1I  ,PL2I  ,1) 

CALL  GLLSC(I) 

IF  ( { S  IGL ( 1)* SIGL (2  ) )  .LT.  0.0)  GO  TC  765 
IF  (SIGL(l)  .LT.  0.0)  GO  TO  760 
KL12  =  BET(1,7) 

51  =  PL{51,I ) 

52  =  PL(  54,1  ) 

IF(SIGL{2)  .LE.  0.0)  GO  TO  757 
S2  =  PL( 54,1 ) 


72 


151 

76C 


765 


77C 

775 


GO  TO  780 

IF  ISIGLCl)  .GE*  0.0)  GO  TO  765 
KL12  =  BET{2, 8) 

S2  =  PL( 55tl  ) 

SI  =  AMIN  UPL  (52,1  )  ,PL(53  ,I  )  ) 

GO  TO  78C 

IF  (  (  SIGL  (  1)=«‘SIGL(2  ))  .GE.  0.0)  GO  TO  780 
IF  (SIGL(l)  .LT.  0.0)  GO  TO  775 
KL12  =  BET(1,8) 

51  =  PL( 51,1) 

52  =  PL(  55,1  ) 

IF  (SIGL{2)  .GE.  0.0)  GO  TO  770 
S2  =  PL(  55,1  ) 

GO  TO  780 


IF  (SIGL(l)  .GE.  0.0)  GO  TO  780 
KLL2  =  BET(2,7) 

S2  =  PL(  5^,1) 

SI  =  AMIN  KPL  (52,  I  )  ,PL(53,I)) 

GO  TO  780 


78C  KL  12  -  KL  12=5^PL(  61  ,I  ) 

51  =  SIGL( i)/Sl 

52  =  S!GL( 2)/S2 

53  =  SIGL(3)/PL  (56,1) 

54  =  (  5i»Sl)-IKL12=«'Sl^S2)^'(S2*S2)+(S3^S3) 

PL( 62,  i)  =  1.0- S4 

51  =  AMINKPL  (52,1)  ,PL(53,I  )  ) 

52  =  PL(  51,1  )’f'Sl 

53  =  PL(54,I  >’i^PL(55,I  ) 

52  =  (  5IGL(  1  )«SIGL  (l)-SIGL(l  )  ^'SIGKZI^  (SI-PKSl  ,  I)  )*S  IGL(  1  )  )/S2 

53  =  <  SIGL(  2)*SIGL  (2)^*  (PL  (55 , 1) -PL(54 ,1)  )^SIGL(  2)  )/S3 
PL(7l,I)  =  1. 0- (SZ^-SB^-I  SIGL(3)^SIGL  (3  )/PL  (56  ,  I)/PL(  56,  I  )  ) ) 

IF  (  I  .LE.  1)  GO  TO  785 

IMl  ==  I-l 

Si  =  SIN(  2.0«TH)-SIN(2.0^PL(14,IM1)  ) 

S3  =  CGS(  2.0^TH)-CGS(2.0«PL(14,  IMl)  ) 

53  =  TS2(  2)-TS2  (1) 

54  =  C.5*(  S1*S3+S2>«‘TS2(3I  ) 

PLC  7C,  I)  =  S4 

PL(63,I)  =  (PL(60,II-ABS(S4))/Pt(60,I) 

78  5  PL(64,n  =  EPSL(l) 

PL(  65,  1)  =  EPSL(2) 

PL(  66,  I)  =  EPSL(3) 

PL( 67,  n  =  SIGL ( 1 )  • 

PL( 68, I)  =  SIGL(2) 

PL( 69,  1)  =  SIGL(3) 

IF  (  .NOT.  CSANB  )  GO  TO  790 
IF  (I  .NE .  2)  GO  TO  795 
79C  CONTINIE 
END  I-LOQP 
GO  TO  80  5 

795  DO  8CC  I  =  62,NPL 
DO  80C  J  =  3,NL,2 
PL(  I, J  )  =  PL(I  ,1) 

JPl  =  J+1 

IF  (JPl  .GT.  NL)  GO  TO  800 
PL(  I, JPl)  =  PL(I,2) 

80C  CONTINUE 
805  CONTINUE 
RETURN 
END 


APPENDIX  C 


SAMPLE  CASES 

Unidirectional  Composite 


TH0RNEL-5O/EPOXY 


NL  .NPL,NPCtNFPE  ,NLC 

8  71  1-^20  I 

EFil,EF22,EF23,  MUF  12  ,NUF23  ,NUF13  »EF  12  ♦  EF23  , EF 1 3  ,  EMI  1 ,  E M2  2  ,  EM 3 3,  NUH  12 ,  NUH  23  NUM  1  3 ,  EH  1 2 »E M  2 3,  E M 1  3 
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01 
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F 
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F 
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F 
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F 
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0. 
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1  HLC 

00 
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00 
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03 
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01 
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01 
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0. lOOOOE  01 

0-500COE 
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01 
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NBS 

0.50000E 

04 

0. 

0. 

MBS 
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0. 

3-D  COMPOSIT6 

STRAIN  STRESS 

TEMPERATURE 

RELAT  IONS 

-  STRUCTURAL  AXES 
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3-0  COMPOSITE  STRESS  STRAIN  RELATIONS  -  STRUCTURAL  AXES 
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COMPOSITE  PROPERTIES  -  VALID  ONLY  FOR  CONSTANT  TEMPERATURE  THROJ 
LINES  1  TJ  31  3-3  COMPOSITE  PROPERTIES  ABOUT  MATERIAL  AXES 
LINES  33  TO  54  2-0  COMPOSITE  PROPERTIES  ABOUT  STRUCTURAL  AXES 
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GH  THICKNESS 


FORCES 

FURCE  DISPLACEMENT  RELATIONS 

DI  SPL 

THERMAL  FORCES 

NX 
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R  ECUCED  STIFENE  SS  .-lATRI  X 

0.1t221E  07  u.i;>476t  05  -0.  n.l547f>E  'J5  0.61569F  05  0.  -0.  ''i.  0.40568E  ')5 


fHSP.  01  SPLACt'^cNT  FORCE  PFLATIOrjS  FORCES 
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STRUCTURAL  AXES 


0 . 1430E-06 

-0.2449E-06 

0.4471 E-37 

0  . 

0. 

0. 

-0 .2449E-06 

0.7909E-06 

-0.2246E-06 

-0  . 

-0. 

-0. 

0.4471E-07 

-0.2246E-06 

O.9353E-0fa 

0  . 

0, 

0. 

0. 

0, 

0. 

0.29  37E-05 

0. 

C. 

0. 

0. 

0. 

0  , 

0. 2937E-05 

0, 

U  . 

0. 

0. 

0  . 

0. 

0.2000E-06 

3-0  COMPOSITE  STRESS  STRAIN  RELATIONS  -  STRUCTURAL  AXES 


-0.3131E-05 
0.  1199E-D4 
0,2g33E-04 
-0. 

-0. 

-0. 


0.15  04E 

Od 

0.4779E 

07 

0.4285E 

06 

-0. 

-0, 

0. 

0.4779E 

07 

0.2875E 

07 

0,4620E 

06 

-0. 

-0. 

0. 

0.4285E 

06 

0.  46  20E 

06 

0.1160E 

0  7 

-0. 

-0. 

-0. 

-0. 

-0. 

-0. 

0.3405E  06 

-0, 

-0. 

-0. 

-0. 

-0. 

-0. 

0.34O5E  06 

-0. 

0. 

0. 

-0. 

-0. 

-0  > 

0.5001E  07 

COMPOSITE  PRUPi-KTIbS  -  VALID  ONLY  FOR  CONSTANT  TEMPERATURE  THROUGH  THICKNESS 
LINES  1  TO  ii  3-0  COMPOSITE  PROPERTIES  AdOUT  MATERIAL  AXES 
LINES  33  TO  54  ?~‘0  COMPOSITE  PROPERTIES  AdOUT  STRUCTURAL  AXES 


1 

RHUL 

0.5i65E- 

-01 

2 

TC 

0.640CE- 

-01 

3 

CC  11 

C. 1504E 

08 

4 

CC12 

C.4779F 

07 

5 

CC13 

C.42  85E 

C6 

6 

CC22 

0.2875t 

07 

7 

CC23 

C.4620E 

06 

8 

CC3  3 

C.  1160E 

07 

9 

CC44 

C.  3405E 

06 

10 

CC55 

a.3  405E 

06 

1  1 

CC66 

C.5001E 

07 

12 

CTEl  1 

-0.3L31E- 

-05 

13 

CTE22 

0. 11 99F- 

-04 

14 

CTE33 

0.29  33E- 

-04 

15 

Kll 

0.21H9E 

03 

16 

HK22 

C.  7  544E 

02 

17 

HK33 

0. 3715E 

01 

18 

hhC 

0.2043E 

00 

19 

EC  11 

0.6992E 

07 

20 

EC  22 

C.  1264E 

07 

21 

EC  3  3 

0. 1069E 

07 

22 

EC  23 

C,  3405E 

06 

23 

EC31 

0. 3405F 

06 

24 

ECU 

C.  50  016 

07 

25 

NUC12 

C. 17126 

01 

26 

NUC21 

0.30  966 

00 

27 

NUC  13 

-0.31266 

00 

28 

NUC31 

-0.4780E- 

■Cl 

29 

NUC23 

0.28396 

CO 

30 

NUC  3  2 

C. 24016 

CO 

31 

ZCGC 

0.32006- 

•01 

32 

a2DEC 

0. 

33 

CC  11 

U. 14886 

08 

34 

CC  12 

0.4600E 

07 

35 

CC13 

0. 

36 

CC22 

0.2691E 

07 

37 

CC23 

C. 

38 

CC  3  3 

C. 500  16 

07 

39 

EC  11 

G.6992E 

07 

40 

EC  2  2 

0. 1264E 

07 

41 

EC  12 

0.5001E 

07 

42 

NUC12 

0.i712E 

01 

43 

NUC2i 

C.3096E 

00 

44 

CSN13 

-  G. 

45 

C5N31 

-  0. 

46 

CSN23 

-  0. 

47 

CSN32 

0, 

48 

CTEl  1 

-0.31  31E- 

05 

49 

CTE22 

0. 11 99E- 

04 

50 

CTE12 

-  0. 

51 

HKU 

C.21 89E 

03 

52 

HK22 

0. 7544E 

02 

53 

HK  12 

-0. 

54 

HHC 

0,2043E 

00 

78 


FORCES 

FORCE  DISPLACEMENT  RELATIONS 

01  SPL 

THERMAL  FORCES 

NX 

0.9525E  06 

0.2949E  06 

0. 

0.1221E-03 

0.4578E-04 

0. 1068E-03 

UX 

-0.1656E  03 

NY 

0.2949E  06 

0.1722E  06 

0. 

0.4578E-04 

0.2289E-04 

O.3015E-O4 

VX 

-0.3427E  03 

NXY 

0. 

0. 

0.3200E  06 

0. 1068E-03 

0.3815E-04 

0. 3052E-04 

VXPUY 

0, 

MX 

0.12216-03 

0.4578E-04 

0.1068E-03 

0.3251E  03 

0.  1007E  03 

0.6390E  02 

WXX 

-0.2235E-07 

MY 

0.4578E-04 

0,  22  89E-04 

0.3815E-04 

0.1007E  03 

0.5879E  02 

0.2260E  02 

WYY 

-0.2980E-07 

MXY 

O.IO60E-O3 

0.3815E-04 

0.3052E-04 

0  .6390E  02 

0.2260E  02 

0.1092E  03 

WXY 

0.2980E-07 

REDUCED  BENDING  REGIDITIES 

0.32513E  03  0.100676  03  0.63900E  02  0.10067E  03  0.58793E  02  0.22595E  02  0.63900E  02  0,22595E  02  0.10924E  03 


RECUCEO  STIFFNESS  MATRIX 

0.95252E  06  0.294946  06  -0.49233E-10  0.29494E  06  0.17225E  06  -0.  1928IE~10  -0.49233E-  10  -0.19281E-10  0.32003E  06 


01 SP.  DISPLACEMENT  FORCE  RELATIONS  FORCES 


UX 

0.2235E-05 

-0.  3827E-05 

0.11326-21 

-0 .3456  6-12 

0,6409E- 12 

-0.7792E-12 

NX 

VX 

-0.3827E-05 

0.  1236E-04 

0.1559E-21 

Q.6225E-12 

-0. 2756E-11 

-0.36806-12 

NY 

VXPUY 

0.1132E-21 

0.  1559E-2 1 

0.  3125E-05 

-0,a006E-12 

-0.5443E-12 

-0.2920E-12 

NXY 

WXX 

-0.3456E-12 

0.6225E-12 

-0.8006E-12 

O.603OE-O2 

-0.1104E-01 

-0.1712E-02 

MX 

WYY 

C. 64096-12 

-0. 2756E-11 

-0.5443E-12 

-0.1104E-01 

0.3631E-01 

-0.1055E-02 

MY 

WXY 

-0.7792E-12 

-0. 3680E-12 

-0.2920E-12 

-0 .1712E-02 

-0.1055E-C2 

0.1037E-01 

MXY 

DISP. 


DISPLACEMENT  FORCE  RELATIONS 


FORCES 


0,i211E-Ol 

0.2235E-05 

-0.3827E-05 

0 .1132E-2 1 

-0.3456E-12 

0. 64096-12 

-0.7792E-12 

0.4833E 

34 

-0.2273E-01 

-0.382  7E-05 

0.1236E-04 

0.1559  E-21 

0.6225E-12 

-0.2756E-  11 

-0.3680E-12 

-0.3427E 

03 

-0.4003E-  10 

0.1132E-21 

0.1559E-21 

0 .3125E-05 

-0.8006E-12 

-0.5443E-12 

-0.2920E-12 

0. 

0.3415E  OC 

-0.3456E-12 

0.6225E-12 

-0.a006E-l2 

0.6833  E-02 

-0.1104E-01 

-0.1712E-02 

O.SOpOE 

02 

-0.5518E  00 

0.6409E-12 

-0.2756E-11 

-0 .54436-12 

-0.1104E-01 

0.3631E-01 

-0.1055E-02 

-0.2980E 

-07 

-0,85616-01 

-0. 7792E-12 

-0.3680E-12 

-0.2920E-12 

-0.1712E-02 

-0.1055E-C2 

0.1037E-01 

0.2980E 

-07 

FOR  THIS  CASE  NBS<  X,Y,XY-MI  IS  5000.  0.  0. 

FOR  THIS  CASE  MBS(  X,  Y  ,XY-M)  IS  50.  0.  0. 

LAYER  PROPERTIES,  ROWS-PROPE  RTY  ,  COLUMNS-L  A  YER 


1 

KV 

0. 

0. 

0. 

0  . 

0. 

0. 

0. 

0. 

2 

KF 

0.  5000E 

00 

0.  5000E 

00 

0.50306 

00 

0.50006 

00 

0.5000E 

00 

0.50006 

00 

0.5000E 

33 

0.53306 

33 

3 

KFB 

0,  5000E 

00 

0.  5000E 

00 

0.50006 

00 

0 .500  06 

00 

0.5000E 

00 

0.5000E 

00 

0.5000E 

03 

0.5300E 

03 

4 

KM 

0.5000E 

00 

0. 5000E 

00 

0.50006 

00 

0 .50006 

00 

0.5000E 

00 

0.5000E 

00 

0.5000E 

03 

0.5303E 

33 

5 

KMB 

0.5000E 

00 

0. 5000E 

00 

0.50006 

00 

0.500  06 

00 

0.5000E 

00 

0.50COE 

00 

0.5000E 

30 

O.SOOOE 

33 

6 

RHOL 

0.  5165E- 

01 

0.5165E- 

-01 

0.5165E- 

-01 

0 ,51656- 

-0  1 

0.5165E- 

01 

0.5165E- 

01 

O.5165E-01 

0.51656 

-31 

7 

TL 

0. 8000E- 

02 

0.  8000E-02 

0.8003E- 

-02 

0.80006- 

-0  2 

0.8000E- 

02 

0.8000E- 

02 

0.  SOnoE 

-02 

0.93036 

-3  2 

8 

DELTA 

0.6586E- 

04 

0.6586E-04 

0.65866- 

-04 

0.65866- 

-04 

0.6586E- 

04 

0.6586E- 

04 

0.6586E 

-04 

0.5586E 

-04 

9 

ILDC 

-O.OOOOE- 

19 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

10 

Z8 

0.4000E- 

02 

0.  1200E 

-01 

0.20006-01 

0.28006- 

-01 

0.3600E- 

01 

0.4400E- 

01 

0.  5200E 

-01 

0.6000E 

-31 

11 

ZGC 

-0.2800E- 

01 

-0.  2000E-01 

-0.12006- 

-01 

-0  ,40306- 

-0  2 

0,4000E- 

02 

0.1200E- 

01 

0.20006 

-01 

0.29006 

-01 

12 

THCS 

0. 

0. 

0. 

0  . 

0. 

0. 

0. 

0, 

13 

THLC 

0.5236E 

00 

-0.5236E 

00 

0.52366 

00 

-0.52366 

00 

-0.5236E 

00 

0.5236E 

00 

-0.5236E 

00 

0.52366 

03 

14 

THLS 

0.5236E 

00 

-0.  5236E 

00 

0.52366 

00 

-0 .52366 

00 

-0.5236E 

00 

0.52366 

00 

-0.52366 

30 

0.5236E 

33 

15 

sell 

0.2549E 

08 

0.25496 

08 

0.25496 

08 

0.25496 

08 

0.2549E 

OB 

0.254  9E 

08 

0.25496 

09 

0.25496 

08 

16 

SC12 

0.4118E 

06 

0.41186 

06 

0.41186 

06 

0,41186 

06 

0.4118E 

06 

0.4118E 

06 

0.41186 

06 

0,4118E 

06 

17 

SC  13 

0.4118E 

06 

0,41186 

06 

0.41186 

06 

0. 41186 

06 

0.41 106 

06 

0,4118E 

06 

0.41 186 

06 

0.41 18E 

06 

18 

SC  22 

0. 1160E 

07 

0. 11606 

07 

0.11606 

07 

0.1160E 

07 

0.1160E 

07 

0,1160E 

07 

0.11606 

07 

0.11606 

37 

19 

SC23 

0.4787E 

06 

0.47876 

06 

0.47876 

06 

0  .4787E 

06 

0.4787E 

06 

0.47876 

06 

0.4787E 

05 

0.4787E 

06 

79 


20 

SC  3  3 

0.11 60E 

07 

0.  1160E 

07 

0. 116D6 

07 

0 . 1160E 

07 

0.1160E 

07 

0. 11606 

07 

0.  U60E 

07 

0.1160E 

07 

21 

SC44 

0.3405E 

06 

0.  3405E 

06 

0.3405E 

06 

0,340  5E 

06 

0.3405E 

06 

0. 3405E 

06 

0.  3405E 

06 

0,34056 

06 

22 

SC  55 

0.3405E 

06 

0.  3405E 

06 

0.3405  £ 

06 

0 .3405E 

06 

0.3405E 

06 

0.340  5E 

06 

0.  3405E 

OS 

0.3405E 

06 

23 

SC66 

0.6339E 

06 

0. 6339E 

06 

0.6339E 

06 

0 .63396 

06 

0.6339E 

06 

0.63396 

06 

n.6339E 

06 

0.S339E 

06 

24 

CTEl  1 

-0.6I38E 

-07 

-0.61  38E 

-0  7 

-0.6  L33  E- 

-0  7 

-0.61386- 

-07 

-0.6138E- 

07 

-0.6138E- 

07 

-0.5138E 

-07 

-O.Sl 3BE 

-0  7 

25 

CTE22 

0.2334E 

-04 

0. 2334E 

-04 

0.2334E- 

-04 

0,2334E 

-04 

0.2334E- 

C4 

0.2334E- 

04 

0.  2334E 

-04 

0.2334E 

-04 

26 

CTe33 

0. 2334E 

-04 

0.  2334E-04 

0.2334E- 

-04 

0 . 2334E 

-0  4 

0.2334E- 

04 

0.2334E- 

04 

0.  23346 

-04 

0.23346 

-04 

27 

HKll 

0. 2906E 

03 

0.  2906E 

03 

0.2906E 

03 

0.2906E 

03 

0.29066 

03 

0. 29066 

03 

0.  29  06E 

03 

0.2906E 

03 

28 

HK2  2 

0.3715E 

01 

0.3  715E 

01 

0.3715E 

01 

C .371 56 

0  1 

0.3715E 

01 

0.3715E 

01 

0.37156 

01 

0.3715E 

01 

29 

h:<33 

0.37i5E 

01 

0.  37156 

01 

0.37L5E 

01 

0 .37156 

01 

0.37156 

01 

0.3715E 

01 

0.3715E 

01 

0.37156 

01 

30 

HCL 

0, 2043E 

00 

0.2043E 

00 

0.2043E 

00 

C  .20436 

00 

0.2043E 

00 

0.2043E 

00 

0. 2043E 

00 

0.20436 

00 

31 

EL  11 

0. 2523E 

08 

0.2528E 

08 

0.2528E 

08 

0 .2528E 

08 

0 .2528E 

03 

0.25286 

08 

n.  252  86 

OB 

0.25286 

08 

32 

EL  22 

0.9557E 

06 

0.9597E 

06 

0.9597E 

06 

0.9597E 

06 

0.9597E 

06 

0.9  59  76 

06 

0.  95976 

06 

0.9597E 

06 

33 

EL  33 

0.9597E 

06 

0. 9597E 

06 

0.9597E 

06 

0 .9597E 

06 

0.9597E 

06 

0.9597E 

06 

0.  95976 

06 

0.9597E 

06 

34 

a23 

0, 34C5E 

06 

0.3405E 

06 

0.3435  E 

06 

0,34D5E 

06 

0.3405E 

06 

0. 34056 

06 

0.  3405E 

06 

0.3405E 

06 

35 

GL  13 

0.6339E 

06 

0.6339E 

06 

0,6339E 

06 

0 .63396 

06 

0  .6339E 

06 

0.6339E 

06 

0. 63396 

06 

0.6339E 

06 

36 

a  12 

0. 6339E 

06 

0.6339E 

06 

0.6339E 

06 

0 ,63396 

06 

0  .6339E 

06 

0.6339E 

06 

0.63396 

06 

0.63396 

06 

37 

NUL  1  2 

0.  2514E 

00 

0. 2514E 

00 

0.2514E 

00 

0.251 46 

00 

0.2514E 

00 

0. 251  46 

00 

0. 25146 

00 

0.25146 

00 

38 

NUL21 

0.9541E 

-02 

0.  9541E 

-02 

0.9541 E- 

-02 

0.954  IE 

-02 

0.9541E- 

02 

0.95416- 

02 

0. 9541E 

-02 

0.9541E 

-0  2 

39 

NUL  1  3 

0.2514E 

00 

0.2514E 

00 

0.2514E 

00 

0.2514E 

00 

0.2514E 

00 

0.251  46 

00 

0.  2514E 

00 

0.25146 

00 

40 

NUL  3  1 

0.9541E 

-02 

0. 9541E 

-02 

0.9541  E- 

-02 

0  .954  IE 

-0  2 

0,95416- 

02 

0.9541E- 

02 

0.9541E 

-02 

0.9541E 

-0  2 

41 

NUL  2  3 

0.40  94E 

00 

0.4094E 

00 

0.4094E 

00 

0  .409  46 

00 

0.4094E 

00 

0.4094E 

00 

0.40946 

00 

0.40946 

00 

42 

NUL  3  2 

0. 4C94E 

00 

0.4094E 

00 

0.4094E 

00 

0 .4094E 

00 

0.40946 

00 

0.409  46 

00 

0.40946 

00 

0.4094E 

00 

43 

SMFK22 

0. 1522E 

01 

0.1522E 

01 

0.1522E 

01 

0 .152  2£ 

01 

0.1522E 

01 

0. 1522F 

01 

0,15226 

01 

0.1522E 

01 

44 

SMFD22 

0 .  1  9 1  8E 

01 

0.  1  91  8E 

01 

0. 1913E 

01 

0 . 191 8E 

01 

0.1918E 

01 

0. 1918E 

01 

0.  1  9186 

01 

0.19186 

01 

45 

SMFS22 

0. 1344E 

01 

0.  1342E 

01 

0.  1408E 

01 

0. 1400E 

0  L 

0.1399E 

01 

0. 1397E 

01 

0.  1399E 

01 

0,13956 

01 

46 

SMFC22 

-O.OOOOE 

-19 

-0.  OOOOE 

-19 

-0.0  00  3  E 

-19 

-O.OOOOE 

-19 

-O.OOOOE- 

19 

-O.OOOOE- 

1  9 

-O.OOOOE 

-19 

-O.OOOOE 

-19 

47 

SMFS12 

0. 3024E 

01 

0.3024E 

01 

0.3024E 

01 

0.30246 

01 

n,3024E 

01 

0.  302  4E 

01 

0.  3024E 

01 

0.30246 

01 

48 

SMFS23 

0.  1396E 

01 

0.13966 

01 

0. 1396E 

01 

0. 13966 

01 

0,1396E 

01 

0.13966 

01 

0.  1396E 

01 

0,13966 

01 

49 

ILMFC 

-O.OOCOE 

-19 

0. 70606 

02 

0.7063E 

02 

0 .7060F 

02 

0.7060E 

02 

0. 7060E 

02 

0.7060E 

02 

0.7060F 

02 

50 

TEMPO 

-0.  30  00E 

03 

-0.  3000E 

03 

-9.3300E 

03 

-0. 30006 

03 

-0.3000E 

03 

-0. 30006 

03 

-0.  3000E 

03 

-0.30006 

03 

51 

L  SCI  IT 

0.9676E 

05 

0.  9676E 

05 

0.9676E 

05 

0  .9676E 

05 

0.96766 

05 

0.9676E 

05 

0.  9676E 

05 

0.9676E 

05 

52 

LSCl  IC 

0. 5333E 

05 

0. 5333E 

05 

0.5333E 

05 

0 .533  3E 

05 

0.5333E 

05 

0.5333E 

05 

0.5333E 

05 

0.5333E 

05 

53 

L  SCI  ID 

0.6578E 

05 

0.6578E 

05 

0.6578  E 

05 

0.65786 

05 

0.6578E 

05 

0.6578E 

05 

0.65T8E 

05 

0.S57RF 

05 

54 

LSC22r 

0. 3  7i4E 

04 

0.3718E 

04 

0, 3544E 

04 

0 .35646 

04 

0.3566E 

04 

0.3572E 

04 

0.3568E 

0  4 

0.3576E 

04 

55 
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0.5722E-05 

VXPUY 

-0.6456F-03 

MX 

C.3662E-03 

0. L907E-05 

0.10916-10 

D. 38726  03 

0.52  83E  01 

-0.  1370E-04 

WXX 

-0.745  IE -0  3 

MY 

0.2384E-05 

0.4387E-04 

0.4657E-09 

0.5283E  01 

0.18756  03 

-0. 5963E-03 

WYY 

-0.7451 E-07 

MXY 

C.1091E-10 

0.3492E-09 

0.5722E-05 

-0 . 1  370  E-04 

-0. 5963E-03 

0.1385E  02 

WXY 

0.8527E-I3 

RECUCEU  BENJING  RcGIUITIES 

O.30723E  03  J.52B25E  01  -0.13o98E-04  0.52825E  01  0.18748E  03  -0.5963  3E-03  -0.  13698E-04  -0.  59633E-03  0.  13B47E  02 


KECUCfcO  STIFFNE  SS  MATRI  X 

0.84135E  Oft  0.l547ft£  05  -0,ft421lE-01  0.i5476E  05  0.84185E  06  -0.27953E  01  -0.6421  1E-C1  -0.27953E  01  0.405S8E  05 


UISP  . 


DISPLACEMENT  FORCE  RELATIONS 


FORCES 


C.ll  38E-05 

-0.21846-07 

0.3756E-12 

-0.1124E-11 

0.2167E- 13 

-0.7195E-18 

NX 

-0.21  846-07 

C.  11  38E-05 

0.8184E-10 

0.1860  6-13 

-0.2783E-12 

-0.7574E-16 

NY 

C. 37566-12 

0. 8184E-10 

0.2465E-04 

-0.2750E-18 

-0. 1128E-15 

-0. 1019E-10 

NXY 

VXPUY 


*yxx 

-0.11246-11 

0.  17146-1  3 

-0 .584  76-18 

0.2583  6-02 

-0.72796-04 

-0. 57926-09 

MX 

WYY 

0. 24696-13 

-0.27836-12 

-0.97466-16 

-0 .7  279  E-04 

0.53366-02 

0.22976-06 

MY 

»<XY 

-0.40576-18 

-0. 85736-16 

-0.10196-10 

-0.5792  6-09 

0.22976-06 

0.72226-01 

MXY 

OISP  . 


DISPLACEMENT  FGRC E  R EL AT  IONS 


FORCES 


0 .5tA4E-U2 
-O.AC62E-03 
-0.34976-07 

0.1292E  00 
-0  .36406-02 
-0.2896E-07 


0. 11 386-05 
-0.21  dAE-07 
0.375bE-l2 

-0.11246-11 
0.2469E-13 
-0.40576-1  8 


-0,2184E-07 

0.11886-05 

0.8184E-10 

0.1714E-13 

-0.2783E-12 

-0.8573E-16 


0  .3T56E-12 
0.8184E-10 
0.2465E-04 

-0.5847E-18 
-0  .9746E-16 
-0  .l019E-ia 


-0.1  124E-11 
0,18606-13 
-0.2750E-  IS 

0.25336-02 
-C.7279E-04 
-0 .57926-09 


0.21676-  13 
-0.27836-12 
-O.l  128E-  15 

-0. 72  7 96-  04 
C  .53366-02 
C.2297E-  C6 


-0.  7195F-13 
-o. 7574F-1 6 
-0.  1019E-1  0 

-0. 5792F-09 
0.  2297E-06 
n. 72226-01 


0,47456  04 
-0.2546E  03 
-0.6456E-03 

0.5000E  ^2 
-0.74516-^7 
0.3527E-13 


FOR  THIS  CASE  N  B  S  (  X  ,  Y  ,  X  Y- M)  IS  5000. 

FOR  this  case  M8S(  X,  Y  ,XY-M)  IS  50. 

layer  properties,  RQWS-PROPERTY  ,  COLUMNS-LAYER 


1 

KV 

2 

KF 

3 

KFB 

4 

KM 

5 

KMB 

6 

RHOL 

7 

TL 

8 

DELTA 

9 

ILOC 

10 

Z3 

1  1 

ZGC 

L2 

THCS 

13 

THLC 

14 

THLS 

15 

SC  11 

16 

sc  12 

17 

SC  13 

18 

SC  22 

19 

SC  23 

20 

SC  33 

21 

SC44 

22 

SC  55 

23 

SC66 

24 

CTEl  1 

25 

CT622 

26 

CT633 

27 

HK  11 

28 

HK2  2 

29 

HK3  3 

30 

HCL 

31 

6L  11 

32 

EL  22 

33 

EL  33 

34 

GL  23 

35 

GL  13 

36 

GL  12 

37 

NUL  1  2 

38 

NUL2  1 

39 

NUL  1  3 

40 

NUL  31 

41 

NUL  2  3 

42 

NUL  32 

43 

SMFK22 

44 

SMF0  22 

45 

S.MFS22 

46 

SMFC22 

47 

SMFS12 

48 

SMFS23 

49 

ILMFC 

50 

TEMPO 

51 

LSCllT 

52 

LSCliC 

53 

LSCl  ID 

54 

LSC22T 

55 

LSC22C 

56 

LSC12 

57 

LSC23 

53 

LSCC23 

59 

LSCC  13 

60 

LSCOF 

0. 

0.  5000E  00 
0.  500J6  00 
O.SOCJh  00 
0.50  COE  00 
0.51656-01 
0.  80  006-02 
0  •  6  5  8  66  -  0  4 
-0.00006-19 
0.  40  006-02 
-0.2  8006-01 
0. 

0. 

0, 

0.25496  03 
0.41186  06 
0.41 186  06 
0.  1 1606  07 
0.47876  06 
0.11606  07 
0.34056  06 
0.  34056  06 
0.63396  06 
-0.61386-07 
0,23346-04 
0.  2  6  346-04 
0.29C6E  03 
0.3715E  01 
0.3  715E  01 
0.2C43E  00 
0.2528E  08 
0.9597E  06 
0,95976  06 
0.  34056  06 
0.6339E  06 
0.63396  06 
0,25146  00 
0.95416-02 
0.2514E  00 
0.95416-02 
0.40  946  00 
0.4C94E  00 
0,1522E  01 
0.  191  36  01 
0.  13696  01 
-0. OC 006-19 
0.30246  01 
0.  13966  01 
-0,00006-19 
-0.30006  03 
0.96766  05 
0.53336  05 
0.6578E  05 
0.3644E  04 
0.17526  05 
0.254  76  04 
0.18666  04 
-O.OCOOE-19 
-0. OOOOE-19 
-O.OOOOE-19 


0. 

0.  50006  00 
0.5000E  00 
0.50006  00 
0,50006  00 
0.51656-01 
0. 30006-02 
0.65866-04 
0. 

0.  120UE-01 
-0.  20006-01 

0. 

0.  15716  01 
0.1571E  01 
0.25496  08 
0.41  186  06 
0.4118E  06 
0. 11606  07 
0.4787E  06 
0.1160F  07 
0.  34  05  6  06 
0.34056  06 
G. 63396  06 
-0.61386-07 
0.23346-04 
0. 23346-04 
0.2906E  03 
0.37156  01 
0.37156  01 
0.20436  0.1 
0.25286  08 
0.9597E  06 
0.  95  976  06 
0.34056  06 
0.63396  06 
0.63396  06 
0.25146  00 
0.  95416-02 
0.25146  00 
0.  95416-02 
0.40946  00 
0.40946  00 
0.15226  01 
0.1  9136  01 
0.  1332E  01 
-0.  00006-19 
0.30246  01 
0. 13966  01 
0.70606  02 
-0.  30006  03 
0.  96  76  6  05 
0.  53336  05 
0.65786  05 
0,36106  04 
0.17366  05 
0,25476  04 
0.18666  04 
-0. 00006-19 
-0.00006-19 
0. 10526-01 


0. 

0.50006  00 
0.5000E  00 
0.50036  00 
0.5000E  00 
0.51656-01 
0.80006-02 
0.65866-04 
0, 

0.20006-01 

-0.1200E-OL 


0. 

0. 

0. 


0.25496 

08 

0.41186 

06 

0.41136 

06 

0. 11606 

07 

0,47876 

06 

0.11606 

07 

0.34056 

06 

0.3405E 

06 

0.63396 

06 

-0.  t>1336- 

■0  7 

0.23346- 

-04 

0.23346- 

-04 

0. 2906  E 

03 

0. 37156 

01 

0.37156 

01 

0.20436 

00 

0.25286 

08 

0.9597E 

06 

0.95976 

06 

0.34)56 

06 

0.63396 

06 

0.63396 

06 

0.25146 

00 

0.95416- 

-02 

0.25146 

00 

0.95416- 

-02 

0.40946 

00 

0.40946 

00 

0.15226 

01 

0.19136 

01 

0.13576 

01 

-O.OOOOE- 

-19 

0.30246 

01 

0.13966 

01 

0,70606 

02 

-0.3JOOE 

03 

0.96766 

05 

0. 53336 

05 

0.65786 

05 

0.36796 

04 

0.17696 

05 

0.25476 

04 

0.18666 

04 

-0.00006 

-19 

-0.00006 

-19 

0.10526 

-01 

0  . 


0.50006 

no 

0 .50006 

0  0 

0 .50006 

0  0 

0 .50006 

CO 

0 .51656- 

0  1 

0.80006- 

■0  2 

0 .653  66- 

■0  4 

0  . 

0 .28006- 

■0  1 

-0  .40006- 

■0  2 

0  . 

0,157  IE 

01 

0.157  IE 

01 

0 .25496 

08 

0.41186 

06 

0.41186 

06 

0.11606 

07 

0  .47876 

06 

0 .11606 

07 

0.34056 

06 

0.340  56 

06 

0  .63396 

06 

-0.61386- 

-0  7 

0 .23346- 

-04 

0.23346- 

-0  4 

0.290  66 

03 

0  .37156 

0  1 

0 .371 56 

0  1 

0 .20436 

00 

0  .252  86 

08 

0.959  76 

06 

0.95976 

Oo 

0 .340  56 

06 

0 .63396 

06 

0 .63396 

06 

P.2514E 

00 

0.954  16 

-0  2 

0.251  46 

00 

0  .954  16 

-02 

0  .40946 

00 

0 .40946 

00 

0 .152  2F 

01 

0.19186 

01 

0 .13826 

0  1 

-O.OOOOE 

-19 

0.30246 

01 

0  .  13966 

01 

0.70606 

02 

-0.30006 

0  3 

D ,96766 

05 

0 .53336 

05 

0.6578E 

05 

0.36106 

04 

0. 17366 

05 

0.25476  04 
0.18566  04 
-0.00006-19 
-0.00006-19 
0 .10526-01 


0. 

0.50C0E  00 
0.50006  00 
0.50  COE  00 
0.50006  00 
0. 51656-01 
0 .800nE-U2 
0.65866-04 
0. 

O.36COE-01 
0  .40  006-0? 

0. 

0.15716  Cl 
0,15716  01 
0.?549E  08 
0.4118E  C6 
0.41186  06 
0.11606  07 
0.47876  06 
C,1160E  07 
0.34056  06 
C.34C56  06 
0.63396  C6 
-0.61386-07 
0.23346-04 
0.23346-04 
0.2Q06E  03 
0.37156  01 
0.3715E  01 
0.20436  00 
0.25286  08 
0.95976  06 
0.95976  06 
0.34056  06 
0,63396  06 
0.63396  06 
0.25146  00 
0. 95416-02 
0.25146  00 
0.95416-02 
0,40946  00 
0.4094E  00 
0.15226  01 
0.19186  01 
0.13826  Cl 
-0  .00006-19 
0.30246  01 
0.13966  01 
0.70606  02 
-0.30006  03 
0.96766  05 
0.53336  05 
0.65786  05 
0.361  IE  04 
0.17366  05 
0,25476  04 
0.18666  04 
-O.OOOOE-19 
-0.00006-19 
0,10526-01 


0  . 

0.500  06  0  0 
0.50006  00 
0 . 50C  06  ^ 
0.50006  00 
0. 51656-01 
0.80006-02 
0.658  66-04 
0  • 

0.44006-01 
0. 1 2006-01 

0. 

0. 

0. 

0.25496  ^8 
0.41186  0  6 
0.41  1  86  0  6 
0.  11606  07 
0.47876  06 
0.11606  0  7 
0.34056  06 
0,34056  06 
0.  6339E  <^6 
-0.61386-07 
0, 23346-04 
0. 23346-04 
0.  29066  03 
C.3715E  01 
0,37156  01 
0.20436  00 
0,252  86  C8 
0.9597E  06 
0.9597E  06 
0.34056  06 
0.63396  06 
0,63396  06 
0.2514E  00 
0.95416-02 
0.251  46  00 
0.9  54  15-02 
0,409  46  00 
0,40946  OP 
0.15226  01 
0.191  86  01 
0.13376  01 
-0,00006-19 
0.30246  01 
0.13966  01 
0.7060E  02 
-0.3000E  0? 
0.9676E  05 
0.5333E  05 
0.65786  05 
0.37316  04 
0.17946  05 
n.2547E  04 

0.18666  04 
-O.OOOOE-  19 
-0,00  006-1  9 
0. 10526-01 


0. 

0.5000E  00 
0.50006  00 
0.50006  00 
0.50006  00 
0.  51 656-01 
0.  80006-0? 
0.  65866-04 
0, 

0,  5?00f^-91 
0. 20096-01 

0. 

0.  1  5716  01 
9,  1  5716  01 
0.25496  08 
0.  41  1  86  05 
0.41  1  86  06 
0.11606  07 
0.47376  06 
9.1169E  07 
0.34056  05 
0.34056  05 
0.63396  06 

-  0,61  386-07 
0.  233^6-04 
0. 23346-04 
0,29065  03 
9.  37156  01 
0.37155  01 
0,20436  00 
0.25286  03 
0.  95976  06 
0.  95976  05 
0,34056  06 
0.63396  05 
0.63396  06 
0.2514F  00 
0.  95415-02 
0.25145  DO 
0.  95415-02 
0.4094E  no 
0.40946  00 
9.15226  01 
0.  1  91  85  01 
0.  13826  01 

-0,00006-19 
0.  30245  01 
0.13965  01 
0.79606  02 

-  0.  3  0006  03 
0,96766  05 
0.53336  05 
0.65786  05 
0.3610E  04 
9.  1  7366  05 
0.  25476  04 
0.18666  04 

-  9,  00006-19 
-0.00006-19 

9.1  0526-01 


0. 

0.50006  DO 
•0.5  00  0E  OO 
0.5000F  00 
0.5009E  9n 
0.51  55E-01 
0.3000F-02 
0.55866-04 
0. 

0.5000E-0 1 
0.28906-01 

C  , 

0. 

0. 

0.2549E  08 
0.41186  06 
0.41 186  06 
0.11606  07 
0.47875  06 
0.11696  07 
0.34056  06 
0.34  05  6  9f, 
0.63396  06 
-0.513  86-0  7 
0.23346-94 
0.23346-94 
0.29066  03 
0.37156  01 
9.37156  01 
9.20436  00 
0.2528»='  08 
0.95976  06 
0.95976  06 
0.3405E  06 
0.6339E  06 
0.53395  06 
0,25146  00 
0.95416-02 
0.25146  00 
0.95416-92 
0.40946  00 
0.40946  00 
0.15226  01 
0.19136  01 
0.13246  01 
-0.00006-1  9 

0.3024E  01 
9.1396F  01 
0.7060F  02 
-0.30006  03 
0.95766  05 
0.53336  05 
0.55786  05 
0.37586  04 
0.19126  05 
0.25476  04 
0.1866E  04 
-0. 00906-19 
-0.99006-19 
0.10526-01 


86 


0.1371E  01 
•0.15UE  01 
-O.OOOOE-19 
0.2027E~02 
■0.3043E-03 
■0.3416E-07 
0.5131E  05 
0,6929E  04 
0.2165E-01 
O.OOOOE-19 
0.1267E  01 


0.1371E  01 
-0.6581E  01 
O.IOOOE  01 
-0.3334E-03 
0.3061E-02 
0. 5927E-07 
■0.  7299E  04 
0.95956  04 
0.3757E~0l 
■0.  3564E-07 
•0.2659E  01 


0.1371E  01 
-0.1208E  01 
O.IOOOE  01 
0.4094E-02 
-•0.3625E-03 
-0.3462E-07 
0.1057E  06 
0.7372E  04 
-0.2195E-01 
-0.1009E-07 
-0.2548E  01 


0.1371E  01 
>0  .1002E  02 
O.IOOOE  01 
-0.3917E-03 
0  ,5128E-02 
0.7531E-07 
-0.8739E  04 
0.1157E  05 
0.4774E-01 
-0.418  6E-07 
-0.3783E  01 


0.1371E  01 
-0.1185E  02 
O.IOOOE  01 
-0.4208E-03 
0.6161E-02 
0.8333E-07 
-0.8199E  04 
0.1256E  05 
0.52826-01 
0.4249E-08 
-0.4372E  01 


0.1371E  01 
■0.1650E  01 
O.IOOOE  01 
0.7194E-02 
■0.4499E-03 
-0.35326-07 
0.1851E  06 
0.8038E  04 
0.2239E-01 
0.8689E-07 
0.6462E  01 


0.1371E  01 
-O.I638e  02 
O.IOOOE  01 
-0.4790E-03 
0.  822BE-02 
0. 9936E-07 
-0.1102E  05 
0.1453E  05 
0.62  98E-01 
-0.  5031E-07 
-0.57040  01 


0,1371E  01 
-0.2525E  01 
O.IOOOE  01 
0.9261E-02 
-0.5082E-33 
-0,3578E-07 
0.2379E  06 
0.8482E  04 
-0.2263E-0 1 
-0.1S35F-06 
-0.1041E  02 


rH0RNEL-50/£P0XY 


NL  ,NPL,NPC,.NFPE  ,NLC 
8  71  54  1420 


Bidirectional  Composite  Residual  Stresses  Only 


PIE 

0.31416E  01 


O.IOOOOE 

0.57000E 

07 

06 

0,  iOOOOE 
0,36000E 

07 

00 

0,20000E 

0,36000E 

00 

00 

0.25000E 

0.36000E 

00 

00 

0.20000E 

0. 

00 

0.13000E 

0. 

07 

0.70000E 

0. 

06 

0.13000E  07 

0.57000E 

0.2i)000E 

0. 

01 

0.40000E 

0. 

01 

0.20000E 

0. 

01 

0. 

O.IOOOOE 

01 

0. 

O.IOOOOE 

01 

0. 

O.IOOOOE 

01 

0.23560E 

0. 

01 

0. 

0, 

0. 

0. 

0. 56000E- 

■05 

0.56000E' 

-05 

0.42800E- 

■04 

0.42a00£- 

-04 

0.5800C6 

0.22500E 

02 

00 

0.58000E 

02 

0.  17000E 

00 

0.12500E 

01 

0. 12500E 

01 

0.125COE 

Cl 

Q.25000E 

00 

0. 

0. 

O.IOOOOE 

01 

0.10500E 

01 

0.10500E 

01 

THCS,RHOF,RHCM,  DIAF 

0-  0.5900CE-01  0.44300E-01  0,260006-03 

KVL 

0.  0.  0.  0. 

KFL 


^C^.SCCOOE  00  0.500006  00  0.50000E  00  0.50000E  00  0.50000E  OO  0.50000E  00  0.5000CE  00  0.50000E  00 


0. 900008  02  0. 


0,900008  02  0.90000E  02  0. 


0.90000E  C2  0. 


PTEMP^°°^~°^  0.  80000E-02  0.80000E-02  0.80000E-02  0.80500E-02  0.805C0E-02  C.80000E-02 

-0.30000E  C3  -0.30000E  03  -0.30000E  03  -0.30000E  03  -0.30000E  03  -0.300006  03  -0.30000E  03  -0.30000E  03 

o“;So“  ?:j330o1o2  S.-fooSol??  o“:}oo°od“I  o“  O.-lSof  O.lOOOOEOl  0.V650CE-D1  O.lOOOOEOl 


0.23CGCE  06  0,21000E  05  C.20000E-01  0.50000E-01  0,45000E-01  0,45000E-01 


87 


3-0 


COMPOSITF  STRAIN  STRESS  TEMPERATURE  RRIATIONS  -  STRUCTURAL  AXES 
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Bidirectional  Composite  with  Bending-Stretching  Coupling;  Residual  Stresses  Only 
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3-D  COMPOSITE  STRAIN  STRESS  TEMPERATURE  RELATIONS  -  STRUCTURAL  AXES 
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COMPOSITE  PROPERTIES  -  VALID  ONLY  FOR  CONSTANT  TEMPERATURE  THROUGH  THICKNESS 
LINES  1  TO  3L  3-D  COMPOSITE  PROPERTIES  ABOUT  MATERIAL  AXES 
LINES  33  TO  54  2-D  COMPOSITE  PROPERTIES  ABOUT  STRUCTURAL  AXES 
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FOR  THIS  CASE  N  B  S(  X  ,  Y  ,  X  Y- M)  IS  0.  0.  0. 
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-O.OCOOE-19 
-0.1172E-02 
0.2473E-02 
-0.  5351E-07 
-0.  2737E  05 
0.882  7E  04 
-0.33926-01 
-0. 00006-  19 
-U.2624E  01 


0. 

0.5000E  00 
0.50006  00 
0. 5000E  00 
0,  50006  00 
0.  5165E-01 
0.  8000E-02 
0.65866-04 
0. 

0. 1200E-Oi 
-0.2000E-01 

0. 

0. 

0. 

0.2549E  03 
0.41136  06 
0,41186  06 
0. 1160E  07 
0.4787E  06 
0.1160E  07 
0.3405E  06 
0.3405E  06 
0.6339E  06 
-0,6138E-07 
0.2334E -04 
0.23346-04 
0.  2906E  03 
0.3715E  01 
0.3715E  01 
0,20436  00 
0.2528E  08 
0. 9597E  06 
0.  9597E  06 
0.  34  05E  06 
0.63396  06 
0.6339E  06 
0.25146  00 
0.  9541E-C2 
0.2514E  00 
0. 95416-02 
0.40946  00 
0,4094E  00 
0.1522E  01 
0.i913E  01 
0.1384F  01 
-0.00006-19 
0.3024E  01 
0.1396E  01 
0.7060E  02 
-0.  300QE  03 
0.9676E  05 
0.  5333E  05 
0.6578E  05 
0,3607E  04 
0.1734E  05 
0.2547E  04 
0.1866E  04 
-O.OOOOE-19 
-0.  OOOOE-19 
0.  1052E-0I 
0.1371E  01 
-0.5461E  01 

o.ioooe  01 

-0.6510E-03 
0,  1953E-02 
-0,49696-07 
-0.1480E  05 
0.  8452E  04 
-0.  3150E-01 
-0.47526-08 
-0.2190E  01 


0. 

0.5000E  00 
0.5000E  00 
0.5000E  00 
0.5000E  00 
0.5165F-01 
0.8000E-02 
0.6535E-04 
0. 

0.2000E-01 
-0. 1200E-01 

0. 

0, 

0. 

0.25496  08 
0,4118E  06 
0.4UBE  06 
0,1160E  07 
0.47876  06 
0.1160E  07 
0.34056  06 
0.3405E  06 
0,63396  06 
-0.613BE-07 
0.2334E-04 
0.2334E-04 
0.29066  03 
0.3715E  01 
0,3715E  01 
0,2043F  00 
0.25286  OR 
0.9597E  06 
0.9597E  06 
0.34056  06 
0.6339E  06 
0.6339E  06 
0.2514E  00 
0.9541 6-02 
0.2514E  00 
0.95416-02 
0.4094E  00 
O.4094E  00 
0.1522E  01 
0. 19136  01 
0.1332E  01 
-0.00006-19 
0,3024E  01 
0,13966  01 
0.70606  02 
-0.30006  03 
0.96766  05 
0.5333E  05 
0.6578E  05 
0.3612E  04 
0.17376  05 
0.2547E  04 
0.19666  04 
-0.00006-19 
-O.OOOOE-19 
0.I052E-01 
0.13716  01 
-0.4100E  01 
O.IOOOE  01 
-0. 1302E-03 
0.14326-02 
-0.45876-07 
-0.17296  04 
0.80776  04 
-0.2903  6-01 
-0.40906-09 
-0.18296  01 


0  , 

0.5000E  00 
0.50006  00 
0.5030E  00 
0.50006  00 
0 .51656-0  I 
0.80006-02 
0.658  66-04 
0  * 

0.28006-01 
-0 .4000E-02 

0  . 

0  . 

0  . 

0.2549F  08 
0.4118E  06 
0.411B6  06 
0.U6nE  07 
0.47376  06 
0.11606  0  7 
0.34056  06 
0.34056  06 
0.63396  06 
-0.61386-07 
0 .23346-04 
0 .233  46-0  4 
0.29066  03 
0.3715E  01 
0.3715E  01 
0.2043E  00 
0.25286  08 
0.9597E  06 
0.9597E  06 
0.3405E  06 
0.6339E  06 
0,6339E  06 
0.2514E  00 
0.954  16-02 
0.25146  00 
0.954  LF-OZ 
0.4C94E  00 
0.40946  00 
0.15226  01 
0.19196  01 
0 . 13796  0  1 
-0.00006-19 
0.30246  01 
C  .  13966  0  1 
0.70606  02 
-0.30006  03 
0.96766  05 
0.53336  05 
0.65786  05 
0.36196  04 
0.17406  05 
0.254  76  0  4 
0.186  66  04 
-0 .00006-19 
-0.00006-19 
0 . 10526-01 
0. 1371E  01 
-0.320  16  0  1 
O.IOOOE  01 
0.39056-03 
C.911  2E-03 
-0.42066-07 
0.11346  05 
0.77026  04 
-0  .26666-0 1 
0 . 1839E-08 
-0.  1540E  0  1 


0. 

0.5000E  00 
0.50COE  00 
0.50006  on 
0.50C0E  00 
0.51656-01 
Q. 80006-02 
0. 65866-04 
0. 

0.36006-01 
0  .40  00E-02 

0. 

0.15716  Cl 
0.15716  01 
0.25496  08 
0.41186  06 
0.41186  06 
0.11606  07 
0.47876  06 
0.11606  07 
0.3405E  06 
0.34056  C6 
0.63396  06 
-0.61386-07 
0.23346-04 
0.2334E-04 
0.29066  03 
0.37156  01 
0.37156  01 
0.20436  00 
0.25286  03 
0.95976  06 
0.95976  06 
0.3405E  06 
0.63396  06 
0.6339E  06 
0.2514E  00 
0.9541E-02 
0.2514E  00 
0,9541E-02 
0.40946  00 
0.4094E  00 
0.1522E  01 
0.1918E  01 
0.1379E  01 
-O.OOCOE-19 
0.30246  01 
'^.1  3966  01 
0.70606  02 
-0.30006  03 
0.96766  05 
0  .5333F  05 
0.65796  05 
0.3619E  04 
0.17406  05 
0.2547E  04 
0.18666  04 
-o.onooe-19 
-Q.OOCOE-  19 
0.10526-01 
0.1371F  Cl 
-0.32C1E  01 
0.10006  01 
0.39056-03 
0.91  12E-03 
0,42066-07 
0.1134E  05 
0.77C2E  04 
0.2666E-01 
0  .35816-08 
-0.1540E  01 


0. 

0.50CCE  00 
0.50006  CO 
0.500C6  00 
0.500CE  0  0 
0.51656-01 
0.80006-02 
0. 65866-04 
0. 

0.4400E-01 
0. 12006-^1 

0. 

C.1571E  01 
0.1571E  01 
0.2549E  ^^8 
0.41186  06 
0.41186  06 
0.1160E  07 
0.47876  06 
0,1  1606  07 
O.3405E  06 
0,3405E  06 
0.6339E  06 
-0,6138F-07 
0, 2334E-04 
0. 233'iE-04 
0  .  2  9  0  66  0  3 
0.371  56  01 
0.371  56  01 
0.2043E  00 
0.25286  08 
0.9597E  06 
0.95976  06 
0.340  56  06 
0.63396  06 
0.6339E  06 
0.2514E  00 
0.9541E-02 
0.  251  46  00 
0,95416-':'2 
0,40946  00 
0.409  4E  0  0 
0.15226  01 
0.191  86  91 
0.13826  01 
-0.00006-1  9 
0.302  46  01 
0.1396E  01 
0.70606  02 
-0.30CC6  03 
0.96766  05 
0.53336  05 
0.65786  05 
0,361  26  04 
0.17  3  76  0  5 
0.25476  04 
0.1866E  04 
-O.OOOOE-19 
-C.OOOOE-l  9 
0. 10526-01 
0.  13716  01 
-0.4100E  01 
O.IOOOE  01 
-0.13026-03 
0.14  3  26  -  02 
0.45876-07 
-0.17296  04 
0.00  7  76  04 
0.29C8E-01 
-0.30706-09 
-0.  13296  01 


0. 

0.50006  00 
0.50006  00 
0.50006  00 
0.  5000E  00 
0.51656-01 
0.  80006-02 
0. 6586E-04 
o. 

0.  52096-01 
0.  2000E-01 

0. 

0.  15716  01 
0.  15716  01 
0.25496  08 
0.411PE  9^, 
0.4118E  06 
0.  1  1606  07 
0. 4787E  06 
0,  1160E  07 
0.34056  06 
0,  34056  06 
0.  6339E  06 
-0.  61  396-07 
0.  2334E-04 
9.  2334E-04 
0.  29066  03 
0.3715E  01 
0.3715E  01 
0,20436  CO 
0.  252  86  08 
0.9597E  06 
0.  95976  06 
0.  34056  06 
^.63396  06 
0.  6339E  06 
0.  251  46  DO 
0. 954) E -02 

n.25i4E  no 

0.95416-02 
0.  4n94E  00 
0.40946  00 
0.15226  01 
0.  1  91  86  9L 
0.13846  01 
-o.nooo^-19 
0.  3024E  01 
0.1396E  01 
0,70606  02 
-0,  300CE  03 
0.  96766  05 
0.53336  05 
0.65786  05 
0.  3607E  04 
0.17346  05 
0.2547E  04 
0.186FE  04 
-O.OOOOE-19 
-0,00006-1  9 
0.1052E-01 
n.i371E  01 
-0.  5461E  01 
O.IOOOE  01 
-0.6510E-03 
0.  1953E-02 
0.4969E-07 
-0.14aCF  05 
0. 8452E  04 
0.  31506-01 
-0,2927E-08 
-0.219OE  01 


0. 

0.50006  00 
0.50O0E  oo 
0.50006  00 
O.5000E  00 
C. 51656-01 
0. 30006-02 
0.5586E-04 
0. 

0.6090r-01 

0.2500E-0) 

0, 

0.15716  01 
0.1571E  01 
0,25496  OB 
0.411FE  06 
0,41 186  06 
0.11696  07 
0,47876  06 
0,11606  07 
0.34956  06 
0.34056  06 
0.6339E  06 
-0.61 3PE-07 
0. 23346-04 
0.23346-04 
0,2906^  03 
0.37156  01 
0.3715F  01 
0.20436  OO 
0.25786  03 
0.9597E  06 
0,9597P  06 
0.3495F  96 
0.63396  06 
0.63306  06 
0.25146  00 
0.95416-02 
0,25146  0^ 
0.95416-02 
0.40946  0" 
0.4094F  CO 
0.15226  01 
0.19186  01 
0,13866  01 
-0.00006-19 
0.30246  01 
0.13966  01 
0.70606  02 
-0.3000F  08 
0.96766  05 
0,53336  05 
0.65786  05 
0.36016  04 
0,17316  05 
0,25476  04 
0.13666  04 
-O.OO^nf-iq 
-0,00006-19 
0.10526-01 
0.13716  01 
-O.7038E  01 
0.10006  01 
-0.11726-02 
0.2473E-02 
0.535 lE-07 
-0,27876  05 
0.88276  04 
0.33926-01 
-9,57936-08 
-0.26246  01 
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TABLE  I.  -  SUMMARY  OF  DETAILS  FOR  PREPARING  INPUT  DATA  CARDS.  (See  also  tables  II  to  IV. ) 
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TABLE  n.  -  MULTILAYERED  FIBER  COMPOSITE  ANALYSIS  INPUT  DATA  SAMPLE 


TH0RNEL-5O/EPOXY 

8  71  54  1420  1 


.  50  OOOE+0  8 

.  200  OOE+00 

. 57  OOOE  +  06 

.OOOOOE+00 

.  1 00  00  E+07 

. 1 30  OOE+07 

. 5  70  00E+  06 

. 0  00  OOE+06 

.  1 OOOOE  +0  7 

.  7000  0E+  0  6 

. 36000E+00 

. OOOOOE+00 

.  2  00  OOE  +0  0 

.  1  30  00E  +  07 

.  360  OOE  +  00 

.  25  OOOE+O  0 

.5  7  0  OOE  +0  6 

.  360  0  OE  +  00 

. 4000  OE  +  01 

.  00  0  OOE+00 

.  000  OOE+00 

.  10  OOOE+0 1 

.  20  0  00E+  01 

.  000  OOE  +  00 

.  OOOOOE  +  0  0 

.  1 OOOOE  +  01 

.4000  OE+0  1 

.  1000  OE+0 1 

. 0  000  0E+  00 

. OOOOOE+00 

.  20000  E+  01 

.  0  OOOOE+00 

.  00000  E  +0  0 

, 0  OOOOE  +  00 

. OOOOOE+00 

. OOOOOE+0  0 

.  1 000  OE  +  00 

.0  OOOOE+OO 

-  .  55  OOOE  -06 

. 5  6000E-05 

.  56  OOOE-05 

.  4  28  0  OE  -  04 

.  4280  OE  -04 

.  428  0  OE  -04 

.  58  00  OE  +  03 

. 12500E+01 

.  00  OOOE-f  00 

.  5  8  0  0  OE  +  02 

.  125  OOE  +  Ol 

.22  5  00  E+  00 

. 580  OOE  +  0  2 

.  25000E+00 

. 17000E+00 

,  OOOOOE+0  0 

.  125  OOE  +  Ol 

.OOOOOE+OO 

.  10  OOOE+O  1 

.  100  OOE  +  Ol 

. 10500E+01 

.  10500E  +  0  1 

.51416E+0  1 
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TABLE  n.  -  Continued.  MULTILAYERED  FIBER  COMPOSITE  ANALYSIS  INPUT  DATA  SAMPLE 


OOOOOE+00 

.  05  900E+00 

.  044  30E+00 

. 0  00  26E+  00 

OOOOOE+00 

0  0  00  0E+  0  0 

.  000  OOE  +  00 

.  00  OOOE  +  00 

. 00  0  00E  +  00 

.  .  0000  OE+00 

.OOOOOE  +  00 

. OOOOOE  +00 

50  OOOE+00 

5000  0E+  00 

. 5  00  00E+  00 

.50000E+00 

. 5  000  OE+0  0 

.  5000  OE  +  00 

. 50000E+00 

. 50000E+00 

OOOOOE  +00 

45000E+02 

..45000E+02 

. 45000E+02 

-  . 45000  E+  02 

.OOOOOE+00 

.  900  00E+  02 

.  9  OOOOE  +0  2 

0080 5 E+ 00 

008  05E  +  00 

.  00  805E  +  00 

.  00805E  +  00 

. 008  0  5E  +  00 

. 00  805E  +  0  0 

. 0  08  05E  +  00 

.  00  80  5  E+0  0 

3  0  00  OE  +  0  3 

30  OOOE  +  03 

~ .  300  OOE+03 

- . 30000E+03 

- . 30000E+03 

- . 30000E+03 

- . 30000E+03 

-  .  30000E  +  0  3 

TABLE  n. 


Concluded.  MULTILAYERED  FIBER  COMPOSITE  ANALYSIS  INPUT  DATA  SAMPLE 


830  OOE  -f  00 

16  5  OOE+02 

50000E+00 

10OOOE+O1 

.  1 00  OOE  +  00 

. 1  0000E+  0  1 

. 1 3300E+02 

.  2  600  0  E  +  00 

.  lOOOOE  +  01 

. 3 1900  E+05 

-  2  70  OOE+00 

. 04650E+00 

. lOOOOE+01 

.  1  70  0  0E+  00 

.  lOOOOE+0 1 

.  lOOOOE  +  01 

2300  OE  +06 

045  00  E  +  OO 

.  2  15  00E+  05 

.  02  000  E+00 

.050  00E+  00 

.  04  5  0  OE  +  00 

5000  OE  +  04 

. OOOOOE+00 

.  0  000  OE  +0  0 

lOOOOE+03 

.0000  OE+0  0 

.  00  OOOE  +  0  0 

00  OOOE  +  00 

OOOOOE+00 

. OOOOOE+00 

.00000  E+00 

. 0  0000E+  00 

.  00  OOOE  +  00 

TABLE  m.  -  INPUT  DATA  FOR  BORON/ALUMINUM  COMPOSITE 


boron/aluminum 

8  71  54 

1  1 

•60000E.08 

.60000E.08 

.60000E.08 

.200UUE.00 

.20000E.00 

•20000E.00 

•OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.10000E.08 

.10000E.08 

.10000E.08 

.3300OE.00 

.33000E.00 

.33000E.00 

.OOOOOE.OO 

•OOOOOE.OO 

.OOOOOE.OO 

.40000E.01 

.20000E.01 

.40000E.01 

.20000E.01 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.  lOOOOE.Ol 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.  lOOOOE.Ol 

.lOOOOE.Ol 

.lOOOOE.Ol 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

.28000E-05 

.28000E-05 

.28000E'-05 

.12900E-04 

.12900E--04 

. 12900E-04 

.22300E.02 

.22300E.02 

.22300E.02 

.31000E.00 

.10040E.04 

.10040E.04 

.10040E.04 

.23O00E.00 

•OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

.2250OE.OO 

. lOOOOE.Ol 
.31416E.01 

.lOOOOE.Ol 

.  lOOOOE.Ol 

.lOOOOE.Ol 

F 

F 

F 


F 


.OOOOOE.OO 

.08500E.OO 

.09a00E.00 

.004U0E.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE. 00 

0. OOOOOE. 00 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

O.5OOO0E.OO 

O.5OO0OE.OO 

0.50000E.OO 

o.5ooaoE.oo 

0.50000E.OO 

0.50000E.00 

0.50000E.00 

0.50000E.00 

.OOOOOE.OO 

.45000E.02 

--.45000E.02 

.90000E.02 

.90000E.02 

-.45000E.02 

.45000E.02 

.OOOOOE.OO 

.00500E.O0 

.00500E.OO 

.00500E.00 

.005UUE.00 

.OO5OOE.00 

.00500E.OO 

.00500E.00 

.OO5OOE.0O 

-0.90000E.03 

-0.90000E.03 

~0.90000E.03 

-0.90000E.03 

-0.90000E.03 

-0.90000E.03 

-0.90000E.03 

-0.90000E.03 

.56000E.00 

.  lOOOOE.Ol 

.313OOE.O0 

.46200E.00 

. 3OOOOE.00 

.29200E.01 

.lOOOOE.Ol 

.  lOOOOE.Ol 

•lOOOOE.Ol 

•lOOOOE.Ol 

.34300E.O0 

.83300E.01 

.52000E.05 

•lOOOOE.Ol 

•lOOOOE.Ol 

.lOOOOE.Ol 

.46000E.06 

.52000E.05 

.00520E.00 

.0O52OE.OO 

.O0905E.00 

.00905E.00 

.50000E.04 

.OOOOOE.OO 

.OOOOOE.OO 

• 10000E.03 

.OOOOOE.OO 

•OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE .00 

0. OOOOOE. 00 

0. OOOOOE. 00 

TABLE  IV.  -  INPUT  DATA  FDR  BORON/EPOXY  COMPOSITE 


BORON/EPOXY 
8  71  54 

•60000E.08 
.2OO00E.00 
.56000E.06 
•OOOOOE.OO 
.40000E.01 
•OOOOOE.OO 
•OOOOOE.OO 
•iOOOOE.Ol 
.28000E-05 
.32000E-04 
.22300E.02 
.17000E.01 
•OOOOOE.OO 
•IOOOOE.Ol 
.31416E.01 
F 
F 
F 
F 

•OOOOOE.OO 
.UUOOOE.OO 
0. OOOOOE.OO 
0.500O0E.O0 
0.50000E.OO 
.OOOOOE.OO 
-.45000E.02 
.00500E.00 
.0050OE.O0 
-0.30000E.03 
-0.30000E.03 
.84000E.OO 
.16500E.02 
.12200E.01 
.IOOOOE.Ol 
.46000E.06 
.05300E.00 
.50000E.04 
.iOOOOE.03 
0. OOOOOE.OO 
0. OOOOOE.OO 


1  1 
.60000E.08 
•OOOOOE.OO 
.56000E.06 
•OOOOOE.OO 
.20000E.01 
•OOOOOE.OO 
.OOOOOE.OO 
•IOOOOE.Ol 
.28000E-05 
.32000E-04 
.22300E.02 
.17000E.01 
.22500E.00 
.IOOOOE.Ol 


.08500E.00 
.OOOOOE.OO 
0. OOOOOE.OO 
0.50000E.00 
0.50000E.00 
.45000E.02 
.45000E.02 
.00500E.00 
.00500E.00 
-0.30000E.03 
-0.30000E.03 
•IOOOOE.Ol 
. IOOOOE.Ol 
.13300E.02 

.25000E.05 

.OOOOOE.OO 
.OOOOOE.OO 
0. OOOOOE.OO 


.60000E.08 
.OOOOOE.OO 
.35000E.00 
.OOOOOE.OO 
.40000E.01 
. IOOOOE.Ol 
.OOOOOE.OO 
.OOOOOE.OO 
.28000E-05 
.32000E-04 
.22300E.02 
.25000E.00 

. IOOOOE.Ol 


.04400E.00 
.OOOOOE.OO 
0. OOOOOE.OO 
0. 50000E.OO 
0.50000E.OO 
-.45000E .02 
.OOOOOE.OO 
.00500E.00 
.OO500E.OO 
-0.30000E.03 
~0.30000E.03 
.90000E.OO 
. IOOOOE.Ol 
.31900E.05 

.02700E.00 

•OOOOOE.OO 
.OOOOOE.OO 
0. OOOOOE.OO 


.20000E.OO 

•OOOOOE.OO 

.350O0E.0O 

.20000E.01 

.OOOOOE.OO 

.OOOOOE.OO 

.OOUUOE.OO 


.31000E.00 

.UOOOOE.OO 

. lOOOOt.Ol 


.00400E.00 

.OOOOOE.OO 

0.500U0E.00 

.90000E.02 

.O05O0E.00 

-0.30000E.03 

.15000E.01 

.120OOE.OO 

.IOOOOE.Ol 

.07000E.OO 


0. OOOOOE.OO 


.200O0E.00 

.56000E.06 

.3500OE.00 

.OOOOOE.OO 
.OOOOOE.OO 
. IOOOOE.Ol 
.OOOOOE.OO 


.17000E.01 

.OOOOOE.OO 


.OOOOOE.OO 

0.50000E.OO 

.90000E.02 

.00300E  .00 

-o.aooooE.os 

.10300E.01 

.IOOOOE.Ol 

.IOOOOE.Ol 

.O53OOE.O0 


0. OOOOOE.OO 


u.  u.  u.  u. 


TABLE  V.  -  INPUT  DATA  FOR  E-GLASS/EPOXY  COMPOSITE 


E-GLASS/EPOXY 
8  71  54 

.10600E.08 
.22000E.OO 
.50000E.06 
•OOOOOE.OO 
.40000E.01 
.OOOOOE.OO 
•OOOOOE.OO 
.lOOOOE.Ol 
.28000E-05 
.32000E-04 
.75,000E.0l 
.15000E.01 
•OOOOOE.OO 
•lOOOOE.Ol 
.31416E.01 


•OOOOOE.OO 
•OOOOOE.OO 
0. OOOOOE.OO 
0.50000E.OO 
0.50000E.00 
-.45000E.02 
-.45000E.02 
.00800E.00 
.00800E.OO 
-0.30000E.03 
-0.30000E.03 
.82000E.00 
.16500E.02 
•llOOOE.Ol 
.lOOOOE.Ol 
.36000E.06 
.03500E.00 
.50000E.04 
.10000E.03 
0. OOOOOE.OO 
0. OOOOOE.OO 


204  1 

.10600E.08 
.OOOOOE.OO 
.50000E.06 
.OOOOOE.OO 
.20000E.01 
•OOOOOE.OO 
•OOOOOE.OO 
. lOOOOE.Ol 
.28000E-05 
.32000E-04 
.75000E.01 
. 15000E.01 
.22500E.00 
.lOOOOE.Ol 


.09000E.OO 
.OOOOOE.OO 
0. OOOOOE.OO 
O.5OOOOE.O0 
0.50000E.OO 
.45000E.02 
.45000E.02 
.00800E.00 
.00800E.OO 
-0.30000E.03 
-0.30000E.03 
.lOOOOE.Ol 
.lOOOOE.Ol 
.13300E.02 

.25000E.05 

.OOOOOE.OO 
.OOOOOE.OO 
0. OOOOOE.OO 


•10600E.08 
•OOOOOE.OO 
.35000E.00 
•OOOOOE.OO 
.40000E.01 
• lOOOOE.Ol 
.OOOOOE.OO 
.OOOOOE.OO 
.28000E-05 
.32000E-04 
.75000E.01 
.25000E.00 

.90000E.00 


.04000E.OO 
.OOOOOE.OO 
0. OOOOOE.OO 
0.50000E.00 
0.5OO00E.00 
.OOOOOE.OO 
.OOOOOE.OO 
.00800E.00 
.00800E.OO 
-0.30000E.03 
-0.30000E.03 
.55OOOE.00 
. lOOOOE.Ol 
.31900E.05 

.02000E.00 

.OOOOOE.OO 
.OOOOOE.OO 
0. OOOOOE.OO 


.220a0E.00 

.OOOOOE.OO 

.35OO0E.OO 

.20000E.01 

.OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 


.17000E.OO 

•OOOOOE.OO 

.90000E.00 


.36000E-03 

•OOOOOE.OO 

0.50000E.OO 


.00800E.OO 

-0.30000E,03 

.8600OE.O0 

.33000E.00 

•lOOOOE.Ol 

.0500QE.00 


O.OOOOOE.OO 


.2200OE.0O 

.50000E.06 

.35O0OE.0O 

.OOOOOE.OO 

.OOOOOE.OO 

•lOOOOE.Ol 

.OOOOOE.OO 


.15000E.01 

•OOOOOE.OO 


.OOOOOE.OO 

0.50000E.00 


.00800E.00 

“0.30000E.03 

.82000E.00 

•lOOOOE.Ol 

•lOOOOE.Ol 

.03500E.OO 


O.OOOOOE.OO 


TABLE  VI.  «  INPUT  DATA  FOR  S-GLASS/EPOXY  COMH3SITE 


S-GLASS/EPOXY 
8  71  54 

•12400E,08 
.22000E,00 
o50000E.06 
•OOOOOE.OO 
«40000E-01 
•OOOOOEeOO 
ft0000CE«00 
olOOOOE.Ol 
«28000E-05 
e32000E~04 
.75000E<»01 
•17000Eo01 
*00000E«00 
.IOOOCE.0.1 
*31416Ec0i 
F 
F 
F 
F 

cOaoooE.oo 

^OUOOOE.OO 
0 ,OOOOOE«UO 
0.50O0OE.00 
0*50000E.OO 
*OOOOOE.OO 
-*45000Eo02 
*01290E.00 
.01290E*00 
-0o30000E  .U3 
-0  »30000E  oU3 
•lOOOOE.Ol 
. 16500E.02 
p 17000E»01 
• lOOOOE.Ol 
o360C0E -06 
«O3500E.00 
•50000E.04 
.10000E.03 
O.OOOOOE.OO 
OoOOCOOEoOO 


204  1 

«12400Eo08 
.OOOOOEeOO 
<,50000E-06 
•OOOOOEcOO 
o20000E-01 
•OOOOOEeOO 
«00000E*.00 
« lOOOOE.Ol 
•28000E-05 
o32000E'-04 
•75000Eo01 
•17000Ec01 
»22500E*00 
a iOOOOEoOi 


.^09000Ec00 
•OOOOOE-OO 
OeOOOOOE-00 
OcSOOOOE.OO 
Oc50000E«.00 
.45000E.02 
o450C0E«02 
o01290Ec00 
»01290E»00 
-0o30000E»03 
-0c30000Eo03 
-10000E«01 
. lOOOOEoOl 
* 13300E.02 

o25000E.05 

o00000E*00 

qOOOOGEcOO 

O.OOOOOE.OO 


. 12400Eo08 
c00000E«00 
.35000E-00 
o00O00E«00 
.40000E-01 
..  iOODOEaOl 
-OOOOOE-OO 
^.OOOOOEcOO 
«.28000E-05 
.32000E'“04 
c75000E*01 
o25000Ec00 

•90000E.OO 


. 04000E .00 
.OOOOOE.OO 
O.OOOOOE^OO 
0.50000E.00 
O.5O000E.00 
“.45000Eo02 
.OOOOOE.OO 
«01290Ec00 
.01290E«0G 
'-0-30000Ea03 
-0o30000E«03 
.66000E<>0C 
.  10000E.-0i 
. 31900E.05 

.02000E.OO 

•OOOOOE .00 
oOOOOOEcOO 
OoOOOOOEoOO 


.22OOOE0OO 
.OOOOOE. 00 
.35000E.00 

.20000E.01 
.OOOOOE. 00 
.OOOOOE. 00 
.OOOOOE. 00 


.  170OOE.OO 
.OOOOOE. 00 

.90000E.00 


.36000E-03 
.OOOOOE. 00 

0.50000E.00 

.9000CE.02 

.01290E.00 

~0.300UUE.03 

. 13000E.01 
.55000E.00 
. lOOOOE.Ol 

.O5OOOE0OO 


O.OOOOOE.OO 


.22OOOE0OO 

r50000E*06 

.35000E.00 

.OOOOOE. 00 
.OOOOOE. 00 
. lOOOOE.Ol 
.OOOOOE. 00 


.17000E.01 

cOOOOOE.OO 


.OOOOOEcOO 

O.5OOOOE0OO 

.90000E.02 

.O129OE.O0 

-0.3a000£.03 

.llOOOEoOi 

.lOOOOE.Ol 

.lOOOOE.Ol 

.O35OOE0OO 


O.OOOOOE.OO 


TABLE  Vn.  -  INPUT  DATA  FOR  THORNE L -2 5/ EPOXY  COMPOSITE 


THORNEL-25/EPOXY 

8  71  54  1440  1 

.25000E.08  .20000E.07 

.20000E.OU  .20000E.07 

.54600E,06  .54600E.06 

.OOOOOE.OO  .OOOOOE.OO 

♦40000E.01  .20000E.01 

•OOOOOE.OO  .OOOOOE.OO 

.OOOOOE.OO  .OOOOOE.OO 

.lOOOOE.Ol  .lOOOOE.Ol 

--.55000E-06  .56000E-05 

.42800E-04  .42a00E-04 

.58000E.03  .58000E.02 

.12500E.01  .12500E.01 

•OOOOOE.OO  .22500E.00 

.lOOOOE.Ol  .lOOOOE.Ol 

.31415E.01 
F 
F 
F 
F 

.OOOOOE.OO  .05200E.00 

•OOOOOE.OO  .OOOOOE.OO 

0. OOOOOE.OO  0. OOOOOE.OO 

0.5O000E.O0  0.50000E.00 

0.50000E.00  0.5OO0OE.OO 

•UOOOUE.OO  .45000E.02 

-.45000E.02  .45000E.02 

.01300E.00  .01300E.00 

.01300E.00  .01300E.00 

-0.30000E.03  -0.30000E.03 

-0.30000E.03  -0.30000E.03 

.lOOOOE.Ol  .lOOOOE.Ol 

.16500E.02  .lOOOOE.Ol 

.49000E.00  .13300E.02 

•lOOOOE.Ol 

.18000E.06  .21000E.05 

.04500E.00 

.50000E.04  .OOOOOE.OO 

.10000E.03  .OOOOOE.OO 

O.OOOOOE.OO  0. OOOOOE.OO 

0. OOOOOE.OO 


.20000E.07  •200U0E.OO 

.80000E.06  .20UO0E.07 

.36000E.00  .360OUE.0O 

.OOOOOE.OO 

.40000E.01  .200U0E.01 

•lOOOOE.Ol  .OOOUOE.OO 

.OOOOOE.OO  .OOOOOE.OO 

.UOOOOE.OO  .UOOOUE.OO 

.56000E-05 
.42000£~04 

.58000E.02  .170O0E.00 

.25000E.00  .OOOOOE.OO 

. 10500E.01  .10500E.01 


.04430E.00  .00029E.00 

•OOOOOE.OO  .OOOOOE.OO 

O.OOOOOE.OO 

0.50000E.00  0.50000E.00 

0.30000E.OO 

-.45000E.02  .90000L.02 

•OOOOOE.OO 

.01300E.00  .01300E.02 

.013OOE.OO 

~0.30000E.03  -0.300U0E.03 

-0.30000E.03 

.30000E-01  .480UOE.00 

•lOOOOE.Ol  .120U0E.00 

.31900E.05  .lOOOOE.Ol 

.02000E.00  .O50OOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

O.OOOOOE.OO  O.OOOOOE.OO 


.25000E.00 

.54000t.06 

.36000E.00 

.lOOOOE.Ol 

•OOOOOE.OO 

•lOOOOE.Ol 

.OOOOOE.OO 


.12300E.01 

•OOOOOE.OO 


0.50000E.OO 

.90000E.02 

.01300E.02 

-0.30000E.03 

.38OOOE.00 
•lOOOOE.Ol 
. lOOOOE.Ol 

.04500E.00 


O.OOOOOE.OO 


TABLE  Vm.  -  INPUT  DATA  FOR  THORNE L -40/ EPOXY  COMPOSITE 


TH0RNEL“40/EPOXY 


8  71  54 

1440  1 

.4UOOOE.U8 

.  11000E.07 

.ilOOOE.07 

.20000E.00 

.  15000E.07 

.B0000E.06 

.50000E.06 

.50000E.06 

. 35000E.00 

.OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

.4000aE.01 

.20000E.01 

.40000E.01 

.OOOOOE.OO 

•OOOOOE.OO 

. lOOOOE.Oi 

•OOOOOE.OO 

•OOOOOE.OO 

•UOOOOE.OO 

.lOOOOE.Ol 

. lOOOOE.Ol 

•OOOOOE.OO 

-.55000E-06 

. 56000E-05 

. 56000E-05 

.42800E-04 

.42800E-04 

.42800E-“04 

.58000E.03 

. 58000E.02 

.58000E.02 

.12500E.01 

. 12500E.01 

.25000E .00 

.OOOOOE.OO 

. 22500E.00 

.lOOOOE.Ol 

.31416E.01 

. lOOOOE.Ol 

.  10500E.01 

F 

F 

F 

F 


•OOOOOE.OO 

.05600E.00 

.04430E.00 

.OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0.50000E.OU 

0.50U00E.00 

0.50000E.OO 

0.50000E.00 

0.5000OE.00 

0.50000E.OO 

•OOOOOE.OO 

.45000b. 02 

-.45000b. 02 

-.45000E.02 

.45000E.02 

. OOOOOE.OO 

.OO90OE.OO 

.00900E.00 

.00900E.00 

.0090OE.00 

.00900E.00 

.00900E.00 

-0.30000E.03 

-0.30000E.03 

-0.30000E.03 

-0.30000E.03 

-0.30000E.03 

“0. 30000E.03 

.84000E  .00 

. lOOOOE.Ol 

.08500E.OO 

.16500E.02 

.  lOOOOE.Ol 

.lOOOOE.Ol 

.5U000E .00 

. 13300E.02 

.31900E.05 

•lOOOOE.Ol 

.25000E.06 

. 21000E .05 

.02000E.00 

.04500E.OO 

. 5  000  OE .04 

.OOOOOE.OO 

•OOOOOE.OO 

.10000E.03 

.OOOOOE.OO 

•OOOOOE.OO 

0. OOOOOE.OO 

0. OOOUOE. 00 

0. OOOOOE.OO 

0 .OOOOOE.OO 

•200uut .00 
.150U0t.07 
.350U0t.00 

.200UUE.0i 
•OOOUOE.OO 
•OOOuOE .00 
.UOOuOE.OO 


.17000E.0O 
•OOOUOE .00 

. 10300E.01 


.00027E.00 

•UOOUOE.OO 

0  -bOOuOE.OO 

.900UUt .02 

.009U0E  .00 

-0.300uaE.03 

.460U0E.00 
.OSOOuE.oO 
. lOOouL.Ol 

.0500QE .00 


0 .UOOuuL.OU 


.230O0E.0O 
.30000E.06 
. 35000E.OO 

.OOOOOE.OO 
.OOOOOE.OO 
. lOOOOE.Ol 
•OOOOOE.OO 


.12500E.01 

.OOOOOE.OO 


•OOOOOE.OO 

0.30000E.OO 

.90000E.02 

.OO90OE.O0 

-0.30000E.03 

.27000b. 00 
. lOOOOE.Oi 
•lOOOOE.Ol 

.O45O0E.OO 

O.OOUOOE.OO 


Lu  U-  U.  li_ 


TABLE  DC.  -  INPUT  DATA  FOR  THORNEL-50/EPOXY  COMPOSITE 


THORNEL-50/EPOXY 


8  71  54 

1420  1 

.50000E.08 

•  10000E.07 

.10000E.07 

.20000E.00 

.250OOE.O0 

•20000E»00 

•13000E*07 

.70000E.06 

.130U0E.07 

.57000E.06 

.57000E,06 

•OOOOOE.OO 

.57000E.06 

•OOOOOE.OO 

.36000E.00 

•OOOOOE.OO 

.36000E.00 

.36000E.00 

.40000E.01 

.20000E.01 

.40000E-0i 

.20000E.01 

.OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•23560E.01 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•lOOOOE.Ol 

•lOOOOE.Ol 

-.55000E-06 

.42800E-04 

•lOOOOE.Ol 

.56000E--05 

.42800E-04 

.OOOOOE.OO 

•56000E~05 

.42800E-04 

•OOOOOE.OO 

•OOOOOE.OO 

.58000E.03 

.58000E-02 

.38000E.02 

.i7000E.OO 

.12300E.01 

.12500E-01 

•OOOOOE.OO 

.  12500E.01 
.22500E.00 

.25000E.00 

•OOOOOE.OO 

•OOOOOE.OO 

•lOOOOE.Ol 

.31416E.01 

•lOOOOE.Ol 

.10500E.01 

.10500E.01 

•OOOOOE.OO 

.05900E.00 

.04430E.OO 

.00026E.00 

•OOOOOE.OO 

•OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0.50000E.00 

0.50000E.OO 

0.50000E.00 

0.50000E.00 

O.5OOOOE.00 

0.50000E.00 

0.50000E.00 

0.30000E.00 

•OOOOOE.OO 

.45000E.02 

-.45000E.02 

.9000UE.02 

.90000E.02 

-.45000E.02 

.45000E.02 

•OOOOOE.OO 

.00805E.00 

.00805E.00 

.00800E.00 

•ooaooE.oo 

.00800E.00 

.00805E.00 

.00805E.00 

.0O8OOE.OO 

-0.3U000E.O3 

~0.30000E.03 

-0.30000E.03 

~0.30000E.03 

-0.30000E.03 

-0.30000E.03 

-0.30000E.03 

-0.30000E.03 

.830OOE.OO 

•lOOOOE.Ol 

.26000E.00 

.27000E.00 

.17000E.00 

.16500E.02 

•lOOOOE.Ol 

•lOOOOE.Ol 

•  04650E  *00 

. lOOOOE.Ol 

.SOOOOE.OO 

.lOOOOE.Ol 

.13300E.02 

.31900E.05 

•lOOOOE.Ol 

.lOOOOE.Ol 

.23000E.06 

.04500E.00 

.21000E.05 

.02000E.OO 

.050OOE.00 

.045O0E.00 

.50000E.04 

.OOOOOE.OO 

.OOOOOE.OO 

•10000E.03 

•OOOOOE.OO 

•OOOOOE.OO 

o  o 

.  . 

o  o 
o  o 
o  o 
o  o 
o  o 
m  m 
•  • 
o  o 
o  o 

0. OOOOOE.OO 

0. OOOOOE.OO 

O.OOOOOE.OO 

O.OOOOOE.OO 

TABLE  X.  -  INPUT  DATA  FOR  MODMOR-l/ EPOXY  COMPOSITE 

MOOMOR-I /EPOXY 

8  71  5410000  1 


•60000E.08 

•09000E.07 

.09000E.07 

.200OOE.00 

.25000E.OO 

.20000E«00 

.  11000E.07 

.70000E«06 

.11000E.07 

.50000E.06 

♦50000E.06 

o50000E.06 

.35000E.00 

.35000E.00 

.35000E.00 

•OUOOOE.OO 

. OOOOOE.OO 

.OOOOOE.OO 

.40000E.U1 

.20000E.01 

.40000E.01 

.200O0E ,01 

.OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

. lOOOOE.Ol 

•OOOuOE.OO 

.OOOOOE.OO 

.OOOOuE.UO 

•OOOOOE.OO 

.uOOOOE.OO 

•OOOOOE . 00 

.  lOOOOE.Ol 

•lOOOQE.Ol 

.  lOOOOE.Ol 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

-•55000E-06 

.56000E-05 

.56000E“05 

•42800E-04 

.42800E-04 

.42800E-04 

•58000E.03 

.58000E.02 

.  58000E.02 

.  170OUE.00 

.  12500E.01 

.12500E.01 

.  12500E.01 

.25000E .00 

•UOOowE .00 

.OOOOOE.OO 

•OOOOOE.OO 

.22500E.OO 

. lUUQUE .01 

•lOOOOE.Ol 

.  10500E.01 

.  1  O^ooii  .  01 

.31416E.01 

F 

F 

F 


.OOOOOE.OO 

.OOOOOE.OO 

.07200Eo00 

.OOOOOE.OO 

.04430E.0U 

•OOOOOE.OO 

.00030E .00 
•  aOOuot .  00 

.OOOOOE.OO 

0 .OOOOOE . 00 
0.50000E  .00 

0. OOOOOE. 00 

0 . 50000E .00 

O.OOOOOE .00 

0. 30000E .00 

0  .  bOOooE . 00 

0.50000E.00 

0.50000E  .00 
•OOOOOE.OO 

0.50000E .00 
.45000E.02 

O.50O00E.00 

-.45000E.02 

.90000E.02 

.90000E.02 

~.45000E .02 
.01i90E.O0 

.45000E.02 
.01  190E.00 

.OOOOOE.OO 
.01190t .00 

.01 19UE.00 

.0li90E.O0 

.01190E.OO 
-0.30U00E  .03 

.01190E.00 

-0.30000E.03 

.01190E.00 

-0.30000Eo03 

-0.300uOL.03 

-0.30000t.03 

-0.30000E.03 

.lOOOOE.Ol 

-0.30000E.03 

.lOOOOE.Ol 

-0.30000E.03 
. 5U000E.0O 

.9  lOoi^L  .  00 

.50000E.OO 

.16500E.02 

.lOOOOE.Ol 

.lOOOOE.Ol 

.092uoE  .00 

.lOOOOE.Ol 

.65000E .00 

.13300E.02 

.  31900E.05 

. lOOuOE  .01 

.lOOOOE.Ol 

. lOOOOE.Ol 
.25000E.06 

.21000Ec05 

.02000E .00 

.050uoE .00 

.04500E.OO 

.04500E  .00 
.50000E  .04 
. lOOOoE.03 

0  .OoOOuE .00 

.UOOOOE.OO 

.OOOOOE.OO 

0. OOOOOE. 00 

.OOOOOE.OO 

.OOOOOE.OO 

0 . OOOOOE .00 

O.OQOuuE.OO 

0  .OOOOOE .00 

O.OOOOOE.OO 


TABLE  XI.  -  INPUT  DATA  FOR  MODMOR-H/EPOXY  COMPOSITE 
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.38000E.08 

.11000E.07 

.11000E.07 

.20000E.OO 

.25000E.00 

.20000E.U0 

.15000E.07 

.80000E.06 

.I50o0t.07 

.50000E.06 

•50000E.U6 

.50000E.06 

.35000E.00 

.35OU0E.00 

.35000E.00 

•OOOOOE.UO 

.4U000E.01 

•OOOOOE.OO 

.20000E.01 

.OOOOOE.OO 

.40000E.01 

.200UOE.01 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

. lOOOOE.Ol 

.OOOOOE.OO 

.OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

.lOOOOE.Ol 

.lOOOOE.Ol 

•lOOOOE.Ol 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

-.55000E-06 

.42800E-04 

.58000E.03 

.56000E-05 

.42800E-04 

.58000E.02 

.56000E-05 
.42800E-04 
• 58000E.02 

.  170O0E.00 

.12500E.01 

.12500E.01 

.  12500E.01 

.25000E.00 

•OOOuOE.OO 

•OOOOOE.OO 

.OOOOOE.OO 

•lOOOOE.Ol 

.22500E.00 

•lOOOOE.Ol 

.10500E.01 

.10500E.01 

.31416E.01 

.OOOOOE.OO 

.06300E.00 

.O443OE.O0 

.00030E.00 

F 

F 

F 


•OOOOOE.OO 

.OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

•OOOOOE.OO 

o  o 

.  . 

o  c 
o  o 
o  o 
c  o 
m  m 
•  • 

c  o 
o  o 

0.5OOOOE.0O 

0.50000E.OO 

0.30000E.00 

0.5O000E.00 

0 .30000E.00 

0.50000E.00 

•OOOOOE.OO 

-*.45000E.02 

.45000E.02 

.45000E.02 

-.45000E.02 

•OOOOOE.OO 

.900U0E.02 

.90000E.02 

.01190E.00 

.O119OE.00 

.01190E.00 

.01190E.00 

.01190E.00 

.01190E.00 

.0119OE.00 

.01190E.00 

-0.3UOOOE.U3 

-0.30000E.03 

-0.30000E.03 
-0. 30000E.03 

-0.30000E.03 

“0.30000E.03 

-0.300UOE.03 

-0.30000E.03 

.84000E.00 

.lOOOOE.Ol 

.70000E.00 

.13700E.01 

.80000E.00 

.16500E.02 

•lOOOOE.Ol 

•lOOOOE.Ol 

.160O0E.00 

.10000E..01 

.70000E.OO 

•lOOOOE.Ol 

.13300E.02 

.31900E.05 

.lOOOOE.Ol 

.lOOOOE.Ol 

.35000E.06 

.04500E.00 

.50000E-04 

.10000E.03 

.21000E.05 

•OOOOOE.OO 

•OOOOOE.OO 

.02000E.00 

•OOOOOE.OO 

•OOOOOE.OO 

.05OUUE.00 

.O45O0E.O0 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

0. OOOOOE.OO 

TABLE  Xn.  -  INPUT  DATA  FOR  BERYLLIUM/EPOXY  COMPOSITE 

BERYLLIUM/EPOXY 
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.44000E.08 
•lOOOOE.OO 
.52000E,06 
•UOOOOE.UO 
.40000E.U1 
.OOOOUE.OO 
•OOOOOE.OO 
. lUOOUE.Ol 
.64000E-05 
.3200UE-U4 
• 10440E.04 
.17000E.01 
.OOOOOE.OO 
.lOOOOE.Ol 
.31416E.01 
F 
F 
F 
F 

.OOOOOE.OO 
.OOOOOE.OO 
0. OOOOOE.OO 
0.50000E.U0 
0.50000E.OO 
.OOOOOE.OO 
-.45000E.02 
.OO50OE.O0 
.0O500E.OO 
-0.30000E.03 
-0 .30000E.03 
.lOOOOE.Ol 
. 16500E.02 
.90000E.00 
.lOOOOE.Ol 
.13300E.06 
.05300E.00 
.50000E.04 
. 10000E.03 
0. OOOOOE.OO 
0. OOOOOE.OO 


1  1 
.44000E.08 
.OOOOOE.OO 
.52000E.06 
. OOOOOE.OO 
.20000E.01 
.OOOOOE.OO 
.OOOOOE.OO 
.lOOOOE.Ol 
.64000E-05 
.32000E-04 
.10440E.04 
. 17000E.01 
,22500E.00 
.  lOOOOE.Ol 


.05700E.00 
.OOOOOE.OO 
0. OOOOOE.OO 
0.50000E.00 
0.50000E.OO 
.45000E.02 
.45000E.02 
.0O5O0E.00 
.00500E.OO 
-0.30000E.03 
-0.30000E.03 
.lOOOOE.Ol 
. lOOOOE.Ol 
. 13300E.02 

.25000E.05 

.OOOOOE.OO 
.OOOOOE.OO 
0. OOOOOE.OO 


.44000E.08 
.OOOOOE.OO 
. 35000E .00 
. OOOOOE.OO 
.40000E .01 
. lOOOOE.Ol 
.OOOOOE.OO 
.OOOOOE.OO 
.64000E-05 
.  32000E-04 
. 10440t .04 
.25000E.00 

. lOOOOE.Ol 


.04400E.00 
. OOOOOE.OO 
0. OOOOOE.OO 
0.80000E.OO 
0.5OOOOE.OO 
-.45000E.02 
.OOOOOE.OO 
.00500E.0C 
, 00500E.00 
-0.30000E.03 
-0. 30000E .03 
. 53000E.OO 
. lOOOOE.Ol 
. 31900E.05 

.02700E.00 

.OOOOOE.OO 
.OOOOOE.OO 
0. OOOOOE.OO 


. lOOOOE.OO 
.OOOOOE.OO 
.3500QL.00 

. 20UU0E.01 
.OOOuOE .00 
.OOOuoE.OO 
.jQoooE.OO 


.450OOE .00 
•OOOuOE.OO 

.lOOOOE.Ol 


.00300E .00 
.OOOOOE.OO 

0 .500O0E.00 

.90000E.02 

.OOboOt.OO 

-0 . 3COOOE.03 

. 14000E .01 
.05200E.00 
.lOOOOE.Ol 

.070OOE .00 


0. OOOOOE.OO 


.lOOOOE.OO 
.52000E.06 
.  35000E.00 

.OOOOOE.OO 
•OOOOOE.OO 
. lOOOOE.Ol 
.OOOOOE.OO 


.17000E.01 

.OOOOOE.OO 


•OOOOOE.OO 

O.5OO00E.00 

.90000E.02 

.0C300E.00 

-0. 3CCOOE.03 

. 10500E.G1 
.lOOOOE.Ol 
. lOOOOE.Ol 

.053OOE.O0 


0. OOOOOE.OO 
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TABLE  Xm.  -  INPUT -OUTPUT  FORMAT  IDENTIFICATION  FOR  INPUT  DATA 


Card  first  entry 

Format  statement  number 
(see  compiled  listing) 

Comments 

Output 

heading 

Read 

Write 

THORNEL  50/EPOXY 

— 

4 

4 

Composite  system 

NL 

11 

5 

10 

Integers 

EFll 

70 

35 

37 

Fiber  and  matrix  elastic  constants 

EMU 

are  read  in  one  statement 

VCF 

41 

35 

37 

Correlation  coefficients  for  thermoelastic 

VAF 

40 

35 

37 

Fiber  thermal  coefficients  of  expansion 

VAM 

45 

35 

37 

Matrix  thermal  coefficients  of  expansion 

CHK 

55 

35 

37 

Constituent  heat  conductivities  and  capacities 

BTA 

60 

35 

37 

Correlation  coefficients  for  conductivities 

PIE 

65 

35 

37 

Constant  n 

TLINP 

80 

75 

75 

Boolean  for  thickness 

CSANB 

85 

75 

75 

Boolean  for  bending  symmetry 

BIDE 

87 

75 

75 

Boolean  for  interply  layer  effects 

RINDV 

88 

75 

75 

Boolean  for  load  conditions 

THCS 

90 

35 

37 

Load  angle,  densities,  equivalent  fiber  diameter 

KVL 

95 

35 

37 

Ply  void  content 

KFL 

100 

35 

37 

Ply  fiber  content 

THCL 

105 

35 

37 

Ply  orientation  angle 

TL 

110 

35 

37 

Ply  thickness 

PTEMP 

111 

35 

37 

Ply  temperature 

BET 

115 

35 

37 

Adjustment  factors  for  limit  conditions 

LSC 

120 

35 

37 

Limit  conditions  -  stress,  strain 

NBS 

130 

35 

37 

Load  conditions  -  membrane  forces 

MBS 

131 

35 

37 

Load  conditions  -  bending  moments 

DISVl 

132 

35 

37 

Displacements 

111 


Figure  1.  -  Typical  multilayered  fiber  composite  and  some  basic  definitions. 


114 


Modmor-II/epoxy 


igure  4,  -  Composite  systems  for  which  input  data  are  supplied. 


j  =  KDNj 


Boolean,  real,  array,  and 
string  declarations 


Array  dimensions  and 
integer  declarations 

- H 

Read  and  write 
\  array  dimensions  / 

Read  and  write 
\  input  data  / 


Kvtj- 

Pai 

Pi2j 

‘r' 

’aj 

9cs- 

Pa2j 

orF 

’a3j 

Pn4j 

'Pa2j 

^Pa3f 

Pa3j 

=  Pa3j 

Trim 

Pn4j 

Trim 

End  j  -  LOOP 

i  , 

Call  GECL 

. r 

J  =  2(1)N;  I 


If  ply  properties  are  supplied, 
read  them  in  at  this  location 


|Call  GACD3| 
Call  GPCFD2 
>1  m  =  1(1)N^ 

, . "T 

Call  COMPSA 

zzrz: 

Write  load 
condition 


Write  ply 
properties 


Pjgj  =  0.0186  12 

(6j  +  6i.i)2 

0.0093 

(6  stored  in  P^) 


End  ]-  LOOP 


AdditionaI\No_ 
data  / 


Read  and  write 
additional  input  data 

- ! - - 

_ j _ 

j  Subroutine  INVA 
I 

Subroutine  GLLSC 

! 

Subroutine  GACD3 

I 

Block  data 

I _ 

Subroutine 

- 1 - 

_ I _ 

Subroutine  GPH 
_ \ _ 

Subroutine  GECL 

- ! - 

_ I _ 

Subroutine  GSMF 
Subroutine  CO/VIPSA 


Figure  6.  -  Code  MAIN  PROGRAM  flow  chart 
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